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Abstract. We study microlocal limits of plane waves on noncompact Riemannian mani- 
folds (A/, g) which are either Euclidean or asymptotically hyperbolic with curvature —1 near 
infinity. The plane waves E{z, ^) are functions on M parametrized by the square root of 
energy z and the direction of the wave, ^, interpreted as a point at infinity. If the trapped 
set K for the geodesic flow has Liouville measure zero, we show that, as 2 — > +00, E{z,^) 
microlocally converges to a measure fi^, in average on energy intervals of flxed size, [z, z + 1], 
and in ^. We express the rate of convergence to the limit in terms of the classical escape 
rate of the geodesic flow and its maximal expansion rate — when the flow is Axiom A on 
the trapped set, this yields a negative power of z. As an application, we obtain Weyl type 
asymptotic expansions for local traces of spectral projectors with a remainder controlled in 
terms of the classical escape rate. 



For a compact Riemannian manifold (M, g) of dimension d whose geodesic flow is ergodic 
with respect to the Liouville measure fiL, quantum ergodicity (QE) of eigenfunctions [Sh, 
Ze87, CdV] states that any orthonormal basis (ej)jgN of eigenfunctions of the Laplacian 
with eigenvalues zj, has a density one subsequence (cj^) that converges microlocally to 
in the following sense: for each symbol a S C°° (T* M) of order zero, 

(QPfe,. (a)ejfc,ej,)^2(M) ^ .1 . [ adfii- (1-1) 

^^L{!^*M) Js*M 

Here S*M stands for the unit cotangent bundle, Op;j(a) denotes the pseudodifferential op- 
erator obtained by quantizing a (see Section 3.1), and we put hj = zj^. The proof uses the 
following integrated form of quantum ergodicity [HeMaRo] : 



h'^ ^ (Opfe(a)ej,ej)^2 - ^ [ adfiL 

'1 '-<Zj<h 



asft^O. (1.2) 



See Appendix D for a short self-contained proof of this result using the methods of this paper. 

In the present paper, we consider a non-compact complete Riemannian manifold (M, g) 
and show that generalized eigenfunctions of the Laplacian on M known in scattering theory 
as distorted plane waves or Eisenstein functions, converge microlocally on average, similarly 
to (1.2), with the limiting measure //^ depending on the direction of the plane wave ^ - see 
Theorem 1 . We also give estimates on the rate of convergence in terms of classical quantities 
defined from the geodesic flow on M - see Theorem 2. 

Our microlocal convergence of plane waves is similar in spirit to the QE results (1.1) 
and (1.2). However, unlike the case of QE where ergodicity of the geodesic flow is essential, 
our result is based on a different phenomenon, roughly described as dispersion of plane waves. 
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This difference manifests itself in the proofs as foUows: instead of averaging an observable 
along the geodesic flow as in the standard proof of quantum ergodicity, we propagate it. See 
Section 2 for an outline of the proofs of Theorems 1 and 2. 

Geometric assumptions near infinity. The manifold M has dimension d = n + 1. For 
our results to hold, we need to make several assumptions on the geometry of (M, g) near 
infinity and on the spectral decomposition of its Laplacian A. They are listed in Section 4 
and we check in Sections 6 and 7 that they are satisfied in each of the following two cases: 

(1) there exists a compact set Kq C M such that (M \ go,Ko) is isometric to M"^"*^ \ 
-6(0, Rq) with the Euclidean metric for some Rq > 0; here -6(0, Rq) denotes the ball 
centered at of radius Rq, 

(2) (M, g) is an asymptotically hyperbolic manifold in the sense that it admits a smooth 
compactification M and there exists a smooth boundary defining function x such 
that in a collar neighbourhood of the boundary dM, the metric has the form 

dx"^ + h{x) 

9 = 2 • (l-^) 

where h{x) is a smooth 1-parameter family of metrics on dM for x £ [0,e). We 
further assume that g has sectional curvature —1 in a neighbourhood of dM. 

In case (1), we call {M,g) Euclidean near infinity, while in case (2), we call it hyperbolic near 
infinity. Case (2) in particular includes convex co-compact hyperbolic quotients r\B["+^ - 
see Appendix A. Other possible geometries are discussed in Section 2.1. 

Distorted plane waves/Eisenstein functions. Let A be the (nonnegative) Laplace- 
Beltrami operator on M. In the study of the relation between classical dynamics and high 
energy behavior it is natural to use the semiclassically rescaled operator h'^A, with h > 
small parameter tending to zero. 

The operator h'^A has continuous spectrum on a half- line [cq/i^, oo) (here cq is for the 
Euclidean and for the hyperbolic case), parametrized by distorted plane waves (or 

Eisenstein functions in the hyperbolic case) Eh{X,^) G C°°(M), satisfying for A G M, 

{h^A-X^ -coh^)Eh{X,0 = 0. (1.4) 

Here ^ lies on the boundary dM of a compactification M of M. We can think of an element 
of dM as the direction of escape to infinity for a non-trapped geodesic; then ^ is the direction 
of the outgoing part of the plane wave Eh{X,(^) at infinity. 

For instance, in the case of manifolds Euclidean near infinity, cq = 0, dM = §" is the 
sphere, and for m near infinity, 

where Sine is incoming in the sense that it satisfies a Sommerfeld radiation condition, or 
equivalently, that it lies in the image of C^(M"'^^) under the free (incoming) resolvent 
Ro{X/h) of the Laplacian on the Euclidean space M"+^. These conditions provide a unique 
characterization of Eh{X,^). We can also write £'/i(A,^) = E{X/h,S,), where E{z,S) is the 
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nonsemiclassical plane wave, and rewrite the results below in terms of the parameter z, as 
in the abstract. 

We will freely use the notions of semiclassical analysis as found for example in [Zw], and 
reviewed in Section 3. We denote elements of the cotangent bundle T*M by (m, z^), where 
m ^ M and u G T^M. The semiclassical principal symbol of /i^A is equal to p{m, v) = 
where \i'\g is the length of G T^M with respect to the metric g. Therefore, the plane wave 
Efi should be concentrated on the unit cotangent bundle (see [Zw, Theorem 5.3]) 

S*M := {{m,u) G T*M \ \u\g = 1}. 

If (7* : T*M ^ T*M denotes the geodesic flow, then the Hamiltonian flow of p is e*^^ = (7^*. 

Semiclassical limits of when the trapped set has measure zero. In scattering 
theory trajectories which never escape to infinity play a special role as they can be observed 
only indirectly in asymptotics of plane waves. The incoming tail (resp. outgoing tail) r_ C 
S*M (resp. r+ C S*AI) of the flow is defined as follows: a point (m, ly) lies in r_ (resp. 
r+) if and only if the geodesic g^{'m,u) stays in some compact set for t > (resp. t < 0). 
The trapped set K := r+ n r_ is the set of points (m, i^) such that the geodesic g^{m, v) lies 
entirely in some compact subset of S*M. 

Our first result states that if /Ul(K) = 0, then plane waves Eh{X,^) converge on average 
to some measures supported on the closure of the set of trajectories converging to in M: 

Theorem 1. Let {M,g) be a Riemannian manifold satisfying the assumptions of Section 4 
and suppose that the trapped set has Liouville measure iil{K) = 0. For Lebesgue almost every 
^ G dM, there exists a Radon measure fi^ on S*M such that for each compactly supported 
h-semiclassical pseudodifferential operator A G 'I'^(M), we have as /i — ?■ 0, 



{AEhiX, 0,Eh{>^, O)l^(m) - [ <^(^) dfi^ 

JS*M 



_ ^0. (1.5) 



The measure has support 



supp(/if ) C {(m, z/) G 5*M I lim g^(m,u) = (^}, (1-6) 

t— s>+oo 

and disintegrates the Liouville measure in the sense that there exists a smooth measure d^ on 
dM such that, if is the Liouville measure generated by ^/p = \v\g on S*M, then 



Ifi^dC = fJ.L- (1-7) 
JdM 



IdM 

The limiting measure is defined in Section 4.3. Implicit in (1.7) is the statement that 
for any bounded Borel S C S*M, we have n^{S) G L^{dM). 

In the case when WF fi(A) n r_ = (in particular when g is non-trapping), we actually 
have a full expansion of {AEh,Efi) in powers of h, with remainders bounded in L^^{dM x 
[l,l + /i]) -see (5.14). 

The now standard argument of Colin de Verdiere and Zelditch (see for example the proof 
of [Zw, Theorem 15.5]) shows that there exists a family of Borel sets A{h) C dM x [1, 1 + h] 
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such that the ratio of the measure of A{h) to the measure of the whole dM x [1, 1 + /i] goes 
to 1 as /i —7- 0, and for each A £ 'if^{M) as in Theorem 1 with (7{A) independent of h, 

{AEh{X,0,Eh{X,0)LHM) ^ / c7{A)dfi^ uniformly in (A,^ S A{h). (1.8) 

Js*M 

This statement can be viewed as an analogue of the quantum ergodicity fact (1.1), though 
as explained above, it is produced by a different phenomenon. 

Estimates for the remainder. We next provide a quantitative version of Theorem 1, 
namely an estimate of the left-hand side of (1.5). We define the set T{t) of geodesies trapped 
for time t > as follows: let Kq be a compact geodesically convex subset of M containing a 
neighborhood of the trapped set K, then (see also Section 5.2) 

r{t) := {(m, ly) G S*M \ m e Kq, g\m, v) G Kq}. (1.9) 

A quantity which will appear frequently with some parameter A > is the following inter- 
polated measure 

r(/i,A):= sup h^-^ ^lL{T[^K-^\\ogh\)) , (1.10) 

0<6»<1 

where /i > is small. This converges to as /i — )• when iii{K) = and it interpolates 
between h (when ^ = 0) and the Liouville measure of the set of geodesies that remain trapped 
for time A~^| log/i| (when 9 = 1). When the measure fJiLiTif)) decays exponentially in t, as 
in (1.14), r{h,A) can be replaced by simply 0{h) + ^l{T{A~^\ logh\)). The 0{h) term here 
is natural because of the influence of the subprincipal part of the operator A. 

We next define the maximal expansion rate as follows (see also (3.17)): 

Amax := limsupy^log sup \\dg\m,i/)\\. (l-H) 

We can estimate the left-hand side of (1.5) in terms of the (interpolated) measure of the set 
of all trajectories trapped for the Ehrenfest time. If we pair with a test function in ^ instead 
of taking the L| norm, then the estimate becomes stronger, corresponding to the set of all 
trajectories trapped for twice the Ehrenfest time: 

Theorem 2. Let {M,g) be as in Theorem 1. Take Aq > Amax- Then for each compactly 
supported h-semiclassical pseudodijjerential operator A £ '^^(M) and for each f £ C°°{dM), 



{AEh,Eh) - [ a{A)dfi^ 

J S*M 

[_ fiO ({AEh, Eh) - I a{A) d^,^ di 

JdM \ JS*M J 



= 0(r(/i,2Ao)), (1.12) 



= 0{r{h,Ko))- (1.13) 

Li([l,l+h,]) 



The proof of Theorem 2 actually gives an expansion of {AEh,Eh) in powers of h, with re- 
mainder ij,l{T{A~^\ log/i|)) instead of r(/i, A) - see (5.35) and the proofs of Propositions 5.11 
and 5.13. This full expansion is cumbersome to write down, therefore we only do it for the 
trace estimates (1.16) below. 
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Remainder in terms of pressure. When the trapped set K has Liouville measure and 
is uniformly partially hyperbolic in the sense of Appendix B.l, we estimate using [Yo] 

piL{T{t)) = 0(e*(^(-^")+^)), (1.14) 

for each e > 0, where P{J'^) < is the topological pressure of the unstable Jacobian - see 
Appendix B.l. When K is a hyperbolic basic set (Axiom A flow), then P(J") < by [BoRu], 
and the remainders in (1.12) and (1-13) are then polynomial in h: 

r(/i,A) = 0(/i + /i-(^(^")+^)/^), 

and one can get rid of e here by slightly changing Aq . In the special case where g has constant 
sectional curvature —1 near K, the bounds in (1.12) and (1.13) become 
and 0{h + /i"~^~), respectively, where K has Hausdorff dimension dim/f(i^) = 25 + 1. See 
Appendix B.2 for details. 

In all cases, if K is nonempty, then it has Minkowski dimension at least 1; since g^^'^{T{t)) 
contains an e"^"*/^ sized neighborhood of K., we have 

/iL(r((2Ao)-i|iog/i|)) > /^"/^ ^JiL[T{^.-,^\\o^h\))>h^. (1.15) 

Local Weyl asymptotics for spectral projectors. It is possible to express the spectral 
measure of /i^A in terms of the distorted plane waves (see (4.5)), and using (1.13), we obtain 
an expansion in powers of h for local traces of spectral projectors up to an explicit remainder. 
We only write it here for the case where the flow is partially uniformly hyperbolic with 
P(J") < 0, but a more general result with the Liouville measure of r(Aoi|log/i|) holds - 
see Theorem 4 in Section 5.3. Below, we fix a quantization procedure Op/j on M mapping 
compactly supported symbols to compactly supported operators. 

Theorem 3. Let (M, g) be as in Theorem 1, let Aq > Amax o.'ri'd assume that the trapped set 
K is uniformly partially hyperbolic with iii{K) = and that the topological pressure P{J^) 
of the unstable Jacobian on K is negative. Then there exist differential operators^ Lj of 
order 2j on T*M , with Lq = 1, such that for each compactly supported zeroth order classical 
symbol a, we have for each s > and N 

T¥(Op^(a)l[o,,](/i^A)) = (27r/i)-"-i^/i-'' J Ljodfi^ + h-''0{h~^ + h^) (1.16) 

J=o \u\l<s 

where is the standard volume form on T*M and l[o^s](^^^) denotes the spectral projector 
of h^A onto the frequency window [0,s]. The remainder is uniform in s when s varies in a 
compact .subset of (0, oo). 

In particular, if g has constant sectional curvature —1 near K and the Hausdorff dimension 
of K is given by 26 + 1, then the remainder in (1.16) becomes 0{h~^~), for N large enough. 

Applications. In a separate paper [DyGu], we show that Theorem 3 implies new asymp- 
totics for the spectral shift function (or scattering phase) with remainders in terms of P( J") 



In this paper, the symbols Lj will denote different operators in different propositions. 
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when the trapped set has Liouville measure and the manifold is EucUdean near infinity 
with uniformly partially hyperbolic geodesic flow near K. 

Previous works. Let us briefly discuss the history of Quantum Ergodicity (QE) and explain 
its relation to the present paper. The original QE statement was proved by Shnirelman [Sh], 
Zelditch [Ze87], and Colin de Verdiere [CdV] in the microlocal case, by Helffer-Martinez- 
Robert [HcMaRo] in the semiclassical case (with the integrated estimate using an 0{h) 
spectral window like in the present paper, rather than the 0{1) window used in the microlocal 
case), and by Gerard-Leichtnam [GeLe] and Zelditch-Zworski [ZcZw] for manifolds with 
boundary (ergodic billiards). Quantum ergodicity for boundary values and restrictions of 
eigenfunctions to hypersurfaces was studied by Hassell-Zelditch [HaZc], Burq [Bu05], Toth- 
Zelditch [ToZclO, ToZcll], and by Dyatlov-Zworski [DyZw]. 

The first result on noncompact manifolds, namely for embedded eigenvalues and Eisenstein 
functions on surfaces with cusps, was proved by Zelditch [Zc91]. For the special case of 
arithmetic hyperbolic surfaces, a stronger statement of Quantum Unique Ergodicity (QUE), 
saying that the whole sequence of eigenstates microlocally converges to the Liouville measure, 
was proved by Lindenstrauss [Li] and Soundararajan [So] for Hecke-Maass forms and by Luo- 
Sarnak [LuSa] and Jakobson [Ja] for Eisenstein functions. For further information on the 
topic, the reader is directed to the recent reviews [No, Sa, Ze09]. 

As remarked above, our result differs from the above works in that it uses dispersion of 
plane waves instead of the ergodicity of the geodesic flow. This dispersion phenomenon was 
used to study microlocal limits of plane waves on convex co-compact hyperbolic quotients 
satisfying 5 < n/2 by Guillarmou-Naud in [GuNa], and on surfaces with cusps at complex 
energies by Dyatlov [Dy2]. Both [GuNa] and [Dy2] guarantee microlocal convergence of the 
Eisenstein functions that is uniform in A and ^, rather than the (weaker) L\ ^ estimates 
of the current paper; these statements are formally similar to QUE, while our statement 
is formally similar to QE. In [GuNa], uniform in A and ^ estimates are possible because 
Lagrangian states, when propagated by the Schrodinger group U{t), would disperse faster 
than they fail to be approximated semiclassically, a phenomenon similar to the one studied 
by Nonnenmacher-Zworski [NoZw]. In fact, it is plausible that the result of [GuNa] is true 
when the condition 5 < n/2 is replaced by the negative pressure condition of [NoZw]. As 
for [Dy2], the energy being away from the real line makes the measure corresponding to Eh 
exponentially increasing, rather than invariant, along the flow; then the result of propagation 
of by U{t) is multiplied by e"*^* for a certain v > 0, and decays in as t — >• +oo (it 
is then more correct to say that this paper relies on damping of plane waves rather than 
dispersion). 

We see that the uniform convergence in [GuNa] and [Dy2] is possible because one has 
better control on the propagated Lagrangian states. Such better control is directly related 
to having a polynomial bound on the scattering resolvent. In the less restricted situation of 
our paper, however, it is not clear if such a bound would hold; therefore, we need to average 
in A and ^ to pass to trace (or, strictly speaking, Hilbert-Schmidt norm) estimates, just as 
in the proof of Quantum Ergodicity. 
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The expansions for local traces of the spectral measure as in Theorem 3 were studied by 
Robert-Tamura [RoTa] for nontrapping perturbations of the Euclidean space, yielding a full 
expansion in powers of h in that setting. 

2. Outline of the proofs 

In this section, we explain the ideas of the proofs of Theorems 1 and 2, in the case of 
manifolds Euclidean near infinity. We also describe the structure of the paper. 

We start with Theorem 1. Take t > 0; we will use limf_>.+oo lim/i-!>o limits, therefore 
remainders that decay in h with constants depending on t will be negligible. Since Eh is a 
generalized eigenfunction of the Laplacian (1.4), we have 

Eh{X,0 = e-''''"/^'''^U{t)Eh{X,0. (2.1) 

Here U{t) = g**'^^/^ is the semiclassical Schrodinger propagator, quantizing the geodesic flow 
g^. Since Eh does not lie in L'^{M), we cannot apply the operator U{t) to it; however, (2.1) 
can be made rigorous, with an 0{h°°) error, by using appropriate cutoffs - see Lemma 3.10. 
We will not write these cutoffs here for the sake of brevity. 

Take a compactly supported and compactly microlocalized semiclassical pseudodifferential 
operator A on M; then by (2.1), 

{AEh,Eh) = {AU{t)Eh,U{t)Eh) = {A-' Eh, Eh), (2.2) 

where A~* := U{—t)AU{t) is a pseudodifferential operator with principal symbol (t{A) o g^^. 
(It is not compactly supported, but we ignore this issue here.) We now use the following 
decomposition of plane waves: 

Eh = XoE^h + El El = ex«-, Ei = -Rhi\)Fh, Fh := {h^A - X^)xoEI 

Here E^ is the outgoing part of the plane wave, defined in a certain neighborhood of infinity 
and solving (1.4) there, while xo is a cutoff function equal to 1 near infinity and supported 
inside the domain of E^; then 

Fh = [h^A,xo]El 

is compactly supported and we can apply to it the semiclassical scattering resolvent Rh{X). 
Note that here Rh{X) is the incoming resolvent; in particular, it is bounded — )■ for 
Im A < 0. (The situation is more complicated in the case of manifolds hyperbolic near infinity, 
in particular the domain of E^ and the cutoff xo will depend on ^.) For A = 1 + 0{h), the 
function Fh is microlocalized inside the set 

:= {{m,u)\me supp(dxo), = C S*M. 

In general, we cannot expect the resolvent Rh{X) to be polynomially bounded in h, and 
thus cannot determine the wavefront set of Ej^. However, we will show the following weaker 
propagation of singularities statement: the function 
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Figure 1. A phase space picture of the main argument. The right side of 
each picture represents infinity; xo = 1 in the hghter shaded region and dxo is 
supported in the darker shaded region, while (/9 = 1 to the left of the vertical 
dashed line. The horizontal dashed lines on the right represent the wavefront 
set of E\] they terminate at the solid arrows, which denote the set W^. 



where Kq C M is a sufficiently large compact set, is polynomially bounded in h and for 
each {m,v) G WF/i(S^), the geodesic g^{m,v) is either trapped as t — ?■ +oo or passes 
through for some i > 0. For the case of manifolds Euclidean near infinity, this statement 
follows directly from the explicit formula for the scattering resolvent on the free Euclidean 
space; for manifolds hyperbolic near infinity, we use the microlocal properties of the resolvent 
established in [Vail]. See assumption (A6) in Section 4.2, Section 6.2, and Proposition 7.4 
for details. 

If A and 1 — Xo are both supported in the ball of radius i?, let if G C^{M) be independent 
of t and equal to 1 in the ball of radius i? + 1. Then we write 

A-' = * + A-^\ := A-V, := " 

Now, each {m,v) G WF/i(Aj"*) has the following properties: |m| > -R + 1, and for {m' ,u') = 
g'^{m,i'), \m'\ < R. (See Figure 1.) Therefore, the geodesic g^{m,v) escapes to infinity for 
s > and never passes through ; it follows from the discussion of the wavefront set of 
in the previous paragraph that 

\\A^'El\\L2=0{h^{l + \m\L2i^K,)))- 

Therefore, we can write 

{AEh,Eh) = {A-{\oElx^El) + {A^'Eh,Eh) + 0{h^{l + \\Eh\\l,^j,^^))). (2.3) 

The first term on the right-hand side is explicit, as we have a formula for E^] we can calculate 
for Lebesgue almost every ^ and A = 1 + (D{h), 

Ihn lim{A^\oE'fl{X,0,XoE'dX,0) = [ ad^x^. (2.4) 

It then remains to estimate the second and third terms on average in A and ^. For this, we 
use the relation (4.5) of distorted plane waves to the spectral measure of the Laplacian to 
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get for any bounded compactly supported pseudodifferential operator B, 

^~i^^^(^'Olli^^^^^(MxaMx[i,m]) < C/^'^pl[i,(m)^](/^^A)||^s- (2.5) 

Here HS denotes the Hilbert-Schmidt norm. One can estimate the right-hand side of (2.5) 
uniformly in h - see Lemma 3.11 and the proof of Proposition 4.5. Then h~^\\Eh\\'j^2(^Xo)^ 
when integrated over As [1, 1 + /i] and is bounded uniformly in h; this removes the third 
term on the right-hand side of (2.3). 

Finally, the average in A,^ of the second term on the right-hand side of (2.3) can be 
bounded, modulo an Ot{h) remainder, by the norm ||<7(Ao*)IIl2(5*m) of the restriction 
of the principal symbol of A^^ to the energy surface S*M, with respect to the Liouville 
measure. Now, cj(j4q*) = {(^{A) o g^^)(p converges to zero as t — )• +oo at any point which 
is not trapped in the backwards direction. Since we assumed iil{K) = 0, by the dominated 
convergence theorem ||o"(^o*) 11/^2(5* j^,/) converges to zero as t — )• +c«; this finishes the proof 
of Theorem 1. 

For the estimate (1-12) in Theorem 2, we need to take t up to the Ehrenfest time: 

t = te:= Ao ' log(l//i)/2, 

replacing the limt_>+oo lim/i_j.o limit in the argument of Theorem 1 by just the lim/i_j.o limit, 
but with t depending on h. The operator A^^ is then still pseudodifferential, though in a 
mildly exotic class. To avoid a quantization procedure uniform at infinity, we give an iterative 
argument, propagating A for a fixed time for ~ log(l//i) steps, applying t-independent cutoffs 
and removing the microlocally negligible terms at each step. The proof then works as before, 
with the term {A^^EhjEh) bounded by the Liouville measure of the support of the full 
symbol of Aq*, which depends on h and is contained in g^{T{t)), where T{t) is defined 
in (1.9); this proves (1.12). The interpolated quantity r{h,A) from (1.10) appears because 
of the subprincipal terms in (2.4). 

For (1.13), we have to propagate to twice the Ehrenfest time: t = 2te. The operator A~^ 
is not pseudodifferential, but we can use (2.1) to write 

{A^'Eh, Eh) = {U{-t/2)AU{t/2) ■ U{t/2)ipU{-t/2)Eh, Eh). (2.6) 

The operators U{—t/2)AU{t/2) and U{t/2)ipU{—t/2) are both pseudodifferential in a mildly 
exotic class; multiplying them, we get a pseudodifferential operator whose full symbol is 
supported inside g^^'^{T{t)), and thus (2.6) can be estimated by the Liouville measure of this 
set, giving the remainder (1.13). 

A problem arises when trying to get a rate of convergence in (2.4) for t up to twice the 
Ehrenfest time. We are unable to propagate the Lagrangian state E^{X, ^) pointwise in ^ and 
A for time t, therefore we do not get an estimate in (1.13). However, for / G C°^{dM) 
we can still approximate the integral 

f_f{0{Ai\oElxoEl)d^ (2.7) 

JdM 
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as follows. Define the operator 

n^(A) := / /(e)(X0i?^(A,O) ® {xoE'i{X,0)d^. 

JdM 

Here (8) denotes the Hilbert tensor product; that is, if u,v G C°°{M), then u v is the 
operator with the Schwartz kernel 

Ku(g)v{m,rn) = u{m)v{m'). (2.8) 

We can show that for X G iI''^°™p(M) satisfying certain conditions, XH^^X* is a Fourier 
integral operator associated to the canonical relation 

{(m, V] m , u') I (m, v) G S*M^ {m , v') = g'^{m, v) for some s £ (—To, Tq)}, 

for a fixed Tq > depending on X. (For comparison, for the spectral measure of h'^A 
we would have to formally take all possible values of s, which would destroy any hope on 
microlocally approximating it when the geodesic flow is chaotic.) We can then write 

XUj{X)X* = (2^/i)" e-'^''/^^^^U{s)Bsds, 

J -To 

where Bg is a smooth family of pseudodifferential operators, compactly supported in s G 
(—To, To) - see Lemma 5.12. We then write the integral (2.7) as 

TT{U{-t)AU{t){l - (^)n|^(A)) = TV e~'^''/^'^^^U{-t)AU{t){l - ip)U{s)Bs ds 

J —To 

/To 
g-a2s/(2h)f^(_^/2)AC/(t/2) • C/(t/2)(l - ip)U{s)B,Ui-s - t/2) ■ U{s) ds. 
-To 

The operators U{-t/2)AU{t/2) and U{t/2){1 - ip)U{s)BsU{-s -t/2) are pseudodifferential 
in a mildly exotic class; thus their product is also pseudodifferential and (bearing in mind 
that s varies in a bounded set) one gets a microfocal expansion for (2.7) through a local trace 
formula for Schrodinger propagators - see Lemma 3.12 and Proposition 5.13. 

2.1. Other possible geometric assumptions. Our results should be true for asymptot- 
ically hyperbolic manifolds without the constant curvature assumption near infinity. The 
main difficulty here is constructing a good semiclassical parametrix for the Eisenstein func- 
tion Eh{X,(^) near G dM; this can be done by WKB approximation, and the phase is 
a Busemann function 4>^{m) near ^, however one would need a good understanding of the 
regularity of 4>^{m) as m — )• ^. This is in a way related to the high-frequency parametrix of 
[McSBVa] in the non-trapping setting. For the asymptotically Euclidean of asymptotically 
conic ends, this might be more complicated as we would need a parametrix of Eh{X,S^) in 
a large neighbourhood of ^ G dM, essentially in a region with closure containing a ball of 
radius 7r/2 in dM. In particular, the Lagrangian supporting the semiclassical parametrix of 
EhiX,^,) would not a priori be projectable far from ^, which would make the construction 
more technical. We leave these questions for future research. 

The convergence result in Theorem 1 should be true in the case where M has a boundary, 
for instance M = M"+i \ n with n a piecewise smooth obstacle. In fact, it should be 
straightforward to check that the method of proof applies when combined with the idea 
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of [ZcZw], based on the fact that the region in phase space near the boundary where the 
dynamics is comphcated is of Liouville measure (since we assume = 0). To get a 

good remainder in that setting would be more involved since one would need to care about 
the amount of mass of plane waves staying in the regions near the boundary where the 
dynamics is complicated, as we propagate up to Ehrenfest time. A reasonable case to start 
with is that of strictly convex obstacles. 

2.2. Structure of the paper. In Section 3, we review certain notions of semiclassical anal- 
ysis and derive several technical lemmata; in particular, in Section 3.2, we review the local 
theory of semiclassical Lagrangian distributions and Fourier integral operators and in Sec- 
tion 3.3 we study microlocal properties of Schrodinger propagators, including the Hilbert- 
Schmidt norm bound (Lemma 3.11). In Section 4, we formulate the general assumptions on 
the studied manifolds and derive some immediate corollaries; Section 4.1 contains the geo- 
metric assumptions and the definition of the trapped set and Section 4.2 contains the analytic 
assumptions on distorted plane waves. In Section 4.3 we construct the limiting measures /ng 
and in Section 4.4 we prove averaged estimates on Eisenstein functions. 

In Section 5, we give the proofs of our main theorems. Section 5.1 contains the proof 
of Theorem 1, Section 5.2 contains the proof of the estimate (1-12) in Theorem 2, while 
Section 5.3 contains the proof of the estimate (1-13) in Theorem 2. Section 5.3 also contains 
the Tauberian argument proving an expansion of the local trace of a spectral projector 
(Theorem 4). Sections 6 and 7 study the Euclidean and hyperbolic near infinity manifolds, 
respectively, and show that the general assumptions of Section 4 are satisfied in these cases. 

Appendix A provides a formula for the limiting measures in the case of a convex co-compact 
hyperbolic quotient, which generalizes the limiting measure of [GuNa] to the case 6 > n/2. 
Appendix B discusses the classical escape rate, in particular explaining (1.14). Appendix C 
gives a self-contained proof of Egorov's theorem up to the Ehrenfest time (Proposition 3.9). 
Finally, Appendix D contains a short proof of (a special case of) quantum ergodicity in 
the semiclassical setting, which is simpler than that of [HeMaRo] because it does not rely 
on [DuGu, PcRo]. 

3. Semiclassical preliminaries 

In this section, we review the methods of semiclassical analysis needed for our argument. 
Most of the constructions listed below are standard: pseudodifferential operators, wavefront 
sets, local theory of Fourier integral operators, and Egorov's theorem. However, Section 3.3 
contains the propagation result for generalized eigenfunctions (Lemma 3.10) and a Hilbert- 
Schmidt norm estimate in an 0{h) spectral window (Lemma 3.11), which the authors were 
unable to find in previous literature. 

We will also need Egorov's theorem up to the Ehrenfest time (Proposition 3.9); while 
several versions of this fact are available, we could not find a detailed proof for the case of 
manifolds and when the Ehrenfest time is defined via the maximal expansion rate of the 
flow. For this reason, and also because we insert cutoffs in between the propagators, we give 
a proof of Proposition 3.9 in Appendix C. 
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3.1. Notation. In this subsection, we briefly review certain notation used in semiclassical 
analysis. The reader is referred to [Zw] (especially Chapter 14 on semiclassical calculus on 
manifolds) or [DiSj] for a detailed introduction to the subject. 

The phase space. Let M be a d-dimensional manifold without boundary. We denote points 
in M by the letter m and elements of the cotangent bundle T*M by (m, u), where v G T.^M. 
Following [Vail, Section 2], we consider the fiber-radial compactification T* M of T*M. The 
boundary of T*M, denoted by dT* M and called the fiber infinity (unlike [Vail], we do 
not use the notation S*M for fiber infinity — we reserve it for the unit cotangent bundle 
{{I'lg = 1} C T*M), is associated with the cosphere bundle over M and the interior of T* M 
is associated with T*M. Take some smooth inner product on the fibers of T*M; if |z/| is the 
norm of a covector (m, z/) G T*M generated by this inner product and (u) = y^Y+~\u\^, then 
(z^)~^ is a boundary defining function on T* M . 

We will mostly use compactly microlocalized operators, for which the fiber-radial compact- 
ification is not necessary. However, it will come up in the elliptic estimate (Proposition 3.2) 
and in the proof of the propagation of singularities result for plane waves on asymptotically 
hyperbolic manifolds (Proposition 7.4). 

Symbol classes. For any A; G M and any p G [0, 1/2), we consider the symbol class S^{M) 
defined as follows: a smooth function a(?7i, v] h) on T*M X [0, ho) hes in 5^(M) if and only 
if for each compact K C M and each multiindices a, /3, there exists a constant Ca/3K such 
that for h small enough, 

sup \d^di!a{m,u;h)\ < C„^^/i-''(H+l/?l)(^)fc-l/3|. (3.1) 

meK, ueT^M 

These classes are independent of the choice of coordinates on M. Note that we do not fix the 
behaviour of the symbols as m — >• 00. The important special case is p = 0, which includes 
the classical symbols studied in [Vail]. The class Sq{M), denoted simply by S^{M), would 
be sufficient for the convergence Theorem 1. The classes with p > will be important 
for obtaining the remainder estimate of Theorem 2; these classes arise when propagating 
symbols in for short logarithmic times, as in Proposition 3.9. 

Since plane waves are microlocalized on the cosphere bundle, away from the fiber infinity, 
we will most often work with the classes S'p°™'', consisting of compactly supported functions 
satisfying (3.1); we have 5p°™^ C for all k. 

PseudodifTerential operators. Following [Zw, Section 14.2], we can define the algebra 
^'p(M) of pseudodifferential operators with symbols in S^{M). (The properties of the sym- 
bol classes required for the construction of [Zw, Section 14.2] are derived as in [Zw, 
Section 4.4]; see also [GrSj, Chapter 3].) As before, denote = ^'q. Since our symbols 
can grow arbitrarily fast as m — 00, we do not make any a priori assumptions on the be- 
havior of elements of ^'p near the infinity in M. However, we require that all operators 
A G ^''^(M) be properly supported] namely, the restriction of each of the projection maps 
T^m^T^m' : M X M — )• M to the support of the Schwartz kernel KA{m,m') of ^ is a proper 
map. Then each element of "^^{M) acts Hf^ ^^^{M) Hl~^^^{M), where H^ ^^^iM) denotes 
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the space of distributions locaUy in the semiclassical Sobolev space Hf^ (see for example [Zw, 
Section 7.1] for the definition of semiclassical Sobolev spaces). 

We have the semiclassical principal symbol map 

a : ^^(M) ^ S^p{M)/h^-^PS^p-\M) 

and its right inverse, a non-canonical quantization map 

Op, : S'p{M) ^ M/^(M). 

The standard operations of pseudodifferential calculus with symbols in have an 0{h^~'^P) 
remainder instead of the 0{h) remainder valid for the class Sq. More precisely, we have for 
A £ ^'^(M) and B £ ^'^'(M), 

a{A*)=^) + 0{h'-^P)^^.-^^^^y 

a{AB) = a[A)a{B) + ©(/i^-^O^.+.^.i^^^), 

a{[A,B]) = -ih{a{A),a{B)] + 0{h''^^-^P^)^,^,,^,^^^^^. 

Here {•, •} stands for the Poisson bracket. The 0{-) notation is used in the present paper in 
the following way: we write u = Oz{F)x if the norm of the function, or the operator, u in 
the functional space X is bounded by the expression F times a constant depending on the 
parameter z. 

Wavefront sets. If A : C°°(M) — )• C°°{M) is a properly supported operator, we say 
that A = C'(/i°°)^-<x) if A is smoothing and each of the C°°{M x M) seminorms of its 
Schwartz kernel is 0{h°°). For each A £ "^piM), we have A = Op,(a) + 0{h°^)^-ac for 
some a £ S^{M). Define the semiclassical wavefront set WF/i(A) C T M of j4 as follows: 
a point {m,v) £ T* M does not lie in WFh(A), if there exists a neighborhood U of {m,v) 
in T* M such that each (m, z^)-derivative of a is C'(/i°°(i/)^°°) in [/ n T*M. The notion 
of the wavefront set does not depend on the choice of the quantization procedure. We 
have WFh{A) = ff and only if A = C'(/i°°)vp-oc and WFhiA*) = WFh{A), WFh{AB) C 
WFh{A) n WFh{B). For A, B £ ^'^(M), we say that A = B microlocally in some open set 
U C T*M, if WFhiA -B)r\U = $. 

Operators with compact wavefront sets are called compactly microlocalized; those are 
exactly operators of the form Op,(a) + C'(/i°°),j,-cx) for some a £ Sp°"^'^. We denote by 
\I'p°™^(M) the class of all compactly microlocalized elements of \I'p(M); as before, we put 
^'™™P(M) = 'I'q°™''(M). Compactly microlocalized operators should not be confused with 
compactly supported operators; that is, operators whose Schwartz kernels are compactly sup- 
ported. That being said, most operators that we use will be both compactly supported and 
compactly microlocalized. 

We will need a finer notion of microsupport on /i-dependent sets, used in the proofs in 
Sections 5.2 and 5.3, for example in Proposition 5.9: 

Definition 3.1. An operator A £ ^'p°™^(M) is said to be microsupported on an h-dependent 
family of sets V{h) C T*M , if we can write A = Op,(a)-|-C'(/i°°),j,-oo, where for each compact 
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set K C T*M , each differential operator on T*M, and each N , there exists a constant 
CaN such that for h small enough, 

sup \d'^a{m,i^;h)\ < CaNh^ . 

{rn,u)GK\V{h) 

Since the change of variables formula for the full symbol of a pseudodifferential operator 
contains an asymptotic expansion in powers of h, consisting of derivatives of the original 
symbol, Definition 3.1 does not depend on the choice of the quantization procedure Op^. 
Moreover, if ^ G vl'™™^ is microsupported inside some V{h) and B £ then AB, BA, 
and A* are also microsupported inside V{h). It follows from the definition of the wavefront 
set that (m, i^) G T*M does not lie in WFh{A) for some A S ^'p°™'', if and only if there 
exists an /i-independent neighborhood U of (m, v) such that A is microsupported on the 
complement of U. Note however that A need not be microsupported on WFh{A), though it 
will be microsupported on any /i-independent neighborhood of WFhiA). Finally, it can be 
seen by Taylor's formula that if ^4 G ^'p°™''(M) is microsupported in V{h) and p' > p, then 
A is also microsupported on the set of all points in V{h) which are at least away from 
the complement of V{h). 

Ellipticity. For A G *p(M), define its elliptic set e\l{A) C T*M as follows: {m,u) E 
ell(74) if and only if there exists a neighborhood U of (m, z^) in T*M and a constant C such 
that |cr(A)| > C~^{u)^ in U r\T*M. The following statement is the standard semiclassical 
elliptic estimate; see [H6III, Theorem 18.1.24'] for the closely related microlocal case and for 
example [Dyl, Section 2.2] for the semiclassical case. 

Proposition 3.2. Assume that P £ ^'p(M), A £ ^'^'(M), and WFh{A) C ell(P). Assume 
moreover that A is compactly supported. Then there exists a constant C and a function 
X £ C^{M) such that for each s G M, each u S Hf^ ioc(-^) '^'^'^ each N , we have 

\\Au\\hi ^ C\\xPu\\^.^y.u+0[h'^)\\xu\\H-N. 

h 

Moreover, if P is a differential operator, then we can take any x such that the Schwartz 
kernel of A is supported in {% / 0} x / 0}. 

Semi-classical wave- front sets of distributions. An /i-dependent family u{h) G V'{M) 
is called h-tempered, if for each open U compactly contained in M, there exist constants C 
and such that 

||n(/i)||^-^(^)<C/i-^. (3.2) 

For a tempered distribution u, we say that (mo, i^o) £ T* M does not lie in the wavefront set 
WF/j(n), if there exists a neighborhood V{mQ, vq) in T* M such that for each A G ^^{M) with 
WFhiA) C V, we have Au = 0{h°°)c^. By Proposition 3.2, (mo, fo) WF/j(m) if and only if 
there exists compactly supported A G ^^{M) elliptic at (mo, i^o) such that Au = 0{h°°)c°°- 
The wavefront set of u is a closed subset of T M; it is empty if and only if n = O{h^)(joa(^j^iy 
We can also verify that for u tempered and A G ^'p(M), WFh{Au) C WFh{A) n WF/j(n). 
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3.2. Semiclassical Lagrangian distributions. In this subsection, we review some facts 
from the theory of semiclassical Lagrangian distributions. See [GuSt, Chapter 6] or [VuNg, 
Section 2.3] for a detailed account, and [HoIV, Section 25.1] or [GrSj, Chapter 11] for the 
closely related microlocal case. However, note that we do not attempt to define the prin- 
cipal symbols as global invariant geometric objects; this makes the resulting local theory 
considerably simpler. 

Phase functions. Let M be a manifold without boundary. We denote its dimension by d; 
in the convention used in the present paper, d = n + 1. As before, we denote elements of 
T*M by {m,v), m £ M, v G T^M. Let ip{m,9) be a smooth real-valued function on some 
open subset Uip of M x R-^, for some L; we call m base variables and 9 oscillatory variables. 
We say that cp is a (nondegenerate) phase function, if the diff'erentials d{dg-^ip), . . . , d{dgj^ip) 
are linearly independent on the critical set 

C^:={{m,e)\deip = 0} CU^. (3.3) 

In this case 

A^ := {(m, 5^(^(m, 6)) \ (m, 9) € C^} C T*M 

is an (immersed, and we will shrink the domain of ip to make it embedded) Lagrangian 
submanifold. We say that if generates A^. 

Symbols. Let p G [0,1/2). A smooth function a{m,9;h) is called a compactly supported 
symbol of type p on U^, if it is supported in some compact /i- independent subset of U^,, and 
for each differential operator d"' on M x M^, there exists a constant Ca such that 

sup|9"a| < C7„/i-''l"l. 

Similarly to Section 3.1, we write a S Sp°^^{U^). For the convergence Theorem 1, we will 
only need the class S'™™^ := 5"™™^; the classes Sp°^^ for p > will be required in the proof 
of the remainder estimates of Theorem 2. 

Lagrangian distributions. Given a phase function ip and a symbol a £ Sp°^^{Uip), consider 
the /i-dependent family of functions 

u{m; h) = h-^/^ [ e*^(™''')/'^a(m, 9; h) d9. (3.4) 

We call u a Lagrangian distribution of type p generated by ^p. Using the method of non- 
stationary phase, we can see that if supp a is contained in some /i-independent compact set 
K CU:^, then 

WFh{u) C {(m, a„(^(m, 9)) \ (m, 9)eC^nK}c A^. (3.5) 
The principal symbol of u is the function 

a^{u) G 5--P(A^) 

defined by the formula 

a^{u){m,dm^{m,9);h) = a{m,9;h), {m,9)eC^. (3.6) 



16 



SEMYON DYATLOV AND COLIN GUILLARMOU 



That a^{u) does not depend on the choice of a producing u, up to an 0{h^^'^P) remainder, 
will foUow from Proposition 3.3 and (3.9). As mentioned above, we will not attempt to define 
the principal symbol independently of the choice of ip. 

Following [GrSj, Chapter 11], we introduce a certain (local) canonical form for Lagrangian 
distributions. Fix some local system of coordinates on M (shrinking M to the domain of 
this coordinate system and identifying it with a subset of W^) and consider 

Kf = {{m, v)\m = -dyF{v), u € Up} C T*M, (3.7) 

where F is a smooth real- valued function on some open set Up C W^, such that the image 
of —dyF is contained in M. Then A^t^ is Lagrangian; in fact, it is generated by the phase 
function m ■ v + F{v), with v the oscillatory variable. One can also prove that each La- 
grangian submanifold that does not intersect the zero section locally has the form (3.7) for 
an appropriate choice of the coordinate system on M. (We will not have to work with La- 
grangians intersecting the zero section in this paper; the corresponding distributions have all 
the properties listed below, except that the normal forms (3.4) and (3.14) have to be written 
differently.) 

If h{v;h) E ^^"^([/ir) and x £ C^{M) is equal to 1 near — 9i,F(supp 6), then we can 
define a Lagrangian distribution by the following special case of (3.4): 

v{m; h) = x{m)h-'^l'^ [ e^(™-^+^(^))/''6(i/; h) dv. (3.8) 
Juf 

We need x to make v G C^{M); however, by (3.5) (or directly by the method of nonstation- 
ary phase), if we choose x differently, then v will change by 0{h°^)c^ . 

If V is given by (3.8), then we can recover the symbol b by the Fourier inversion formula: 

^) ^ (27r)-"'/i-"'/2 [ e-'-^-'^l^vim; h) dm + 0(/i°°)^.(Kd), (3.9) 
Jm 

here ^ denotes the space of Schwartz functions. Note that v = 0{h°°)c^ implies 6(i^; h) = 
0{h°°)c^. Moreover, if v G C^(M) satisfies (3.9) for some h £ 5p°°'P(C/f), then v is given 
by (3.8) modulo C'(/i~)c-- 

Any Lagrangian distribution can be brought locally into the form (3.8): 

Proposition 3.3. Assume that ip is a phase function, and the corresponding Lagrangian A = 
can be written in the form (3.7). For a{m,6;h) G S'p°™''(C/<^) and b{v\h) S S'^p""^^ {U p) , 
denote by Ua and vi, the functions given by (3.4) and (3.8), respectively. Then: 

1. For each a G 5'p°™''([/y,), there exists b G Sp°'^'^{Uf) such that Ua = Vh + 0{h'^)c^ . 
Moreover, we have the following asymptotic decomposition for b: 

b[u-h)= h^LXm,9;h) + Oih^'^''^P^)s^o^p^u^), (3.10) 

0<j<N 

where each Lj is a differential operator of order 2j on U^, and (m, 9) G is the solution to 
the equation {m,dmfi'm,9)) = i—duF{i/),u). In particular, if a^{u) is given by (3.6), then 

a^{u){-d,F{u), V- h) = f^pbiu; h) + ©(/i^-'Os— ^(1/^), (3-11) 
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where f^p is some nonvanishing function depending on ip and the choice of the coordinate 
system on M . Adding a certain constant to the function F , we can make f^p independent 
ofh. 

2. For each b G Sp""^^ (Up) , there exists a G Sp°™'^{U^) such that Vh = Ua + 0{h°°)c^- 

Proof. We follow [GrSj, Chapter 11]. For part 1, we apply the method of stationary phase to 
get (3.9) and then use the Fourier inversion formula. For part 2, we take some oq G S^f^^^{U^p) 
satisfying (3.11) and define Ua^ by the formula (3.4) using the symbol oq. Then v^, = Uao+'^feu 
where Vb-^ has the form (3.8) with the symbol bi = 0{h^~'^^). Repeating this process with 
in place of f^, we can write for each A^, 

0<j<N 

where each Ua^ has the form (3.4) with the symbol aj = 0{h^^^~'^P^), and v^^f has the 
form (3.8) with the symbol 5jv = ©(/i^^^-^p)). jf 

a ~ ^ ■ aj is an asymptotic sum, then 

Vb-Ua = Oih^)c^. □ 

We can now give 

Definition 3.4. Let A C T*M be an embedded Lagrangian submanifold. We say that an 
h- dependent family of functions u{m;h) G Cq^{M) is a (compactly supported and compactly 
microlocalized) Lagrangian distribution of type p associated to A, if it can be written as a sum 
of finitely many functions of the form (3.4), for different phase functions (p parametrizing 
open subsets of A, plus an 0{h°°)c^ remainder. Denote by /p°™^(A) the space of all such 
distributions, and put /'=°"^p(A) := /^"'^(A). 

We can write any u G /p°™^(A) as the sum of Lagrangian distributions associated to a 
given finite open covering of WF^(tt) in A; by Proposition 3.3, u is a sum of functions of the 
form (3.8). Moreover, if is a phase function and u G /p°™^(A(^), then u can be written in 
the form (3.4) for some symbol a, plus an 0{h'^)c^ remainder. The symbol cj(p(ti), given 
by (3.6), is well-defined modulo 0{h^-'^P). 

The action of a pseudodifferential operator on a Lagrangian distribution is given by the 
following proposition, following from Proposition 3.3 and the method of stationary phase: 

Proposition 3.5. Assume that u G /™"P(A) and P G ^^(M). Then Pu G Ip"''^^{A). 
Moreover, 

1. Assume that A = A^ for some phase function ip. Then 

a^iPu) = ct(P)|a, • a^iu) + 0{h'''P)s;on.,^^y 

2. Assume that A = Ap is given by (3.7) in some coordinate system on M. Let b{v]h) 
and b^{v;h) be the symbols corresponding to u and Pu, respectively, via (3.8). Let also 
P = Op f^{p) for some quantization procedure Op/j. Then we have the following asymptotic 
decomposition for b^ : 

b''{u;h) = ^ /i^L,(p(m,z.';/i)6(z.;/i))U,=,,™=_a„^(,) + 0(/i^(i-2''))5co.p(^^), 
0<i<Ar 
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where each Lj is a differential operator of order 2j on M x Up x Up- 

Finally, we give the following estimate of the norm of a Lagrangian distribution, fol- 
lowing from the boundedness of the Fourier transform on : 

Proposition 3.6. Assume that u £ /p°™^(Ai?), where Ap is given by (3.7). Assume that u 
is given by (3.8), with h{u\h) the corresponding symbol. Then 

\\u{m-h)\\L^<C\\h{vM\L^{U^y (3.12) 
Here C is a constant independent of h. 

Fourier integral operators. A special case of Lagrangian distributions are Fourier integral 
operators associated to canonical transformations. Let M, M' be two manifolds of the same 
dimension d, and let k be a symplectomorphism from an open subset of T*M to an open 
subset of T*M' . Consider the Lagrangian 

= {(m, v; m\ -u') \ K{m, v) = (m', u')} C T*M x T*M' = T*{M x M'). 

A compactly supported operator U : V'{M') — )• C^{M) is called a (semiclassical) Fourier in- 
tegral operator of type p associated to k, if its Schwartz kernel Kij{m, m!) lies in /i~'^/^/p°™^(Ak) 
We write U £ /p°™'^(K). Note that we quantize a canonical transformation T*M ^ T*M' 
as an operator P'(M') — >• Cq°{M), in contrast with the standard convention, which would 
quantize it as an operator V'{M) C^(M'). The /j-'^/^ f^^t or is explained as follows: the 
normalization for Lagrangian distributions is chosen so that || 1^11^2 ~ 1, while the normaliza- 
tion for Fourier integral operators is chosen so that ||C^||L2(^f/)_j.^2(A,/) ~ 1. 

After sufficiently shrinking the domain of k and choosing an appropriate coordinate system 
on M' (replacing M' with the domain of this coordinate system and identifying it with a 
subset of M"^), we can find a generating function S{m, u') for k; that is, 

K{m, v) = {m' , i^') <^=^ dmS{m,h'') = u, di^iS{m,u') = m' . (3.13) 

Here {m,v') vary in some open set Us C M x M*^. The phase function S{m,iy') — m' • v', 
with I/' the oscillatory variable, parametrizes A^ and for U £ /p°™''(«;), we can write similarly 
to (3.8), 

Ku{m, m') = h-\{m') [ es(^(™''^')"™'-'^')6(m, u'; h) dv' + C'(/i°°)c-, (3.14) 

for some symbol b £ Sp°"^^{Us) and any x ^ C^{M') such that x = 1 near the set 
d^iS{suppb). The function b is determined uniquely by U modulo 0(^°°)sp°™P(c/s)) simi- 
larly to (3.9). Note that if k is the identity map, then S{m, u') = m - v' and we arrive to the 
quantization formula for a semiclassical pseudodifferential operator. 

Similarly to Proposition 3.5, we have 

Proposition 3.7. Assume that U £ Ip°'^^{k) and P £ ^'^(M'). Then UP £ I™"P(k). // 
moreover k is given by (3.13), b{m,iy';h) and b^{m,iy';h) are the symbols corresponding to 
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U and UP, respectively, via (3.14), and P = Opf^{p) for some quantization procedure Opf^, 
then we have the following asymptotic decomposition for : 

b^{m,v']h) = ^ /i-'Lj(p(m^P)6(m,z^'))b=i.',m'=a,,5(m,^') + ^(^^^^~^'''')s^°'"''(;^s)- 
0<i<Af 

Here each Lj is a differential operator of order 2j on M' Us- In particular, 

b^{m, h) = p{dyiS{m, v'), v; h)b{m, v'\ h) + 0{h^~'^'') s'p°'^'^(Us)- 

A similar statement is true for an operator of the form PU, where P G ^^(M); the terms 
of the asymptotic decomposition have the form Lj{p{fh, u)b{m, v')), where we take fh = m 
and V = dmS{m,iy'). 

3.3. Schrodinger propagators. In this subsection, we assume that {M,g) is a complete 
Riemannian manifold and A = is the corresponding (nonnegative) Laplace-Beltrami 
operator. Let p be the semiclassical principal symbol of h'^A. G ^'^(M); we have p{m,u) = 
where \iy\g is the norm of the covector u £ T^M induced by g. We use the notation 

S*M = p~^{l) C T*M 

for the unit cotangent bundle. The geodesic flow g^ on T*M is related to the Hamiltonian 
flow e^^p of p by the formula g^ = e*^p/^. Note that g^ is a canonical transformation. 
The operator A is essentially self-adjoint on L?'{M) by [Ch] and its domain is given by the 
Friedrichs extension. Let 

U{t) = e^*^^/2 ^ (h2A/2) 

be the semiclassical Schrodinger propagator; it is a unitary operator on L'^(M). The basic 
microlocal properties of U (t) are given by the following 

Proposition 3.8. For each t G M, 

1. (Egorov's Theorem) For each compactly supported A G ^^^^ (M) , there exists com- 
pactly supported A^ G ^'p°™^(-/Vf) such that 

U{t)AU{-t) = A* + 0{h°°)L2^L2. (3.15) 

Moreover, WF;,(A*) C g-\W¥h{A)) and a{A'^) = a{A) o g^ + 0{h^-^P). 

2. (Microlocalization) The operator U{t) is microlocalized on the graph of g~^ in the 
following sense: if A,B G ^^(M) are compactly supported, at least one of them is compactly 
microlocalized, and 

g\WFh{A))nWFh{B) = (!}, (3.16) 

then AU{t)B = 0{h°°)L2^L2. 

3. (Parametrix) If A £ \I/'^°™p(M) is compactly supported, then U{t)A is the sum of a 
compactly microlocalized Fourier integral operator (of type 0) associated to , as defined in 
Section 3.2, and an 0{h°°)£^2^i2 remainder. 

The proofs are standard; part 1 can be found in [Zw, Theorem 11.1] (with the mildly exotic 
classes handled as in Appendix C), part 2 follows directly from part 1, and part 3 is 

proved similarly to [Zw, Theorem 10.3]. The operator U{t)A quantizes 5*, not g^*, because 
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of the convention adopted in Section 3.2 that a canonical transformation r*M — )• T*M' is 
quantized as an operator V'{M') — > C^{M). 

Egorov's theorem until the Ehrenfest time. Proposition 3.8 is valid for bounded times 
t; as t — )• oo, the constants in the estimates for the corresponding symbols will blow up. 
However, it is still possible to prove Egorov's Theorem for t bounded by a certain multiple of 
log(l//i), called the Ehrenfest time. To define this time, we fix an open bounded geodesically 
convex set U C M and define the maximal expansion rate 

Amax := limsup log sup \\dg\m,u)\\. (3.17) 

|t|->-oo 1*1 meU,\i^\g^l, 

Here ||(i5(*(m, z^)|| is the operator norm of the differential 

with respect to any given smooth norm on the fibers of T(T*M) (e.g. the norm induced by 
the metric ^f). 

Since we will work on a noncompact manifold, we introduce cutoffs into the corresponding 
propagators: 

Proposition 3.9. Assume that Xi,X2 G "^^{M) satisfy \\Xj\\i2_^i2 < 1 + 0{h) and are 
compactly supported inside U. Let £e> ^ and take Aq, Aq > such that 

Ao > A'o > (1 + 2ee)A 

max- 

Fix 4o £ I^- Then for each integer 

/ G [0,log(l//i)/(2|to|Ao)], (3.18) 
and each compactly supported A S ^™™p(M) with 

W¥h{A) C £e, ■.= {l-ee<\v\g<l+ £e}, 

the compactly supported operator 

A« := (X2C/(to))'A(C/(-to)Xi)' 

lies in ^'p°™'^(M), modulo an 0{h'^)i^2^i2 remainder, with 

PI = l\to\Ao/\og{l/h) < 1/2. (3.19) 

Moreover, the 5^°™'' seminorms of the full symbol of A^''^ are bounded uniformly in I, in the 
following sense: the order k derivatives of this symbol are bounded by Ch~^P' , where C is a 
constant independent of h and I. The principal symbol of A^'"') is 

l-i 

a(.4(')) = (a(A) o/«) ll{a{Xi)aiX2))og^'o + Oih'-^^^). 
j=0 

The wavefront set of A^^\ for I > 0, is contained in WFh(Xi) Pi WF/j(X2) riSe^- Finally, if 
Ua and Ux are open sets such that WFh{A) C Ua and WFh{Xi) n WF/i(X2) C Ux, then 
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A^') is micro supported, in the sense of Definition 3.1, inside the set 

i-i 

y(0 ■.= g-^'<^{UA)nf]g-^''>{Ux). 

j=0 

The set V^''^ does not depend on h directly, however it depends on I, which is allowed to 
depend on h, and our microlocal vanishing statement is uniform in I. 

Proposition 3.9 is the main technical tool of obtaining the polynomial remainder bound 
of Theorem 2; it is also the reason why the classes appear. Its proof, following the 

methods of [AnNo, Section 5.2] and [Zw, Theorem 11.12], is given in Appendix C. We do not 
impose any restrictions on the set U at this point, however in our actual argument it will 
have to contain a neighborhood of the trapped set - see the beginning of Section 5.2. 

Propagating generalized eigenfunctions. The following fact, similar to [Dy2, Proposi- 
tion 3.3], will be used to propagate the Eisenstein functions by the group U{t): 

Lemma 3.10. Assume that u G C°°{M) solves the equation 

{h^A - z)u = 0, |1 - z| < Ch. 

Let X & C^{M); take t G M and assume that xt G C'^{M) is supported in the interior of a 
compact set C M and satisfies 

'^g(suppx,supp(l - Xt)) > \t\. (3.20) 

Here dg denotes Riemannian distance on M . Then 

= xe-''''^^''^U{t)xtu + 0[h'^\\u\\L2^^K,))LHMy 

Proof. Without loss of generality, we assume that t > 0. For < s < t, define 

Us = x{u-e-''^/^^''^U{s)xtu). 

We need to prove that 

\\ut\\L^=0{h^)\\u\\L^K,y (3.21) 

Since x = XXt-, we have = 0; next, 

2hDsUs = -xe-''''^^''^U{s){h'A - z)xtu 

= -e-'''/^^^^xUis)[h^A,xt]u. 

Let B G ^f'^°"^P be compactly supported inside Kt x Kt, equal to the identity microlocally 
near suppxt H S*M, but microlocalized in a small enough neighborhood of S*M so that 
by (3.20), 

ff'(supp x) n WFh{B) n supp(l - Xt) = 0. 
Note that WF h{[h'^ A, xt]) C supp(l -Xt)- Then by part 2 of Proposition 3.8, 

\\xU{s)[h^A, Xt]Buh2 = 0{h^)\\uh2^K^y 0<s<t. (3.22) 
Moreover, by Proposition 3.2 

\\xU{s)[h'A,xt]{l - B)u\\l2 = 0{h^)\\u\\L2(^j,^y (3.23) 
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Combining (3.22) and (3.23), we get ||9sns||2,2 = C'(/i°°)||u||^2(;^^); it remains to integrate in 
s to get (3.21). □ 

Hilbert— Schmidt norm estimates. We now prove Hilbert-Schmidt norm estimates for the 
product of a pseudodifferential operator with a spectral projector. (See [HoIII, Section 19.1] 
for the properties of Hilbert-Schmidt and trace class operators.) To simplify notation, we 
consider a spectral interval of size h centered at A = 1; similar statement is true for the 
interval [A + cih, A + C2/1] with A > 0, replacing 5"*^^ by XS*M. 

Lemma 3.11. Fix ci,C2 G M and let l^i_^_^^f^^i_^_^^^{h'^ A) be defined by means of spectral 
theory. Assume that A £ ^'p°™''(M) is compactly supported. Then 

/i('^-i)/2||]l[i+,,^,i+e,h](/^'A)^||Hs<C||a(A)b2(5.M) + O(/i'-'0- (3.24) 
Here C is a constant independent of A (ifWFfi{A) is contained in a fixed compact set), 
however the constant in 0{h^~'^'') depends on A. We take the norm of cr{A) on the 
energy surface S*M with respect to the Liouville measure ^l- 

Moreover, ifWFhiA) is microsupported, in the sense of Definition 3.1, in some h-dependent 
family of sets V{h) C T*M, then 

ij^^^^^^^^^^^j(/i2A)^||Hg < CiiL{V{h) n S*Mfl^ + 0{h^). (3.25) 

Here ^L{V{h)r\S*M) denotes the volume ofV{h)r\S*M with respect to the Liouville measure 
on S*M and the constant C depends on a certain Sp°^^ -seminorm of the full symbol of A. 

Proof. Take a function x £ o5^(M) such that x is compactly supported in some interval 
(— T, T) and x does not vanish on [ci, C2] (for example, take nonzero Tp E Co°(M) with ^ > 0, 
then I'i/'l > in an interval [cie,C2e]; set x{^) ■= i^i^^))- Then 

]l[l+,,;,,l+e2/x](/i'A) = Zxiih^A - l)/h), 

where Z is a certain function of h^A and it is bounded on L'^{M) uniformly in h. It then 
suffices to estimate the Hilbert-Schmidt norm of 

B = h^'^-^^/^x{{h^^ - ^)/h)A = (27r)-i/i('^-i)/2 ^ ^(t)e-'^^l'^U{2t)Adt. 

J-T 

Let Aq £ ^'q°™^(M) be compactly supported and equal to the identity microlocally near 
WFfi{A). By part 3 of Proposition 3.8, for each t we have 

U{2t)Ao = U2t + R2U 

where U2t E /c°™P((^2t-j compactly supported Fourier integral operator and i?2t = 
0{h°°)L2^L2- Then 

{U{2t) - U2t)A = 0{h^)us. (3.26) 
Indeed, we can write the left-hand side of (3.26) as the sum of i?2t^ and U{2t){l — Aq)A; 
it remains to note that R2t = C'(/i°°)l2_>/^2, ||j4||hs is polynomially bounded in h, and 
\\il-Ao)A\\ns = 0{h^). 

By (3.26), we can replace U{2t) by U2t in the definition of B. Now, the Hilbert-Schmidt 
norm of B is equal to the L^{M x M) norm of its Schwartz kernel Kb- Using the local normal 
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form (3.14) for Fourier integral operators, we can write Kb, up to an 0{h°°)cg° remainder 
and an appropriate cutoff in the m' variable, as a finite sum of expressions of the form (in a 
fixed coordinate system on M) 

J-T Jwi 

Here S{m, v'; 2t) is a generating function for (7^* and 6 is a certain symbol in Sp°^^. Moreover, 
b admits an asymptotic expansion in terms of the full symbol of A, by Proposition 3.7. The 
fact that S and b can be choosen to depend smoothly on t follows from the proof of part 3 
of Proposition 3.8. 

We can deduce from (3.13) that 

g'^\m, 1^) = (m', z^') =^ dt{S{m, u'; 2t)) = p{m, u). (3.28) 

The equation (3.28) is true for the particular generating function constructed in the proof of 
part 3 of Proposition 3.8. One can add any function of t to the function S and still obtain 
a generating function of (7^*; however, the amplitude b with respect to the new generating 
function will no longer be a symbol, as its derivatives in t will not be bounded uniformly in 
h. It follows from (3.28) that the function 

^>(m, m', v', t) = S{m, v'; 2t) - m! ■ v' - t 

is a nondegenerate phase function (with m,m' as base variables and z^',t as the oscillatory 
variables) and generates the (immersed) Lagrangian 

A = {(m, u; m , —v') \ p{m, u) = 1, 3t G {—T, T) : g^*{rn, v) = (m , i'')}. 

Then (3.27) lies in /p°™^(A). By the local normal form (3.8) of a Lagrangian distribution, 
we can write (3.27), up to an 0{h°°)cg° remainder and an appropriate cutoff in the {ra^m!) 
variables, as the sum of finitely many expressions of the form 

h-'^ [ e'^"'-''+"''-'''+^'-''''''^^/^b{iy,iy';h)dudu', (3.29) 

where F parametrizes some open subset of A by (3.7) and 6 is a symbol in S'p°™''. By 
Proposition 3.7 and Proposition 3.3, we see that the symbol b has the following asymptotic 
expansion in terms of the full symbol a of j4: 

b{v, h)= WL ja{m' , v'; 

o<i<v 

where each Lj is a differential operator of order 2j and m, m! are given by the relation 
{■m,i/,m' £ A; in particular, {m\i'') G S*M. 

We now use Proposition 3.6 to estimate the norm of (3.29); as B is, modulo C'(/i°°)hs5 
a sum of operators with Schwartz kernels of the form (3.29), this would give an estimate on 
the Hilbert-Schmidt norm of B. For (3.24), we can write 6(z^, u'; h) as a multiple of a{m' , v') 
plus an 0{h}~'^'^) remainder and note that {m',u') always lies in S*M. For (3.2-5), we use 
that b = 0{h°°) outside of the preimage of V{h) under the map {u, u') 1— t- (m', v'), and that 
sup \b\ can be estimated by a certain 5p°™^-seminorm of a. □ 
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Local traces of integrated Schrodinger propagators. We give the following version of 
the Schrodinger propagator trace formula in the case where there are no contributions from 
closed geodesies: 

Lemma 3.12. Assume that M is a d- dimensional complete Riemannian manifold and Bs 
is a family of compactly supported pseudodifferential operators in ^'p°™^(M), smooth and 
compactly supported in s £ {—Tq,Tq), where Tq > is fixed. Assume also that all Bg are 
micro supported, in the sense of Definition 3.1, in some h- dependent family of hounded sets 
V{h) C T*M, and the following nonreturning condition holds: 

{m,v) G V{h), \s\ < To ^ d{{m,v),g'{m,v)) > C-^\s\hf. (3.31) 

Here C is some constant and d denotes some smooth distance function on T*M . Let Bs = 
Op^(6(s)) + 0(/i°°)^-oo for some family of symbols b{s, m, u) and some quantization procedure 
Op/j. Then for each N and each A > 0, we have the trace expansion 



/Jo 
^-i\^s/{2h) Tr{Uis)Bs) ds 

(3.32) 

V h^ Ljb{0,m,Xi^)dfiLi'm,u) + O{h^^^-'^P^)c 
.J^AT Js*M 



0<j<N 

where fi^ is the Liouville measure and each Lj is a differential operator of order 2j on 
T*M(,„^) X (— ro,To)s, independent of Bg and depending smoothly on A. In particular, 

Proof. As in the proof of Lemma 3.11, we can reduce to computing the trace of the operator 
with the Schwartz kernel (in some fixed local coordinates) 

K{m,m') = i2TTh)-^ / eT^^'^^"''''''^'^-"''-'''-^^'/%{m,i^' ,s;h) du'ds, 

J -To JR'i 

where S{m, u'; s) is a local generating function for g'^ in the sense of (3.13) and b{m, v' , s; h) 
is a certain symbol in Sp having an asymptotic expansion in terms of the jet of hs at the 
point {dyiS{m, u'; s), v'). The trace of the corresponding operator is 

To 

K{m, m) dm = {27Th)-^ J J erAS(my;s)-m-u'-x^s/2^^^ ^/^ ^. dmdu'ds. 

M -To MxRd 

We now use the method of stationary phase. The stationary points of the phase are solutions 
to the equations g^{m, v') = (m, v') and \i''\g = A; they occur at s = and may also occur for 
> Tj, where > is the injectivity radius of M . For A|s| > r,/2, we see by (3.31) that 
the expression under the integral can be split into two pieces, on one of which the symbol is 
0{h°°) and on the other, the differential of the phase function has length at least C~^hf; by 
repeated integration by parts, the latter integral is 0{h°°). 

It remains to evaluate the contribution of the stationary set {s = 0} n XS*M. The phase 
function is degenerate on these points; however, one can pass to polar coordinates u' = ruj, 
with \uj\g = 1 and r > 0, and apply the method of stationary phase in the {r,s) variables. 
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resulting in the expansion (3.32). See for example the proofs of [Ro, Theoreme V-7 and 
Proposition V-8] or [RoTa, Lemma 3.1] for details of the computation. □ 

4. General assumptions 

In this section, we list the general geometric assumptions on the manifold M and analytic 
assumptions on its Laplacian required for our results to hold. As noted in the introduc- 
tion, they are satisfied in particular if outside of a compact set, M is isometric to either 
the Euclidean space (studied in Section 6) or an asymptotically hyperbolic space of constant 
curvature (studied in Section 7). We also derive some direct consequences of the general 
assumptions, including averaged estimates on plane waves and the existence of limiting mea- 
sures /i^. 

4.1. Geometric assumptions. In this subsection, we specify the geometry of the manifold 
M at infinity. 

Let us introduce some notation and terminology first. On a complete Riemannian manifold 
(M, g) we denote by 5* the geodesic flow of the metric 5, considered as a map on the cotangent 
bundle T*M. Any smooth function / on M can be lifted to a function on T*M; denote by 
/, / S C°°{T*M) the derivatives of / with respect to the geodesic flow: 

f{m,u) := dtf{g\m,iy))\t=o, f{m,u) := dff{g\m,ij))\t=o. 

We denote by S*M the unit cotangent bundle {(m, v) \ |z^|g = 1} C T*M. 

A boundary defining function on a smooth compact manifold M with boundary is a smooth 
function j; : M — )• [0, 00) such that x > on M and x vanishes to first order on dM. 

We make the following assumptions: 

(Gl) {M,g) is a complete Riemannian manifold of dimension d = n + 1. Moreover, there 
exists a compactification of M, namely a compact manifold with boundary M such 
that M is diffeomorphic to the interior of M. The boundary dM is called the boundary 
at infinity; 

(G2) There exists a boundary defining function x on M and a constant £0 > such that 
for any point (m, i^) G S*M, 

if x{m, I') < £q and x{m, v) = 0, then x(m, v) < 0; (4-1) 

(G3) For each (m, i^) G S*M such that x{m) < Eq and x{m,i>) < 0, the geodesic g^{m,i>) 
(projected onto the base space M) converges as t — )• +00, in the topology of M, to 
some point ^^oo{it^, i^) £ dM. The function ^_|_oo depends smoothly on (m, u), and we 
extend it naturally (as the limit of the corresponding geodesic) to a smooth function 
on S*M \ r_, with r_ given in Definition 4.1 below; 

(G4) There exists an open set Uoa C M x dM such that C/qo contains a neighbourhood of 
{(^,^) G Af X dM I ^ G dM} and a smooth real-valued function (j){m,(^) = (j)^{m) on 
Uoo such that |f?m0g(?7i)|g = 1 everywhere and the function 

T(m,0 := {m,dmM^)) ^ S*M, (m,a G Uoo, (4.2) 
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dM dM 

Figure 2. Left: an illustration of (G2), showing a forbidden geodesic. Right: 
an illustration of (G3). The point {m,v) does not escape directly in the 
forward direction, but the point (m', u') does. 

is a diffeomorphism from onto with inverse given by 

where the sets and are defined by 

■■= {KO e Uoo I x{m) < £0, x{T{m,0) < 0}, 
y+ := {{m,i^) e S*M I x(m) < eq, x{m,i^) < 0, (m, ^+oo(m, z^)) e f/oo}; 
(G5) if (m, z^) G F+, then v) £ V+ for all t > 0; 

(G6) if G dM and m,m' G M are such that (m, (m', (^) G Ul^, then d^4)^{m) = 
d^(j)^{m') if and only if r(m,^) and T{m',^) lie on the same geodesic. Moreover, the 
matrix dmd^(j)^{m) has rank n. 

Escaping trajectories and the trapped set.. We now define the incoming/outgoing tails 
F-i- and the trapped set K: 

Definition 4.1. Let 'y{t) be a unit speed geodesic. We say that it escapes in the forward, 
respectively backward, direction, if ^{t) goes to infinity in M as t ^ +oo, respectively t — )■ 
— oo. If ^{t) does not escape in some direction, we call it trapped in this direction. Denote 
by r+ d S*M the union of all geodesies trapped in the backward direction, by r_ the union of 
all geodesies trapped in the forward direction, and put K = n r_ ; we call K the trapped 
set. 

An escaping geodesic could potentially spend a long time in the compact part of the 
manifold. It is helpful to consider geodesies that escape in a straightforward way (with the 
boundary defining function x decreasing along them); they appeared in assumption (G3) for 
instance. 

Definition 4.2. We say that {m,v) £ S*M directly escapes in the forward, respectively 
backward, direction, if x{m) < eq and x{m,u) < 0, respectively x{m,v) > 0. Here £q is the 
constant from (G2). Denote byDS^, respectively - , the set of all points directly escaping 
in the forward, respectively backward, direction. 
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One can verify that T± are closed sets and the trapped set K is compact (see [GeSj, 
Appendix]); in fact, since S*M n {x < eo} C VS^ U VS^, we have K C {x > eq}- 

For the example of M = M""*"^ discussed below, we have F-i- = 0. The point (m, u) lies in 
D£-^- if and only if x{m) < Eq and m ■ v > 0. 

Comments on the geometric assumptions. 

A basic example to have in mind for a manifold satisfying our assumptions is M = ]R"+^ 
with the radial compactification M being a closed ball and the boundary at infinity dM 
equal to the sphere S". We will often use this example to illustrate the somewhat abstract 
assumptions of this section. (A more general version will be considered in Section 6.) 

An important corollary of the assumption (G2) is that for e < Eq, the compact set {x > 
e} C M is geodesically convex; i.e., if 7(i) is a geodesic and 7(ti),7(t2) ^ {x > e}, then 
7(t) G {x > e} for t G [ti, f2]- For the example of M = M"+-^, we can take x = (l + |m,|~^)^^/^, 
where |m| is the Euclidean length of m G W^; the corresponding sets {x > e} are balls centered 
at zero. 

It also follows from (G2) that for (m, v) G the function x{g^{m, v)) is decreasing for 

t > 0. One can show that x{g*{m, v)) — )• as i — )• +oo and thus (7*(m, z^) escapes in the 
forward direction; we do not give a proof of this fact as it follows from the more restrictive 
assumption (G3). Also, if a geodesic 7(t) escapes in the forward direction, then for t large 
enough we have 7(t) G T)£^. For M = M"^^, we have ^+oo("i-5 1^) = G S". 

Assumption (G4) means that for m sufficiently close to the infinity, the covectors z/ such 
that (m,i/) G T)£^ are in one-to-one correspondence with the limit points ^_|_oo(m, i^), and 
the inverse correspondence can be described by a phase function. It follows in particular 
from (G4) that for a fixed ^ G dM, the set of directly escaping points (m, u) such that 
C+ooim,!^) = C and (m, ^) G Uoo is the intersection of 2?<S+ with the Lagrangian 

Ag := {(m, | (m, G ^oo}- (4.3) 

In the model case M = ]R"+^ we can put for any R > 0, Uoo = {(w-,^) | |?7i| > R}, 
and (l)({m) = m ■ ^, so that r is the canonical map from M"'+^ x S" to 5*M"^"'^. Then 
f7+ = {{m, ^) I \m\ > R, m • ^ > 0} and = {{m, u) \ \m\ > R, m ■ u > 0}; the difference 
is that is considered as a subset of M"+^ x S", while is considered as a subset of 

The condition (G6) is required in Proposition 5.12. To explain it, note that under the 
assumption (G4), if (m,^) G and {m{t),u{t)) = g^{T{m,£^)), then 

dt4)i^{m{t))\t=o = dm4>^{rn) ■ dtm{t)\t=o = g{dm(p^im),dm<j)^{m)) = 1. (4.4) 

Therefore, d^(j)^{m) is constant on the geodesic passing through T(m,^). 

4.2. Analytic assumptions. In this subsection, we formulate the analytic assumptions 
on plane waves. Let M be as in the previous subsection, A be the (nonnegative definite) 
Laplace-Beltrami operator on M, and /i > be the semiclassical parameter. We make the 
following assumptions: 
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(Al) There exists cq > (equal to for the Euchdean and to for the hyperbohc case), 
such that for each X > 0, h > and ^ G dM, there exists a function, called distorted 
plane wave, Efi(X,S^;m), smooth in all variables and solving on M the differential 
equation (1.4) in m: 

Here ^ gives the direction of the plane wave, while A corresponds to its semiclassical 
energy; 

(A2) for each < Ai < A2, the Schwartz kernel of the semiclassical spectral projector 

n[Ai,A2] := l[A2+coh2,A2+co/i2](^^^) 
can be written in the form 

n[A„A2]("^,^') = (27r/i)-"-i t\^fn{X/h) [_EhiX,C;m)Eh{X,C;m')d^dX. (4.5) 

JXi JdM 

Here integration in ^ is carried with respect to a certain given volume form on 
dM and fn{z) > is a smooth function of z such that \d^fu{z)\ < Ck{z)~^ for each 
k and /n(-z) — >• 1 as z — )■ c«. 

We now assume that plane waves admit the decomposition 

Eh{X, i- m) = xo(m; 0^°(A, m) + ^^(A, m), (4.6) 

where xo-,E'y^iE\ are respectively a cutoff function, an explicit 'outgoing' part of the wave, 
and the 'incoming' part, satisfying more precisely the following properties: 

(A3) Xoi^^'iO is a function smooth in m G M and ^ G dM, supported inside the set C/qo 

from (G4) and Xo("i)0 = 1 for m sufficiently close to ^; 
(A4) E^{X, ^; m) is a smooth function of A G M and (m, ^) G C/00, of the form 

El{X, m) = eX^«M60(A, ^, ^; h), (4.7) 

where Uoo and 0^ are defined in (G4) and is a classical symbol in h defined for 
A G M and (m,^) G J/oo! that is, is smooth in all variables, including h, up to 
h = 0. We also require that b^{X, ^, m; 0) is independent of A; 
(A5) for A in a fixed compact subset of (0,oo) and Eq defined in (G2), the function 

1 + \\Eh{X,^;m)\\L2(^{^>^^}) 

is /i-tempered in the sense of (3.2); 
(A6) for A in a fixed compact subset of (0,oo), each ^ G dM, and each (m, Az^) G 
WF/i(i?^(A, ^)), we have {m,v) G S**-/!^ and either the geodesic 7(t) = g^{m,v) does 
not escape in the forward direction (i.e. (m, i^) G r_) or there exists t > such that 
7(t) lies in the set 

:= {{m, dm4>£,{m)) \ m G supp(amXo)}- (4.9) 

The constants in the corresponding estimates (in the definition of the wave front set 
of a distribution given in Section 3.1) are uniform in A and ^; 
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(A7) there exists ei G (0,eo) such that for {m,v) G S*M directly escaping in the forward 

direction and x{m) < ei, the point (m, ^^cx>{m, z^)) lies in the set Uoo defined in (G4) 

and Xo = 1 near this point; 
(A8) Let r : — )• be the diffeomorphism from (G4). Then its Jacobian with respect 

to the volume measure dvol(m)(i^ on and the Liouville measure on V^, is equal 

to |6°(l,C,"i;0)p, with 6° defined in (A4). 

For example, for M = M""^^ we put cq = 0, Efi{X,^;m) = e*'^?'™/'^ and use the standard 
volume form on the sphere dM = The equation (4.5) then follows from the Fourier 
inversion formula. 

Let us informally explain how the decomposition (4.6) is constructed and provide a justi- 
fication for assumptions (A3)-(A6), putting for simplicity A = 1. First of all, (A4) implies 
that for any x G Cq^{M), xXo^h^ ^ a function of m, is a Lagrangian distribution associated 
to the Lagrangian Ag from (4.3). In fact, in the cases considered in the present paper, 
solves on its domain the equation (1.4); however, we do not make this assumption here, as in 
more complicated cases (such as asymptotically hyperbolic manifolds of variable curvature) 
might only be an approximate solution to (1.4) in a certain sense. 

If we assume that E'^ solves (1.4) on its domain, then the function 

Fh{X, e; m) = {h^A - A2 - co/i2)(xo(m)ii;0(A, C; m)) 

is equal to [h'^A, Xo]E'jl. Since E'j^ is a Lagrangian distribution associated to A^, the wavefront 
set of Fh is contained in W^. We will now take Ej^ = —Rfi{X)Ffi, where Rh{X) is the incoming 
scattering resolvent, a certain right inverse of h'^A — A^ — cq/i^. Moreover, in our cases Rh{X) 
will be microlocally incoming in the weak sense: if we multiply F^ by a (possibly small) 
constant to make Rh{\)Fh bounded polynomially in h, then each point in the wavefront 
set of Rh{X)Fh, when propagated forward by the geodesic flow, will either converge to the 
trapped set or pass through WF^iF^). Thus, the assumption (A6) should be viewed as a 
direct consequence of the fact that the scattering resolvent is microlocally incoming and of 
propagation of singularities. 

The assumption (A7) looks less natural, but will play an essential role in our proofs, 
in Propositions 5.2 and 5.5. It holds for both Euclidean and hyperbolic infinities, but for 
different reasons. For the hyperbolic infinity, Xo("; '^) is equal to 1 in a small neighborhood of 
^ and one can see that for (m, v) directly escaping in the forward direction and converging 
to 1^, the distance from m to ^ in M is 0{x{m)). This is not true in the Euclidean case; 
however, in that case xo is equal to 1 outside of a compact subset of M (that is, near the 
whole boundary dM, not just near ^). 

The assumption (A8) is required to relate the natural measure arising from the function 
E^ to the Liouville measure. If E^ were equal to E^, then this assumption would simply 
follow by taking the trace in (4.5) with a compactly supported pseudodifferential operator 
and a smooth cutoff function in A. 

4.3. Limiting measures. We now define the family of limiting measures /i^. These mea- 
sures result from propagating the natural measure arising from the 'outgoing' part E'j^ of 
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the plane wave, which is supported on the Lagrangian from (4.3), backwards along the 
geodesic flow. In contrast with [Dy2], where the exponential decay of the measure along the 
flow ensured its convergence, our measures will only be defined for almost every ^. 

We first define the measure /ig on S*M, corresponding to E'^, as follows: for each compactly 
supported continuous function a on S*M, put 

adfl^= I |6°(l,C,m;0)pa(r(m,^)) dvol(m). (4.10) 

S*AI J{m,OGU+ 

The support of /i^ is contained in the Lagrangian from (4.3) and the integral (4.10) 
depends continuously on ^. We see from (A8) that for any continuous function / on dM, 



_f{0 a{m,u)dfl^{m,u)d^_ = f{^+oc{m,u))a{m,u)d^iL{m,u). (4.11) 
dM Js*M Jv+ 

We now want to define the measure by 

adfi^= lim / aog^^djl^, (4-12) 

valid for all compactly supported continuous functions a. To show that the limit exists 
for almost every ^ (chosen independently of a) and for every a, we will use monotonicity. 
By (4.11), (G5), and using the invariance of the function ,^+oo and the Liouville measure 
under the geodesic flow, we see that if a and / are nonnegative, then 

f _f{i) I {a o g'^) dfi^d^ = I f{i+oo{m,v))a{m,v)d^iL{m,v) 
JdM JS*M Jg-Ky^) 

is increasing with t. Therefore, for each ^ the integral 



JS*M 

is increasing in t for any nonnegative a. Moreover, the integral of Ia,t{C) in ^ is bounded by 
a t-independent constant, namely by the integral of a by the Liouville measure. Taking a to 
be an approximation of the characteristic function of each member of a countable family of 
compact sets exhausting S*M, and using the monotone convergence theorem, we see that 
there exists a measure zero set X C dM such that for ^ ^ X, we have for each j and for any 
compactly supported continuous function a, 



lim / {aog *) d/2t < oo. 



This limit is a continuous functional on the space of continuous compactly supported func- 
tions on S*M; therefore, there exists unique Borel measure such that (4.12) holds. More- 
over, we see that the limit (4.12) is uniform in a, as soon as we fix a compact set containing 
suppa and impose a bound on sup^.^ \a\. One also sees immediately (1.6), namely that for 
compactly supported continuous a, 



suppan^+^(C) = ^ / adfi^ = 

JS*M 



as Jg*^j{cL o (7 *) djl^ = for all t. 

We can integrate the measure in ^, getting back the Liouville measure: 
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Proposition 4.3. For each f £ C'^{dM) and each a G C^{S*M) we have 



_f{i) I a{m,i^)dfj,^{m,iy)dC = f{^+ooim,u))a{m,i^)dfiL{m,u). (4.13) 

dM JS*M JS*M\T^ 

In particular, if ^l{T^) = (which will always he the case in our theorems, see (5.2)j, then 
J fi^ d(^ is the Liouville measure. 

Proof. The left-hand side can be written as 

hm / f{^+oo{m,u))a{m,u)dfiL{m,v). 

*^+~ Jg-^V+) 

It remains to use the dominated convergence theorem; indeed, the function under the integral 
is bounded and compactly supported, we have 5^*1 (V"+) C g'^^iV^) for ti < t2, and the 
union of g~^{V^) over alH € M is exactly S*M \ r_, as for every geodesic 7(t) escaping in 
the forward direction and for t large enough, the point j{t) is directly escaping in the forward 
direction and {-/{t),C+oo{'y{t))) £ U°°. □ 



Finally, the following lemma will be useful to relate our measure fj,^ to the one obtained 

iO 
h 



from E'^ in the proofs of Theorems 1 and 2: 



Lemma 4.4. Let ^ ^ X , so that /x^ is well-defined. Let a he a compactly supported continuous 
function on S*M . 

1. fj,^ is invariant under the geodesic flow: for each t GM, 

/ aog^dfj.^= / adfi^. (4.14) 

Js*M Js*M 

2. T/suppa C T)£j^ H {x < ei}, where T>£j^ is given hy Definition 4-2 and ei is defined 
in (A7), then 



\b^{l,^,m;0)xo{'m',0\'^'^{''T^j'^rn4>^{m))dvol{m)= / adfj,^. (4-15) 
{m,^)eUao Js*M 

Proof. 1. Follows immediately from the definition (4.12). 

2. First of all, note that for m in the support of the function a{m,dm(p^im)), we have 
(m, ^) G and Xo{'m']0 = 1 by (A7); therefore, the left-hand side of (4.15) becomes the 
integral of a over the measure fi^ defined in (4.10). By (4.12), it is enough to show that for 
t > 0, 

/ aog~*dfL^= / adfi^. 
Js*M Js*M 

For that, it is enough to show that for each / G C^{dM), 

/ fiO / aog'*dfi^= I I adfi/:. 
JdM JS*M JdM JS*M 

Using (4.11), we rewrite this as 



f{i+oo)adnL= I fi(,+oo)adnL- 
Jv+ 
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This is true as suppa C C g^^iV^)- D 

4.4. Averaged estimates on plane waves. One of the principal tools of the present paper 
are microlocal estimates on the plane waves Eh{X,^) on average in A,^, where A takes values 
in a size h interval. They are direct consequences of (4.5) and the Hilbert-Schmidt norm 
estimate (3.24). More precisely, restricting to the case \ = 1 + 0{h) for simplicity, we have 
the following 

Proposition 4.5. Let x C^{M). Then: 

1. X^[i,i+h] Hilbert-Schmidt operator and there exists a global constant C such that 
for each bounded operator A : L'^{M) — )• Lp'^M), we have 

/i-ipx(m)i^^(A,C;m)||2^^^^(^^^^^^j^_^^^j) < C/i"Pxn[i,i+,] ||^s- (4-16) 

2. The functions x^h cire bounded in on average in the following sense: there exists a 
constant C{x) such that for any h, 

/.-i||x(m)i?,(A,e;m)||^,^^^^^(^^^^^j^^^^,,^ < C{x). (4.17) 

The h^^ pref actor in both cases is due to the fact that we are integrating over an interval 
of size h in X. 

Proof. 1. It follows immediately from (4.5) that 

h-' / /n(A//i)A" / _ixEhi\, 0) ® ixEhi\ 0) d^dX = (2vr)"+i/i"xn[i,i+;.]X- (4.18) 
Jl JdM 

Here ® denotes the Hilbert tensor product, defined in (2.8). The integral on the left-hand 
side converges in the trace class norm, as the Schwartz kernels of the integrated operators 
are smooth and compactly supported. Therefore, x^[i,i+h]X is trace class. Since 

X^[i,i+h]X = {x^[i,i+yi]){x^[i,i+h])* : 

we see that xn[i.i+/i] is a Hilbert-Schmidt operator. Now, multiplying both sides of (4.18) 
by A on the left and A* on the right and taking the trace, we get 

/.-iA"/Vn(A//.)V2^x(m)i?.(A,e;m)||^^^^^^(^^,^^j^_^^,,^ 

= (2vr)"+i/i" Tr((Axn[i,i+,])(^xn[i,i+/,])*) (4-19) 

= (2vr)"+i/."Pxn[i,i+,]|||s- 

2. We would like to use Lemma 3.11 to estimate ||xn[i,i+/i] ||hs (we can put x on the other 
side of the projector in (3.24) by taking the adjoint), however this is not directly possible as 
X is not compactly microlocalized. We thus use that Eh solve the equation (1.4), writing by 
the elliptic parametrix construction (same as for the proof of Proposition 3.2) 

X = B + Qxih^l^^-X^ -CQh^) + Rx (4.20) 

for A € [1, l + /i], where B G ^'=°'"P(M), Qx G "^'"^{M), and Rx G h°°^-'^{M) are compactly 
supported and B is independent of A and equal to x microlocally near S*M. We can also 
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assume that Q\ and R\ are smooth in A. Now, we substitute (4.20) into the left-hand side 
of (4.17) and use the triangle inequality. By (4.16), the term featuring B is bounded by a 
constant times /i^Hi^IIji IIhS' which is bounded uniformly in h by Lemma 3.11. The term 
featuring Qx is zero by (1.4). 

Finally, we show that the term featuring R\ is 0{h°°). This does not follow immediately 
from (4.16), as the operator R\ depends on A. We use the following variant of (4.19): for 
A G [l,l + /i], 

Differentiating in A, we get 

2h-^Re{\^l^fn{\/hf/\d-^R-^)Eh{X)A'''^fu{\/hfl^R-^^^^^ 

We now integrate in A from 1 to \ + h. The integral of the left-hand side is bounded by 
a constant times /i"||^i+/i||hs ~ 0{h°°). The integral of the first term on the right-hand 
side is the quantity we are estimating. Finally, the second term on the right-hand side is 
bounded by a constant times Tr((5j^i?^)n[i i_|_^]i?i)|, which is 0{h°°) uniformly in A, as 
the Hilbert-Schmidt norms of both R^ and d^R^^ are 0{h°°). □ 

5. Proofs 

5.1. Proof of Theorem 1. In this section, we prove the convergence Theorem 1 under the 
following assumption: 

I^l{K) = 0, (5.1) 

where denotes the Liouville measure on S*M and K is the trapped set. First of all, note 
that (5.1) implies 

/"L(r±) = 0. (5.2) 

Indeed, fix e G (0, eo)i where eo is defined in (G2), and take the set = H {x > e}. For 
(m, u) G r+ n {x = e}, we have x(m, i^) < 0; indeed, otherwise (m, u) directly escapes in the 
backward direction and thus cannot lie in F^. It follows that g~*{T^) C F^ for t > 0. Since 
F^ is bounded, and /Ul is invariant under the geodesic flow, we have 

Letting e — > 0, we get (5.2). 

We next note that the averaged bound (4.17) on Ei^^ on compact sets, together with (1.4) 
and the elliptic Proposition 3.2, give the following 

Proposition 5.1. Assume that A G ^"^(M) is compactly supported and WFh{A)nS*M = 0. 
Then 

h-^\\{AE,,{X,0,EH{X,mLl,idM>cli,i+h]) = (5-3) 
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Therefore, it is enough to prove (1.5) for a compactly supported A G \I/^°™p(M) microlo- 
caUzed in an arbitrarily small neighborhood of S*M. 

Take t > 0; we will calculate limits of the form limf_!>+oo lim/i-5>0) thus Ot{h°°) expressions 
(that is, expressions that are 0{h°°) with the constants depending on t) will be negligible. 
Take x £ C'^(M) indepedendent of t and such that A = x^X- We first use that Eh is a 
generalized eigenfunction of the Laplacian (1.4) and apply Lemma 3.10: for each A E [1, l + /i] 
and each ^ G dM, 

xEh = xe-**(^'+^°^')/2"t/(t)xtii;h + Ot{h'^\m\L2iK,))L^. (5.4) 

Here U{t) = e**'^^/^ is the semiclassical Schrodinger propagator and xt S Cq°{M) is sup- 
ported in the interior of the compact set Kt C M and satisfies dg(supp x, supp(l — Xt)) > t- 
We also assume that \xt\ < 1 everywhere and Kt contains {x > Eq}, where Eq is defined 
in (G2). By Proposition 3.8, we can write U{-t)AU{t) = A"* + C't(/i°°)i2^i2, where 
A~* £ ^comp -g compactly supported. Then 

{AEh,Eh) = {A-\tEh,XtEh)+Ot{h'^)\\Eh\\l2^j,^y (5.5) 

We will now write 

A-' = Aq-* + A^\ V - ^"V, ^r* ■■= - ^), (5.6) 

where the norm of the principal symbol of Aq* will decay with t and the operator A^^ 
will be negligible on Ej^. The function ip £ C^{M) is taken independent of t and such that 
suppx C {x > e^} for some e^ and ip = 1 near {x > e^}. We also require that if = 1 near 
{x > El}, where Ei comes from the assumption (A7). 

We first show that the terms in (5.5) featuring both A^* and E^^ are 0{h°°). For that, 
we need to show that the trajectories in WF/i(A^*) C supp(l — cp) H 5*(suppx) satisfy the 
geometric property shown on Figure 1: 

Lemma 5.2. Let t > 0. Assume that {ra,v) G S*M satisfies m G supp(l — (p), hut 
(^~*(m, z^) G supp x- Then: 

(1) {m,L') directly escapes in the forward direction, in the sense of Definition .{.2; 

(2) for each s > 0, g^{m, v) does not lie in the set defined in (4.9), for any ^ G dM . 

Proof. (1) We have x{m) < ei < Eq; therefore, if (m,z/) does not directly escape in the 
forward direction, then it directly escapes in the backward direction; this would imply that 
x{g~^(m, u)) is decreasing in t > 0, which is impossible as x(m) < < x{g~^(m, z^)). 

(2) The point g^(m,i^) directly escapes in the forward direction and x{g^{m,u)) < ei. 
If g^{m,v) G VFg, then by (G4), ^ = ^+oo('Ti, z/), but this is impossible as xo = 1 near 
{g'{m,u),i+^{m,v))hY {Al). □ 

Combining Lemma 5.2 with the microfocal information we have on E\, we get 

Proposition 5.3. If Eh = Xo^h ~^ decomposition (4.6), then for each t>0, 

{AEh,Eh) = {A^\tXoElxtXoEh) 
+{Ao\tEh, XtEh) + Ot{h^{l + WEhWl.^K^))). ^^'^^ 
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where A^^,A^^ are defined in (5.6). 

Proof. By (5.5), it is enough to show that 

{A^\tEh,XtEh) - {A^'xtXoElxtXoEl) = Ot{h'^{l + ||i?h||i2(x^))). 

Given that ||xo-E'^||l2(^j) = 0(1), it suffices to prove 

\\BxtElh2 = Otih'^il + \\E4L2^K^))), 

where B is equal to either A^'' or its adjoint. This in turn fohows from 

\\BxtElh2 = Otih^), (5.8) 

with El defined in (4.8). Take (m,i/) G WFh{BxtEl) C S*M. Then by Proposition 3.8, 

(m, u) G WFh{B) C WFh{A-') D supp(l - ^) C g\WFh{A)) D supp(l - 

Since WFh{A) C suppx, we see that m E supp(l — f) and g^^{m, v) £ suppx; therefore, by 
Lemma 5.2, the geodesic g^{m, v) escapes in the forward direction and does not pass through 
for s > 0. But then by (A6) the point {m,v) cannot lie in WF/i(£'^), a contradiction. 
We showed that the wavefront set of BxtE\ is empty, which implies (5.8). □ 

We now use the averaged estimate (4.16) and the Hilbert-Schmidt norm estimates from 
Section 3.3, to estimate the second term on the right-hand side of (5.7): 

Proposition 5.4. There exists a constant C independent oft such that 

h-^\\{A-^'xtEh,XtEh)\\Ll^^aM^li,i+h]) < C\\{a{A)og-')ip\\L2(^s',j^ + Ot{h). (5.9) 

Here 1^2 (^g* j^^^ is the norm of the restriction of a to S*M with respect to the Liouville 
measure. 

Proof. Take a real- valued function ipi G C^{M) independent of t such that ipi = 1 near 
suppc^. Then the left-hand side of (5.9) is bounded by 

h-^\\{Ao\tEh,^iXtEh)\\Ll^+ h-^\\{il - ipi)A^\tEh,XtEh)\\L^^^^^ 

where the L^, and later L^, norms in ^, A are taken over dM x [1, 1 + /i]. The second term 
here is Ot{h°°) by the bound (4.17) and since (1 — v9i)j4q* = Ot(/i°°)^2_^/,2 is compactly 
supported. The first term can be estimated by applying the Cauchy-Schwarz inequality first 
in m and then in (A,^): 

h~^\\{Ao\tEh,fiXtEh)\\Ll^ < W^o^XtEhh^M) ■ \\fiXtEh\\L^M)\\Ll^ 

< h'^/^\\Ao\tEh\\L2 ^ ■ h-^/^\\ipiXtEh\\L2 ^ . 

Now, h^^^'^\\(piXtEh\\T2 is bounded (independently of t) uniformly in h by (4.17). As for 

m,^,A 

h~^^'^\\A^^XtEh\\r'2 , we can estimate it using (4.16) by a constant times 

m,^,A 

/l"/2p-*Xtn[i,i+;,]||HS. 
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Note that the operator A^^xt £ ^'^°™p is compactly supported and it is compactly microlo- 
calized independently of t. It then remains to apply (3.24) (to the adjoint of our operator); 
by Proposition 3.8, the principal symbol of A^^xt is given by (cr(A) o g~^)ip. □ 

We now use the dynamical assumption that ^l{K) = 0. The function {(j{A) o g^^)ip 
is supported in a t-independent compact set and bounded uniformly in t. Moreover, it 
converges to zero pointwise on S*M\T^ as t — )• +oo. Therefore, by (5.2) and the dominated 
convergence theorem we have {cr{A) o g~^)ip — t- in Lp'{S*M), as i — t- +oo. It then follows 
from (5.7) together with the bound (4.17) and from (5.9) that 



lim limsup/i ^||(yl£;?,, - (A^ *xtXo£^°, XtXo-E^°>|lLi (aMxfi ml) 



0. 



To prove Theorem 1, it now remains to show that 
lim limsup/i"^ (^j^^xtXo^ft, XtXo-E°> - 

i-s>+oo i^^Q 



a 



S'M 



(A) dfl/: 



L](dM) 



(5.10) 



uniformly in A = 1 + 0{h). We first note that by (4.7) the function 

XtXo^^(A, m) = e^<^«(™)xt (Hxo(m, 06°(1, e, m; 0) + Ot{h)L2 

is a compactly supported Lagrangian distribution associated to the Lagrangian from (4.3). 
Therefore, by Proposition 3.5, we find 

A-^\tXoEl{\ = e'^'f'ixtXob\l, m; 0)a(^^*)(m, dmh{m)) + Ot{h)L2. (5.11) 
Therefore, 

{A-^\tX^ElxtX^El) = [ a{AY'){m,dmM^))\XtXob''{l,tm;0)\^ dvol(m) + Oi(/i). 

Now, by Proposition 3.8, a{A^^) = {(t{A) o g^^){l — if). By Lemma 5.2, this function is 
supported in D<S+n{x < ei}, with 2?<?+ from Definition 4.2. Also, Xt = 1 near supp a{A^^). 
Then by part 2 of Lemma 4.4, 



{A^'xtXoElxtXoEl) 
Therefore, (5.10) reduces to 



{a{A)og-'){l-y,)dfi^ + Ot{h). 



lim 

t— >+oo 



S*M 



{a{A)og '){l-^)dfi^ 



a 



S*M 



(A) dfi^ 



0. 



LUdM) 



(5.12) 



(5.13) 



By part 1 of Lemma 4.4 and (4.13), we write the norm on the left-hand side of (5.13) as 



S*M 



a{A){^og^)dfl^ 



< 



LUdM) JS*M 



\a{A){^og^)\dlJLi 



The expression under the integral on the right-hand side is bounded and compactly supported 
uniformly in t and converges to zero pointwise on 5"*^^ \ r„; by (5.2) and the dominated 
convergence theorem, we get (5.13). This finishes the proof of Theorem 1. 
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The nontrapped case. We briefly discuss the situation when WFh{A) n r_ =0. In this 
case, for t large enough (depending on A), for any {m^v) G WFh{A) we have g^{m,i') 
supp if and thus 

A^' = Oih^)L2^L^. 

Then by (5.7) and the bound (4.17), 

{AEh,Eh) = (>lr*X*Xoi?^,XtXoi?^) + 0(/i°°)ii_^(eMx[i,i+/^])- 

The quantity {A^*xtXoE]^,XtXoE^) is calculated in (5.12) up to 0{h). However, since is 
a Lagrangian distribution, one can get by Proposition 3.5 a full expansion of this quantity 
in powers of h; this yields 

{AEf,{X,0,EhiX,0)= E^'/ ^i««'^€ + 0(^'^''')Li (OMx[i,i+/.])' (5-14) 

0<j<N •^■S'*M 

where A = Op/j(a) for some symbol a and some quantization procedure Op/j and each Lj is 
a differential operator of order 2j on T*M, with Lq = 1. 

5.2. Estimates on the remainder. In this subsection, we prove (1-12) and establish an 
approximation fact (Proposition 5.8) used in the proofs of (1.13) and Theorem 4. 

Classical escape rate and Ehrenfest time. Let Kq C Af be a compact geodesically 
convex set containing a neighborhood of the projection of the trapped set K onto M. As 
in (1.9), define the set 

T{t) = {(m,i/) G S*M \ me Kq, g\m,v) G Kq}. 

The choice of Kq does not matter here: if Kq C M is another set with same properties and 
T'{t) is defined using Kq in place of Kq., then there exists a constant Tq > such that for 
each T >Tq and t > 0, 

g'^{T'{t + 2T))(int). (5.15) 

Indeed, assume that (5.15) were false. Then there exists sequences Tj — )• +oo, tj > 0, and 
{mj,Uj) ^ S*M such that g~'^i {mj,Vj) and g^^^'^^ {mj,Uj) both lie in Kq, but for each j, 
either (1) {mj,Vj) Kq or (2) g^i {mj,Vj) Kq. We may assume that case (1) holds for all 
j; case (2) is handled similarly, reversing the direction of the flow and taking g*^(mj, Vj) in 
place of {mj,Uj). Take e > such that Kq C {x > e}; since {x > e} is geodesically convex 
for £ small enough, we have {rrij^Vj) G {x > e}. Passing to a subsequence, we can assume 
that (mj, Vj) — )• (m, u) G 5"*^^ as j — )■ +oo. Now, since g~'^^{mj, Vj) G Kq and Tj — )• +oo, we 
have (m, u) G r+ (indeed, otherwise there would exist s > such that g~^{m, i^) £ {x < e} 
and this would also hold in a neighborhood of {m,u)). Similarly, since g^^'^'^J {mj,Uj) G Kq 
and tj + Tj — t- +oo, we have (m, v) G r_. It follows that (m, z^) G -fC, which is impossible, as 
each {nij, Vj) does not lie in Kq, which contains a neighborhood of K. 

By changing Aq slightly and using (5.15), we see that the choice of Kq does not matter 
for the validity of (1.12) and (1.13); more precisely, if Aq > Aq, then r'(/i, Ag) < Cr(/i, Aq), 
where r' is defined by (1.10) using T' in place of T- Also, the maximal expansion rate A^ax 
defined in (1.11) does not depend on the choice of Kq. 
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We now choose a geodesically convex Kq such that its interior contains the supports of 
all cutoff functions and compactly supported operators used in the argument below. We will 
rely on Proposition 3.9 (with U equal to the interior of Kq); we let Aq > Amax and fix Eg > 
and A'o such that Aq > Aq > (1 + 2ee)A 

max- Define the Ehrenfest time 
:= log(l//i)/(2Ao). (5.16) 
Then when propagating an operator in ^'^°™p microlocalized inside 

£e, := {1 - Se < W\g < 1 + Ee} (5.17) 

with cutoffs supported inside Kq, as in Proposition 3.9, for time t = Hq £ [—te,te], we get a 
mildly exotic pseudodifferential operator in ^'^J^^, where 

Pe := teA'o/ log(l//i) = A[,/(2Ao) < 1/2. (5.18) 

First decomposition of {AE^, Eh). By Proposition 5.1, we may assume that A G \I^'=°™p(M) 
is compactly supported and microlocalized inside the set See defined in (5.17). 

We first establish the following decomposition similar to (5.7): 

{AEh,Eh) = e''^{A{^U{to)y^Eh,Eh) 

+ J2^''^{A{^U{to)y{l - ^)ipt,xoElEh) + 0{h^M{Enf), ^^'^^^ 

uniformly in ^ G dM and AG [1,1 + hi]. Here / = C'(log(l//i)) is a nonnegative integer and 
to > and (/?,(/?io G C^(M), specified below, are independent of j. The quantity Af{Eh), 
defined in (5.22), is related to the norm of Eh on a certain compact set, and is bounded 
on average by (5.23). The real- valued parameter /3 is equal to 

(3 = -to{X'^ + coh^)/{2h) (5.20) 

and will not play a big role in our argument. 

To show (5.19), we start by considering the functions ip,ipi,ip2 G C^(M) such that: 

• < (/?2 < 1 everywhere, 

• Lf = 1 near supp(^2 and ipi = 1 near suppc^, and 

• ip2 = I both near the support of A and near the set {x > ei}, with ei defined in (A7). 

The proof of (5.19) only uses the function ip, however the other two functions will be required 
for the more precise decomposition (5.35) below. 

We now have the following analogue of Lemma 5.2: 
Lemma 5.5. There exists to >0 such that if (m, v) G S*M satisfies 

m G supp(l — 992) and g^*{m,v) G supp<y9i for some t > to, (5-21) 

then: 

(1) (m, z/) directly escapes in the forward direction; 

(2) for each s > 0, g^{m,v) does not lie in the set defined in (4.9) for any ^ G dM; 
and 
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\ g *(m, i^) \ 




Figure 3. An illustration of Lemma 5.5. The functions (p,ipi,{p2 are sup- 
ported to the left of the corresponding dashed lines; the right side of the figure 
represents infinity. 

(3) for each s > to, g^{m,v) supp9?i. 

Proof. (1) Let supp(/?i d {x > e^}. The set V8^ n {x > e^}, where V8- is specified in 
Definition 4.2, is compact; therefore, there exists to > such that for t > to and {m,v) G 
V£^ n {x > e<^}, we have g~*(m, z^) supp(/Ji. 

Now, assume that (m, z^) satisfies (5.21), but it does not directly escape in the forward 
direction. Since {m,i>) E supp(l — ip2), we have x{in) < eq; therefore, {m,v) G D<S_. Then 
x{m) > x(g^^(m, u)) > Eip] therefore, (m, v) G VS- n{x > e^}, a contradiction with the fact 
that g~^{m,u) £ supp(/9i and t > to- 

(2) This is proved exactly as part 2 of Lemma 5.2. 

(3) It is enough to use part (1), take to large enough, and use that the set H {x > e,^} 
is compact. □ 

Take to from Lemma 5.5. Let ipt^ G C^(M) be real-valued and satisfy cig(supp ipi, supp(l — 
fto)) > ^0- Take a compact set Ktg C M whose interior contains suppc/^jg. Put 

M{Eh):=l + \\Eh\\L2iK,,y, (5.22) 
this quantity depends on A and ^ and we know by (4.17) that 

^"'ll-^(^'^)ll4,(aMx[i,m]) = 0(l)- (5-23) 

By (1.4) and Lemma 3.10, we have similarly to (5.4), 

ipEh = e'^ipU{to)iptoEh + 0{h'^M{Eh))L^. (5.24) 

Here (3 is given by (5.20). Iterating (5.24) by writing ipt^ = (p + {1 — ^)fto, we get for 
/ = 0{log{l/h)) (or even for / polynomially bounded in h) 

I 

^Eh = e''^{^U{to)yfEh + Y,e'"'i^U{to)yil - ^)^toEh + 0{h^M{Eh))L^, (5.25) 
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uniformly in ^ E dM and A G [1,1 + /i]. Same is true if (p is replaced by any function 
If' G C^{M) such that (ig(supp 93', supp(l — ifto)) > to. One can also replace U{to) by 
[/(-to). 

We now use our knowledge of the wavefront set of Ej^ to prove the following analogue of 
Proposition 5.3: 

Proposition 5.6. If = Xo^^h ~^ decomposition (4.6), then 

\WU{to){l - f)ftoEl\\L^ = 0{h^Af{Eh)), (5.26) 

uniformly in S, ^ dM and AG [1, 1 + /i]. Same is true if we replace each instance of (p by any 
function in the set {ip, ipi, 922} ■ 

Proof. Recalling the definition (4.8) of Ej^, we see that (5.26) follows from 

\\fU{to){l - ip)iptoElh^ = 0{h^). (5.27) 

We now make the following observation: a point (m, v) G S*M in the wavefront set of E\ 
will make an 0{h'^) contribution to (5.27) unless m G supp(l — ip), but g~^°{m, v) G supp ip; 
however, by (A6) and Lemma 5.5, in this case {m,v) WF/i(£'^). To make this argument 
rigorous, we can write (bearing in mind that E^ is polynomially bounded) 

iptoEl = BEl + 0{hnL^, 

where B G '1''^°™? is compactly supported and such that 

(m, v) G WFft,(i?) n supp(l - ip) g~^° (m, v) ^ supp tp. 

Then the operator (^[/(to)(l — f)B is 0{h^)]^2^i2 by part 2 of Proposition 3.8, which 
proves (5.27). □ 

Using (5.26), we can replace E^ by XqE'^ in each term of the sum (5.25): 



I 

ipEh = e''^{ipU{to)yipEh + Y,e''H^U{to)y{l - ip)ft,XoEl 

i=i 

+0{h'^MiEh))L2. 



(5.28) 



Applying the operator A = ipAip, we get (5.19). 

Properties of propagators up to Ehrenfest time. We will now establish certain prop- 
erties of the cut off and iterated propagators up to the Ehrenfest time te defined in (5.16), or, 
in certain cases, up to twice the Ehrenfest time. The need for these properties arises mostly 
because of the cutoffs present in the argument. Define the Ehrenfest index 

le := [te/to\ + 1 ~ log(l//i). (5.29) 

Lemma 5.7. Assume that ip',(p" G Cq°(M) satisfy \ip'\, \ip"\ < 1 everywhere. Let B G ^''^°™p 
be compactly supported and microlocalized inside the set E^^ defined in (5.17). Then: 
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1. If if" = 1 near supp , then for < j < le, 

iip'ui±to)yBuiTjto) = {^'u{±to)yBiuiTtoWy + oih°^)L2_L2, (5.30) 

iU{±to)^yBU{Tjto) = {U{±toWYB{ip"U{^tQ)Y + Oih^)L2^L2. (5.31) 

2. If Bi,B2 G vI/'^°™P satisfy same conditions as B and moreover W¥ h{Bi)r\W¥ h{B2) = 0, 
then for < j < 21^ (that is, up to twice the Ehrenfest time) 

Bi{^'U{±to)yB{U{Tto)v"yB2 = 0{h^)L2_L2. (5.32) 

Same is true if we replace ip'U{±to) by U{±to)ip' and/or replace U{'^to)ip" by ip"U{'^tQ). 

3. If ip" = 1 near supp (p' and both if" and B are supported at distance more than to from 
supp(l — ifto), then for < j < le (P is defined in (5.20) ) 

e^'^^{ip'Ui±to)yBEh = iip'U{±to)yBiU{Tto)^"yEh + 0{h^M{Eh))L2, (5.33) 
e^'^''{U{±toWyBEh = {U{±to)^'yB{^"U{Tto)y^toEh + 0{h^M{Eh))L2. (5.34) 

Proof. We will repeatedly use Propositions 3.8 and 3.9 and omit the 0{h°°)]^2_^i2 remainders 
present there. 

1. We prove (5.30); (5.31) is proved similarly. Assume that the signs are chosen so 
that (5.30) features ip'U{to). We argue by induction in j. The case j = is obvious. Now, 
assume that (5.30) is true for j — 1 in place of j. Then 

{ip'U{to)yBU{-jto) = v'B' + 0{h^)L2^L^, 

where 

B' = U{to){v'U{to)y''B{U{-to)v"y-'U{-to) 

is a compactly supported operator in ^'p°™^, with pe defined in (5.18). Since supp(/5' n 
supp(l — (p") = 0, we have 

ip'B' = if'B'ip" + Oih°-)L2^L2 

and (5.30) follows. 

2. We again assume that the signs are chosen so that (5.32) features ip'U{tQ). Write j = 
ii+j2i where < ji, j2 < le, and write the left-hand side of (5.32) as U{jito)BiBB2U{—jito), 
where 

B = {^'U{to)y'B{U{-to)ip'T, 
Bi = U{-jito)Bi{ip'U{to)y\ B2 = iU{-toWy'B2U{jito). 

Now, i? is a compactly supported member of ^p°™^. Same can be said about Bi and B2, 
by applying (5.31) and its adjoint (where the role of (p' is played by either 93' or ip" and 
the role of ip" , by a suitably chosen cutoff function). Moreover, if Ui,U2 are bounded open 
subsets of T*M such that WF^iBk) C Uk and C/i n C/2 = 0, then by Proposition 3.9, B^ is 
microsupported, in the sense of Definition 3.1, on the set g-'^^^{Uk)', since these two sets do 
not intersect, we see that B1BB2 = 0{h°^)i2_^i2 as needed. 



42 



SEMYON DYATLOV AND COLIN GUILLARMOU 



3. We once again fix the sign so that C/(io) stands next to ^p' . Formally, (5.33) and (5.34) 
follow by applying (5.30) and (5.31), respectively, to the identity e^^^E^ = U{—jto)Eh. To 
make this observation rigorous, we write by Lemma 3.10 

e'^BEh = BU{-to)iptoEh + 0{h'^M{Eh))L2, 

e'^^"Eh = ip"U{-to)vtoEh + 0{h^M{Eh))L^. 

We now use induction in j. For j = 0, both (5.33) and (5.34) are trivial. Now, assume that 
they both hold for j — 1 in place of j. We then write 

e'^P{^'U{t^)yBE^, 

= e'('{^'U{to)yB{U{-to)^"y-'EH + 0{h^N{En))L^ 

= {ip'U{to)yBiU{-to)ip'y-'Ui-to)ipt,Eh + 0{h^MiEh))L2. 

The operator {ip'U{to)y B{U{—tQ)if"y^^U{—to) is a compactly supported element of ^fp"™^; 
moreover, as j > 1, the wavefront set of this operator is contained in suppc/?'. Since (p" = 1 
near suppc/?', we can replace iptQ by if" in the last formula, proving (5.33). 

We next write 

e'^^{U{t^)f'yBEh = e'P{U{to)f'yB{f"U{-t^)y-\t,Eh + 0{h^N{Eh))L^. 

However, (p' {U{to)ip'y^^ B{ip"U{—to)y~^ is a compactly supported element of ^™™p and its 
wavefront set is contained in suppc^'. Since ip" = 1 near suppip', we can replace ipt(^ by ip" , 
obtaining (5.34): 

e'^f'{Uito)ip'yBEh 
= e'^iU{to)f'yB{f"U{-to)y-'v"Ef, + 0{h^Af{Eh))L2 
= {U{to)v'yB{^"U{-to)y-'ip"U{-to)ftoEh + 0{h^N{Eh))L^. □ 

Second decomposition of {AEf^, Ef^) . We now analyse the terms of (5.19), reducing 
{AEfi, Eh) to an expression depending on the 'outgoing' part E'^ of the plane wave (see (4.6)), 
with remainder estimated by the classical escape rate for up to twice the Ehrenfest time: 

Proposition 5.8. For < I < 2le, 

(Aii;,, ii;,) = ^^(I^xoi^o, xoi^o) + o(Vz.(r(/to)) + ,,MxM+/^])' 

j=i (5.35) 

:= ipt,il - ip2)iUi-to)^iyAifUito)y{l - f)ipt,. 

Here 1^ is defined in (5.29) and T{t) in (1.9); we keep the 0{h°°) remainder to include the 
nontrapping case. 

We will use Lemma 5.7; since it only applies to pseudodifferential operators microlocalized 
inside the set S^^ from (5.17), we take an operator 

Xo £ ^^°'"P(M), WF;,(Xo) C fee, ^0 = 1 near S*M n supp v^i^, (5.36) 
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compactly supported inside Kt^. By (1.4) and the elliptic estimate (Proposition 3.2), we 
have 

ipt,Eh = Xo^t,Eh + 0{h^M{Eh))L2 = ^toXoEh + 0{h^MiEh))L^. (5.37) 
Same is true if we replace E^ by Xo^h^ as by (A4) and the fact that \dm4'^\g = 1, we have 
WFhixoK) C S*M. We also recah that WFtiA) C £e,- 

We start the proof of Proposition 5.8 by estimating the first term on the right-hand side 
of (5.19) for I up to twice the Ehrenfest time, in terms of the classical escape rate: 

Proposition 5.9. There exists a constant C such that for < / < 21^, we have 

^-i(^(^C/(to))Vi?/.,i?/.)llL|_,(9Mx[i,m]) <C^M^(n^*o)) + 0(/i'"). (5.38) 

Proof. We write I = h + I2, where < /i, ^2 < ^e! then 

{A{ipU{to)yipEh,Eh) = {{ipU{to)y'^Eh,{U{-to)ipy^A*Eh). 

Now, by (5.33) 

e''^^i^Uito)y^ipEh = e''^^{^U{to)y^XoipEf, + O(/i°°AA(i?0)L2 
= BlEh + 0{h^MiEh))L2, 

where 

Bi = {ipu{to)y'Xoip{ui-to)^^y\ 

Similarly, by (5.34) (recalling that A = ipAip) 

e-''^^{U{-to)ipy^A*Eh = BfEh + Oih°°MiEf,))L2, 

where 

Bf = iu{-to)ipy'A*{ip^u{to)y'^to- 

Put Bi = {BfyBf; recalling (5.23), it is then enough to show that 

h-^\mEh,Eh)h.^^^9jj,[,^,^^) < Ci,L{T{lto)) + 0{h^). (5.39) 

Now, by Proposition 3.9, the operator Bf is a compactly supported element of ^'p°™^ (modulo 
an 0{h°°)i2^i^2 remainder which we will omit), and it is microsupported, in the sense of 
Definition 3.1, inside the set 5'~'^*"({95i 7^ 0}) (here we only use that suppi^ C {ipi 7^ 
0}). Similarly, Bf G ^'p°™^ is microsupported inside g^'^^°{{^i ^ 0}). Therefore, Bi is 
microsupported on the set 

5/=5~'^*°({^i/0})n5'^*o({<^i/0}). 
Note also that Bi is compactly supported independently of /. 

Now, by taking the convolution of the indicator function of an hf" sized neighborhood of 
5; with an appropriately rescaled cutoff function, we can construct a compactly supported 
operator Bi G such that Bi = BfBi + C'(/i°°)^-cx) and Bi is microsupported inside an 

OihP") sized neighborhood Si of 5/. Using (5.16), (5.18), and the estimate on the Lipschitz 
constant of the flow given by (3.17), we see that for (m, i>) G 5; n S*M, there exists (m, u) £ 
Si n S*M such that d((m, D), (m, i^)) < ChP^ and for Aq > Aq > (1 + 2ee)Amax, 

d(5'i*«(m,i>),5'i*«(m,z^)) < Ce'^*«^o/i''- < Ce^'^^ohP'^ < Ce'^^^^o-^o) 



44 



SEMYON DYATLOV AND COLIN GUILLARMOU 



is bounded by some positive power of h. Here d denotes some smooth distance function on 
T*M. Same is true if we replace 5'^*" with (^~'2*oj therefore, if the compact set Kq used in 
the definition (1.9) of T{t) is chosen large enough, we have 

SinS*M cg'^'^nito)). (5.40) 

Using the Cauchy-Schwartz inequality and (4.16), we bound the left-hand side of (5.39) by 

h-'\mEh,Eh)hr < h-'\\{BiEh,BiEh)hi + 0{h^) 

< h'^\\BiEh\\L2(^M)Ll^ • \\BlEh\\L2(M)Ll^ + ©(/l^^) 
<C(/l"/2||i?,n[i,i+^]||Hs)(/i"/'||S;n[i,i+^]||Hs) + 0(/i'"). 

It remains to use (3.25) (or rather its adjoint). Indeed, both Bi and Bi are bounded in 
uniformly in /, and they are microsupported in Sf, therefore, by (5.40) 

h-'\\{BiEh,Eh)\\Li <c^lL{Slr^s*M) + o{h^)<c^lL{nltQ)) + o{h°^). □ 

^^■^ 

As for the sum in (5.19), we have the following 
Proposition 5.10. For 1 < j < 21^, we have 

e'^^{A{ipU{to)y(,l - ^)^u,XoElEh) = {A^xAxoED + 0{h^N{Ehf), (5.41) 
uniformly in ^ £ dM and A G [1, 1 + /i], with A' defined in (5.35). 

Proof. Since A = (pAip, we can replace E^ by (/'i-E'/i on the left-hand side of (5.41). Writing 
down (5.25) for ipi in place of if and using ip2 in place of ipi in the splitting ipt^ = 921 + (1 — 
fi)^Pto in the last step, we get 

cpiEh = e'^^{^iU{to)yip2Eh + e^^'^((^if/(to))^'(l - V2)y^toEh 

i-i (5 42) 

+ ^e"'f'{^iU{to))Hl - Vi)^toEh + 0{h^N{Eh))L^. 

k=l 

We now substitute (5.42) into the left-hand side of (5.41). The first term gives, after us- 
ing (5.37) to replace (ptoXoEl by Xo^ftoXoE'^^ and ipEh by XoipEh 

{A{^U{to)y{l - ^)n,XoEl, {ipiU{to)y^2Eh) = (BoXoElEh) + 0{h°^M{Eh?), 

where 

So = ^2x*^{u{-tQ)^iyA{^u{tQ)y{i - ^)Xono = o(^~)l2^l2 

by (5.32), as suppv92 nsupp(l — (/?) = 0. 

Next, we use Proposition 5.6 to write the second term of (5.42) as 

e'^^{^iU{t^)y{l - ^2)VtoXoEl + 0{h^M{Ef,))L2; 

therefore, this term gives the right-hand side of (5.41). 

It remains to estimate the contribution of each term of the sum in (5.42), which we can 
write, using (5.37), as e^(j-'=)^(BfcXo^°, Xo^°) + 0(/i°°AA(S/,)2), with 

Bk = iptoX^il - ipi){Ui-to)ipi)''A{ipU{to)yil - ^)Xo^to. 
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We need to show that ||^fc||_L2^^2 = 0{h°°) for 1 < k < j . For that, we consider two cases. 
First, assume that k < Ig. Then we have 

^toX^il - ipi)iUi-to)vi)'AiipUito)fip = 0{h^)L2^L2, (5.43) 

as the supports of 1 — and if do not intersect and the operator in between them is a 
compactly supported element of '^^p™^ (modulo an 0{h°°)i2^]^2 remainder which we will 
omit). Since Sfc is obtained by multiplying the left-hand side of (5.43) on the right by 

U{to)ivU{to)y-'-''{l - ^)Xoipt„ it is also 0{h^)L2^L'^. 

Now, assume that k > le- Take G Cg°(M) equal to 1 near suppyji and such that 
|^i| < 1 everywhere. We write by (5.30) and its adjoint, 

U{{k - Qto)BuU{-{j - IM = BlBlBl + 0(/i~)^2^i2, 

Bl = {^,U{to))'-'^^toX*o{l - (^i)(C/(-to)<^i)'-'% 

Bl = iUi-to)^iy^A{^U{to)y% 

Bl = i^U{to)y-'^{l - ^)Xoipt,iU{-to)^iy-'^. 

Now all Bl, i = 1,2,3, are compactly supported members of 'f™™'^. Let Ui,U2 be two 
bounded open sets such that supp((/3t(j(l — fi)) C Ui, suppc^ C U2, and Ui Ci U2 = 0. 
Since k — le > j — le and by Proposition 3.9, the operator Bj, is microsupported, in the 
sense of Definition 3.1, on the set g~^^~^'^'>^^{Ui), while -B| is microsupported on the set 
g-{k-^e)to(jj^y^ since these two sets do not intersect, we get B^ = 0{h'^) 1^2^12, finishing the 
proof. □ 

Combining (5.19) with (5.38) and (5.41), we finally get (5.35). 

Proof of (1.12). Put / equal to the number le defined in (5.29). By (5.35), it is enough to 
approximate the terms {A^xo^h-' Xo^h)- This is done by the following proposition, relying on 
the Lagrangian structure of and featuring the interpolated escape rate r{h, A) from (1.10): 

Proposition 5.11. Put I = le given by (5.29), and r(h,A) defined in (1.10). Then the sum 
on the right-hand side of (5.35) is approximated as follows: 

^{A^XoElxoEl)= [ a{A)dfi^ + 0{hr{h,2Ao))^i.gjj^,,,^^y (5.44) 

~[ JS*M ^ " 

Proof. By Proposition 3.9, the operator A^ is compactly supported and lies in modulo 
an 0{h°°)i^2^i^2 remainder, where 

Pj = y^Pe, with jto/te < 1 + 0(1) 

with te and pe defined in (5.16) and (5.18), respectively. Next, A^ is microsupported, in the 
sense of Definition 3.1, in the set 

:=<7*«({9^i^0})nf/^*»({(/.i^0}). 
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If the set Ko from the definition (1.9) of T{t) is large enough, then QjnS*M C g^^°{T{jto)); 
by the definition (1.10) of r{h,A), we find 

/ji-itoAe^^(Q^. p < ^(/.^ 2Ao). (5.45) 

By (A4) and Proposition 3.5, we have the following analogue of (5.11): 

Moreover, by part 2 of the same proposition, we see that A^xoK is 0(/i°°) outside of the 
set of points m £ such that {m,dm(t>£,{fn)) £ Qj- By Lemma 5.5, (y{A^) is supported in 
supp(l — </?) n 5*°(supp(^i) C V£j^ n {x < ffi}, with V£j^ from Definition 4.2. Using part 2 
of Lemma 4.4, we then get 

{A^xA XoEl) = [ a{A^) dfi^ + 0{h^~^P^ ^i^{Q^)) + 0{h^), (5.46) 

JS*M 



uniformly in ^ € dM and A € [1, 1 + /i]. Now, we write by (5.45) and Proposition 4.3, 

I I 

r 1 

< r{h, 2Ao) E /i^^-^^^)^*"/*'^ = r{h, 2Ao) E 6-^^"^^-'^^''^^^° < Cr(/i, 2Ao). 
3=1 j=i 

It remains to sum up the integrals in (5.46). We have by Proposition 3.9, bearing in mind 
that ififi = if, {1 — (p){l — ip2) = I — if, A = (pAif, and dg(supp(l — (^jg), supp (^9) > to, 

i-i 

a{A^) = {a{A) o 5-i*o)(i - <^) J] o 5-^*0. 

k=l 

By part 1 of Lemma 4.4, we write 

E / dfi^ = [ cr{A) ^{l-^o 5^*0) JJ ^ o /*o d/i^ 

JS'M 

It remains to note that by Proposition 4.3, 

/_/ \(y{A)\\\^og^'^dl,^di= / \a{A)\\\^og^'^d^iL = 0{^JiL{T{lto))) 
JdM Js*M JS'M 

since the expression under the last integral is supported in T{ltQ). □ 
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5.3. Trace estimates. In this subsection, we prove a stronger remainder bound (1-13) for 
the case when {AEh, E^) is paired with a test function in ^ and obtain an expansion of the 
trace of spectral projectors with a fractal remainder - Theorem 4. 

Expressing via Schrodinger propagators. Our argument will be based on 

the decomposition (-5.35). The remainder in this decomposition is already controlled by 
the escape rate at twice the Ehrenfest time tg, defined in (5.16). However, in the previous 
subsection (see Proposition 5.11), we were only able to estimate the sum in (5.35) for I up to 
the Ehrenfest index /g ~ te/to defined in (5.29). We therefore need a better way of writing 
down the Lagrangian states when coupled with a test function in ^, and such a way is 
provided by 

Lemma 5.12. Let /(^) G C°°{dM) and define for A G (1/2,2), 

n^(A) := / /(e)(xoi?^(A,e)) ^ ixoEliX,0)dC- (5-47) 
JdM 

Here denotes the Hilhert tensor product, see (2.8). Assume that Xi,X2 G ^'^°™p(M) 
are compactly supported and the projections Trs(WFh{Xj)) ofWFh{Xj) onto S*M along the 
radial rays in the fibers ofT*M lie inside T>£j^ f^{x < £i}, with VEj^ defined in (4.2) and ei 
from (AT). Then 

XiM(A)X2* = {2TThT r e-'^'-'/^^''^Uis)Bs{X)ds + 0{h°°)L2^L^, 

J-To 

where To > is independent of h and Bs{\) G \i/'^°'^P(^M) is compactly supported on M, 
smooth and compactly supported in s G {—Tq,Tq), and smooth in A. Moreover, i/^+oo is the 
function defined in (G3), then 

a{Bo{l))\s'M = fi^+ooMXiX*^). (5.48) 

Proof. We write 

XiU}{X)X*= [ f{0{XiXoE'i{X,0)^iX2XoEl{X,0)d^- 

JdM 

By (A4), xoE^{X,^) is a Lagrangian distribution associated to A times the Lagrangian 
from (4.3). By Proposition 3.5, we can write 

X,xoEt{X,Oim) = eX'^eM6^.(A,e,m;/i) + 0(/i-)c-, 

where i;^^ is defined in (G4) and bj is a classical symbol in h smooth in A, ^, m and com- 
pactly supported in m. The symbol bj depends on the operator Xj; in fact, we can make 
suppfoj C T~^{TTs(yVFh{Xj))), with r defined in (4.2). We then write the Schwartz kernel of 
Xinj(A)X|, modulo an 0{h°°)c^ remainder, as 

n(m, m'; A, h) = [_ e^('^«(™)-'^«(™')V(Ofei(A, m; h)b2{X,^;m' ,h) (5.49) 
JdM 

Now, the support of each bj in the (m,^) variables lies in the set defined in (G4). The 
critical points of the phase X{(p^{m) — 0^(m')) are given by d^(j)^{m) = d^(t)^{m'); using (G6), 



48 



SEMYON DYATLOV AND COLIN GUILLARMOU 



we see that /i-"/2n( m, m'; A, h) is a Lagrangian distribution associated to the Lagrangian 

Ax:={{m,iy;m',u')\\u\g = X, 3s e {-To,To) : g%m,iy) = {m' (5.50) 
Here Tq > is large, but fixed. 

Now, take some family Bs{X) G x]/'^°™p(M) smooth and compactly supported in s G 
(—To, To) and define the operator 

UeiX) := (2^/i)" e-'^''/^^'''^U{s)Bs{X) ds. (5.51) 

Following the proof of Lemma 3.11, we see that /i^"/^ times the Schwartz kernel nB(m, m'; A, h) 
of nB(A) is, up to an 0{h°°)i2^i2 remainder, a compactly supported and compactly microlo- 
calized Lagrangian distribution associated to the Lagrangian A^. Moreover, the principal 
symbol of /i~"/-^nB(m, m'; A, h) at (m, u, m' , v') such that g^ijn, u) = (m', v') is a nonvanish- 
ing factor times a{Bs){m' ,u'). Arguing as in the proof of part 2 of Proposition 3.3, we see 
that we can find a family of operators Bs{X) such that 

n(m,m';A,/i) = IlB{m,rn -X^h) + 0{h^)c^ . 

It remains to check that the family Bs{X) can be chosen to depend smoothly on A uniformly 
in h (this is not automatic, as multiplication by e^'^^'^'^ for some function destroys this 
property, but does not change the Lagrangians where our kernels are microlocalized for each 
A). For that, it is enough to note (by Proposition 3.3) that if we consider /i~"/^n and 
h-"^'^UB as La grangian distributions in m,m',X, they are associated to the Lagrangian 

{{m,u,m' ,iy',X,qx) \ \u\g = X, 3s £ {-Tq,Tq) : g'' {m,v) = {m! ,u'), qx = -Xs}, 

where q\ is the momentum corresponding to A. For H, this is true as when T{ra',£^) = 
g^(T{m,(,)), we have (j)^{m) — (f)^{m') = —s by (4.4); for H^, this is seen directly from the 
parametrization (3.27), keeping in mind the factor e-«-^^*/(2^) the definition of 11^. 

Finally, to show the formula (5.48), put A = 1, take an arbitrary Z £ ^i'^°^p^ and compute 
the trace 

Tr{XiU}il)X^Z) = (2^/i)" e-^^/(2/^) Tr(C/(s)S,(l)Z) ds + 0(/i°^). (5.52) 

J —To 

The left-hand side of (5.52) can be computed as at the end of Section 5.1, using the limiting 
measure /j,^; by Proposition 4.3, it is equal to the integral of /(^+oo)c(X|ZXi) over the 
Liouville measure on S*M, plus an 0{h) remainder. The right-hand side of (5.52) can be 
computed by the trace formula (3.32), and is equal to the integral of Bq{1)Z over the Liouville 
measure on S*M, plus an 0{h) remainder. Therefore, 

/ a{Z)f{^+^)a{XiX*2)dfiL= [ a{Z)a{Bo{l))dfiL 

JS*M JS*M 
for any Z; (5.48) follows. □ 

Proof of (1.13). By (5.35), it is enough to approximate the sum in this formula up to twice 
the Ehrenfest time: 
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Proposition 5.13. Fix f G C°^((9M). Put I = 2lf>, where le is defined in (5.29), and 
consider 

SfW ■= E [_f{0{A'xoE'i{X,0,XoEliX,0)dC- (5.53) 

~[JdM 

If i+ooij^ii^) is the limit of g^{m,i') as t ^ +co, for {m,u) G 5*M\r_ (see (G3)), then for 
X G + 

Sf{X)= f f{C+ooMA)dfiL + Oir{h,Ao)). (5.54) 

JS*M 

Here r(/i, A) is defined in (1.10). Moreover, for each k 

sup \d^Sf{X)\ < Ckh-^P% (5.55) 
xeli,i+h] 

where pe is defined in (5.18). 

Proof. First of all, take a compactly supported operator X G \I''^°™P such that WF/j(X)nS'*M 
lies inside the set V£j^ n {x < ei} and for Xq defined in (5.36), 

^u{to){i - ip)^tMi -x) = o{h°^)L2^L2, 

<^lC/(to)(l - 9^2)^*0^0(1 -X) = 0{h^)L2^L^. 

Such an operator exists by Lemma 5.5 (it can be easily seen that in part 1 of this lemma, 
(m, z^) actually lies in the interior of DS^). Then by (5.37), the definition (5.35) of , and 
Lemma 5.12, 

SfW = y2 [_f{0{A'XoXxoE'fliX,0,XxoEl{X,0)d^ + 0{h^) 

■^^ JdM 

= Tt{A^ XoXUj{X)X*) + OQf^) 

^ /"To 

= (27r/i)" ^ / e-*^'^/(2/^) Tx{A^XoU{s)Bs{X)) ds + 0(/i~) 
j=i •'-^^ 

for some fixed Tq > and some family Bs{X) G smooth in s and A and compactly 

supported in s; we can make Bg microlocalized inside the set defined in (5.17). We will 
henceforth ignore the 0{h°°) term. 

Now, take 1 < j <l and put j = ji + j2, where < ii,j2 < le, 32 > 1, and |ji - J2I < 1- 
Using the cyclicity of the trace, we find 

Tv{A^XoU{s)Bs{X)) = TiiU{s)B{BijX)), 

b{ := iUi-to)^iy'AiipUito)y\ 

Bi JX) := iipU{to)ril - ^)ipt,XoU{s)Bs{X)^to{l - ^2){U{-to)mrU{-s). 

Put pj = {jto/te)pe', since ji,j2 < j/2 + 1, by Proposition 3.9 the operator B^ is a compactly 
supported element of ^^"^^^ (modulo an 0{h°°)^2^j^2 remainder which we will ignore). Same 
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can be said about the operator 

Bl^iX) = ((^C/(to)P • (1 - ^)ipt,XoU{s)B,iX)^t,il - V2)U{-s) ■ (C/(-to) • Uis)^iU{-s)y\ 

(The operator U{s)ipiU{—s) is not pseudodifferential because ipi is not compactly microlocal- 
ized, but this can be easily fixed by taking Xq G \j/comp equal to the identity on a sufficiently 
large compact set and replacing U{s)ipiU{—s) by U{s)ipiXQU{—s) in i?2,s(A), with an 0{h°°) 
error.) Therefore, B-IB2 g{X) also lies in moreover, it depends smoothly on s and A, 

uniformly in this operator class. (In principle, we get powers of I ~ log(l//i) when differ- 
entiating in s, due to the {U{—to) ■ U{s)ipiU{s)y^ term, but they can be absorbed into the 
powers of h in the expansion (3.32).) 

We now use the trace formula of Lemma 3.12, writing 

I J' 

Sf{\) = i2nhr E / " e-^^'/(2'^) TriUis)BiBl{X)) ds. (5.56) 



The operator i?2,s(A)-' is microsupported, in the sense of Definition 3.1, inside the set 
c/-J2*o({(^2 / l})n5"(J2-i)to(-{(^^ ^ 0}); by Lemma 5.5, this set lies inside g~^^^''{V£+) and in 
particular does not intersect any closed geodesies, therefore (3.31) holds. The estimate (5.55) 
now follows immediately from (3.32). The power h~'^P'' arises because we integrate over the 
energy surface {{I'lg = A} depending on A; therefore, d^Sf{X) will involve fcth derivatives 
of the full symbol of BIB2 ^(A) in the direction transversal to the energy surface, which are 
bounded by h~^PK The sum (5.56) has / ~ log(l//i) terms; however, our estimate is not 
multiplied by log(l//i) because one can see that the sum of Liouville measures of the sets 
where these terms are microsupported is bounded. 

As for the approximation (5.54), we write (note that we take s = in i?2,s) 

ii-i 

a{B{) = {a{A)og-^-'^)\{^og-'^'\ 

k=l 

h-l 

a{Bl,{X))\s^M = ((1 - ^MBo{X)))og^^^o Yl ^o/*o. 

fc=0 

Since the Liouville measure is invariant under the geodesic flow, the contribution of the 
principal term of (3.32) to Sf{X) for A = 1 is 

. I j-i 

/ a{A) ((1 - ^MBoiX))) o 5^*0 J] ^ o /*o 
J^'^' 3=1 k=i 

Now, by (5.48), a{Bo{l)) = /(^+oo) on the support of the integrated expression; recombining 
the terms as in the proof of Proposition 5.11, we get the right-hand side of (5.54), with an 
fJ-LiTilto)) remainder. The subprincipal terms (and also the difference Sf{X) — Sf{l) for 
A S [1, 1 + /i]) are estimated using the bound on the Liouville measure of the set where 
BIB2 g is microsupported; arguing as in the proof of Proposition 5.11, we see that they are 
bounded by a constant times r(/i, Aq). □ 
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Expansion of the trace of spectral projectors in powers of h. We now use the 

results obtained so far to derive an asymptotic expansion for the trace of the product of 
the spectral projector Ijg ;^2] (P(/i)) with a compactly supported pseudodifferential operator, 
with the remainder depending on the classical escape rate for up to twice the Ehrenfest time. 
Here we denote 

P{h):=h^{A-co), (5.57) 

with the constant cq from (Al). It will also be more convenient for us to use the spectral 
parameter s = in the following corollary and theorem (not to be confused with the time 
variable s used in Lemma 5.12). 

We start with the following consequence of the decomposition (5.35), the bound (5.54), 
and the spectral formula (4.5): 

Corollary 5.14. Take Aq > Amax; with A^ax defined in (1.11), and let T{t) be defined 
in (1.9). For e > 0, let (p ^ C'o°((l - 1 + ^)) (^Qual to 1 near [1 - e/2, 1 + e/2] and for 
s GM, define ips ■= (f ■ l(-oo,s]- -(f ^ > small enough, then for each compactly supported 
A G ^^{M), there exist some functions Sh{s),Qh{s) and some constants C > 0, Ck > such 
that for all s and all k gN 

Tr {Aips{P{h))) = Sh{s) + Qh{s), \d^Sh{s)\ < Ckh-^-'-''/\ 



\Qh{s + u) - Qh{s)\ < Ch-''fiL[rl^^-^ j j +0{h'^) forue[0,h]. 
Proof By (4.5), 

Tr {AipsiPih))) = i27Th)-^-' j A"/n(A//i)(^,(A2) J {AEh{X,0, Eh{X,()) d^dX. 

v^te dM 

Now, note that the decomposition (5.35) (with / = 21^) is still valid in any 0{h) sized interval 
inside (\/l — e, \/l + e), if e is small enough. More precisely, if we write 

Sh{s) := {27rh)-'^-^ J X'' fuiX/h)ips{X^)Si{X) dX, 

where S'i(A) is defined by (5.53) with /(,^) = 1, then we have the expansion (5.58) with Qh{s) 
satisfying the required bound. To estimate the derivatives of Sh{s), it now suffices to use the 
bound (5.54), noting that it is valid for |A^ — 1| < e if e is small enough. □ 



Using the last corollary, we can show the following trace decomposition with a fractal 
remainder: 

Theorem 4. Let P{h) be defined in (5.57), let A = Op;j(a) G '^'^{M) be a compactly 
supported operator, then there exist some smooth differential operators Lj of order 2j on 
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T*M, depending on the quantization procedure Op^, with Lq = 1, such that for any compact 
interval I C (0, oo), all s £ I, all h > small, and all N G N 

N 

Tr(^]l[o,,](P(/i))) ={27rh)-^-^J2'^' I L.adfi^ + h-^0{f,LiTiAo^\logh\)) + h'') 

where the remainder is uniform in s. Here is the standard volume form on T*M ; we have 
fJ-ui = ^s^^/(''^ ~^ ^)'' '"'^fire ojs is the symplectic form. 

Proof. By rescaling h, it suffices to prove the result for |s — 1| < e/2 where e > is obtained in 
Corollary 5.14, we can thus assume |s — 1| < e/2. Let (ps be defined as in Corollary 5.14, and 
■0 G C^((-l + e/2, 1 - e/2)) such that V + = 1 on [0, 1 + e/2]. For s G (1 - e/2, 1 + e/2), 
one has l[o_s] (-P(/i)) = ips{P{h)) + ip{P{h)) and it suffices to study the expansion in h of 
aA,h{s) and Tr{Aip{P{h))) where 

aA,h{s) ■■= Tr {Aips{P{h))) = Tr {Aips{P{h))x) ■ (5.59) 

if X G C^{M) is such that A = X-^X- Since A is compactly supported, one can use the 
functional calculus of Helffer-Sjostrand [DiSj, Chapters 8-9] to deduce that A7p{P{h))x G 
qjcomp^j^-j ^ compactly supported and microsupported pseudodifferential operator with a 
classical symbol'^. Its trace thus has a complete expansion in powers of h (see [DiSj, Th 9.6]): 

N 

TiiA'iP{P{h))x) = (27r/i)-"-i V / L''a dfi^ + 0(/i-'^+^) (5.60) 

where Lj are some differential operators of order 2j and LQa{m, u) = a(m, v)ip{\v\'^g). 

Let us now analyze the function (JA,h- This is a smooth function of s > by the smoothness 
assumption on the Efi(X, ^) in A, it is constant in s for |1 — s| > e, and we know that (JA^his) = 
0{h~'^~^) uniformly in s by Lemma 3.11. Let 9{s) G o?*'(M) be a Schwartz function such that 
6 G Co°(— 7/, ry) for some small r/ > and 9{t) = 1 near i = 0, and let 6h{s) = h~^9{s/h). We 
write 

a'^^l,{s) := dsaA,h{s) = Tr{A^{P{h))dU,{P{h))x) G C^moo)) 

where dIls{P{h)) is the spectral measure of P{h). The operator Aip{P{h))dIls{P{h))x has 
a smooth compactly supported kernel and is trace class. We clearly have cr'j^h^^h G =y (M) 
and by a simple computation, its semi-classical Fourier transform is given by 

[ e-'^'a'Ah*0h{s)ds = Tr(^(^(P(/i))e-*^^W)^(t) 

and thus 

^'A,h*Ohis) = -^ [ e^t*Tr(A^(P(/i))e-i^W)^(t)cit. 



An alternative method in the Euchdean near infinity setting is the functional calculus of HelfTer-Robert 
[HeRo]. 
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Now we can apply Lemma 3.12 with Bu = ^e^'^°^'^^'^9{—u/2)A(f{P{h)); the condition (3.31) is 
satisfied because 6 is supported in a small neighborhood of zero. This shows that, as /i — )• 0, 
we have the expansion (locally uniformly in s) 

for all G N, where 6 is a symbol such that Opf,{b) = lAip{P{h))+0{h'=^), Lj are differential 

' — ' Tl — 1 

operators of order 2j on T*M, smooth in ys, with Lq = s~2~. In particular, one has 
Lo?)(m, -y/sz/) = ^s^^ a{m, ^/sh')ip{s\u\'^). Notice that b is supported in {{m,i>) S T*M \ 
G suppc^} thus Ljb{m, y/sv) is smooth in s G M when {m,v) S S*M. Since (JA,h{s) is 
bounded in s, the convolution (JA,h * (^h{s) is well defined (as an element in L°°(M)) and we 
have for all N and for all s < 2 

(yA,h * Oh{s) = (27r/i)-"-^ (y^^^ / L^Km, Vuv)d^iL{m, u)du + 0{h^^^)\ . 

j^O ./o Js*M ' 

We are going to show that, uniformly in s G M, 

oaA^) - <yA,h * eh{s) = 0{h^) + 0(/i->L(r(Ao 'I log/i|))), (5.61) 
using the decomposition 

<^A,h{s) = Sh{s) + Qh{s) 
defined in (5.58). Since dsSh{s) is a compactly supported symbol we get by integrating by 
parts A'' times 

for alH G M and all S> 1. Thus, taking the Fourier transform we deduce that Sh{s) — S^-*^ 
^/i(s) = 0{h°°) uniformly in s. From (5.58), we obtain by induction that for all s,u 

\u\ X 



\Qh{s + u)- Qh{s)\ < Ch-'^^l + Y j^L(r(Ao ^1 log/iD) + 0{h^) 

then multiplying by 9h{—u) and integrating in n, we obtain (5.61). 
Given (5.61), we have 

N ,.s r. 

o'A,h{s) = {2TTh)~'^~^ 2_,h-' / / Ljb{m^y/uv)d^L{fn,v)du 

Jo Js'M 

+0{h-^fiL{nA,'\ log h\))) + 0(/i-"+^). 

Since the symbol of b is explicitly obtained from a using Moyal product, we can rewrite 
this expression with a instead of b and with some new differential operators with the same 
properties as Lj but supported in {\i^\g G supp^^}; using polar coordinates S*M x on 
T*M, we deduce that there exist some differential operators L'j of order 2j on T*M such 
that 



(5.62) 



Ljb{m,y/uv)d^L{'m,v)du = l L'M{m,v)di.i^{m,v) 
JS'M JmI<s 
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and LQa{m, v) = ip{\iy\^g)a{m^ v). Combining this with (5.62) and (5.60), we obtain the desired 
result where L,- in the statement of the Theorem corresponds now to L'- + L"-. □ 

6. Euclidean near infinity manifolds 

In this section, we assume that (M, g) is a complete Riemannian manifold such that there 
exists a compact set Kq C M such that for 8 := M \ Kq, 

{£,g) is isometric to (R"+^ \ B{0, R), g^nd) 

where R > 0, B{0,R) is the Euclidean ball of center and radius R and ^euci is the Euclidean 
metric. We will check that all the assumptions of Section 4 are satisfied. 

6.1. Geometric assumptions. We let x £ C°°{M) be an everywhere positive function 
equal to x(m) = |m|~^ in £ identified with M"+^ \ B{0,R), and such that x > R~^ in Kq. 
(We take it instead of the function (1 + |mp)^^/^ used in Section 4 for the model case of 
M""*"^, to simplify the calculations and since we no longer need smoothness at zero.) We 
shall use the polar coordinates m = co/x in £, where uj £ S". Assumption (Gl) is satisfied 
by taking the radial compactification of M, i.e. adding the sphere at infinity: the map 

: M"+i \ B{0,R) (0, l/R) x defined by ip{m) = (x(m),x(m) m) is a diffeomorphism 
and the radial compactification of M is obtained by setting M = M U dM where dM := S", 
the smooth structure on M is the same as before on M but we extend it to M by asking 
that ip extends smoothly to the boundary dM and = (0,.^) if ^ G dM = S" (see for 
instance [Mc] for more details). In other words, smooth functions on M are smooth functions 
on M with an asymptotic expansion in integer powers of l/|m| to any order near infinity. 

Assumption (G2) is clearly satisfied for sq := l/2i2 since the trajectories of the geodesic 
flow in X < eo are simply g*{m, v) = (m + tiy, v). A point (m, v) S S*M is directly escaping 
in the forward direction in the sense of Definition 4.2 if and only if x{ra) < eq and m- v > 0. 
Now, (G3) is satisfied with ,^+oo(?tT'; i^) = i' for (m, u) £ DS^. 

For the assumption (G4), we define 

U^ = {x < eo} X dM C M X dM, 

Then t : U^o ^ S*M from (4.2) maps each (m,^) G (M" \ B{0,2R)) x S" to itself as 
an element of 5*(M"' \ B{0,2R)). Assumptions (G4) and (G5) follow immediately. To see 
assumption (G6), we note that for x{m),x{m') < eq and some ^ G S", we have d^4)^{m) = 
d^(j)^{ra') if and only if m — m' is a multiple of ^. 

6.2. Distorted plane waves and analytic assumptions. We recall a few well known facts 
about scattering theory on perturbations of M", we refer to [Mc] for a geometric approach 
and to [MeZw, HaVa] in a more general setting (asymptotically Euclidean case). A plane 
wave for the flat Laplacian on M"+^ is the function, for A G (1/2, 2), 

u[\- m, i) := ce^™-^, C e rri £ M"+\ c G C. (6.1) 
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This is a semiclassical Lagrangian distribution, its oscillating phase has level sets given by 
planes orthogonal to ^. The continuous spectrum of the Laplacian A associated to the metric 
g is the half-line [0,oo). We will take the resolvent of /i^A to be the L^-bounded operator 

R^{\) ■= (h'^A - X^y^ in Im(A) < 0. 

This admits a continuous extension to {A 7^ 0,Im(A) < 0} as a bounded operator from i^comp 
to Lf^^. For A > we call RhiX) the incoming resolvent and Rh{—X) the outgoing resolvent 

The distorted plane wave is defined for ,^ G (see [MeZw, HaVa]) by 

EhiX,^;m) ■.= 2iXh(^Y lim [(x')-"/'e>^i?,(A; m, C/x')], (6.2) 

where Rh{X;m,m') is the Schwartz kernel of Rh{X) and € £. This is a smooth func- 
tion of (m,^) G M X S", and in the case of M = W^~^^ it is given by (6.1) with c = 1 
(see [Mc, Chapter 1]). We shall use the notation Eh{X,(,) for the C°°{M) function de- 
fined by m I— )• Eh{X,^;m) and we notice that (/i^A — X'^)Efi{X,^) = in M. One has 
Eh{X;m,^) = Eh{—X;m,^) since Rh{X)* = R{—X) for A G M, and the decomposition of the 
spectral measure in terms of these functions in given as follows: by Stone's formula, the 
semiclassical spectral measure is given by 

dUh{X) = —(RhiX) - Rh{-X)) dX for A G (0, 00) (6.3) 
vr 

in the sense that F{h^A) = F{X^)dIlh{X) for any bounded function F; by combining this 
with the Green's type formula of [HaVa, Lemma 5.2], we deduce that 

dUhiX;m,m) = A"(27r/i)-""^ [ Eh{X; m, OEhiX, m', ^ d^dX. 

Here corresponds to the standard volume form on the sphere S". The assumptions (Al) 
and (A2) are then satisfied. In fact, using [HaVa], one can define distorted plane waves and 
verify assumptions (Al) and (A2) for the more general case of scattering manifolds. 

Outgoing/incoming decomposition. We now construct the decomposition (4.6) of E^ 
into the outgoing and incoming parts and verify assumptions (A3)-(A8). Take xo S C°°{M) 
(thus constant in ^) supported in {x < eq} and equal to 1 near {x < £o/2}, so that assump- 
tions (A3) and (A7) hold, where we put ei := eo/2. We next put 

E'fl{X,^;m) :=e^'"-«, x{m) < so, 
so that (A4) holds with 6" = 1 and (A8) follows. We then claim that 

Eh = XoE'i + El, (6.4) 

where 

El := -Rh{X)Fh. Fh{\0 = (h'^A - X''){xoE'iiX,0) = [/i^A, xo]i?^(A, ?). 

We can apply Rh{X) to Fh{X,^) as the latter lies in C^(M); in fact, suppF/j C {eo/2 < a; < 
eo}- To show (6.4), note that the incoming resolvent Rh{—X) satisfies 

RhWxi = xoRlWxi - Rh{X)[h^A,xo]RlWxi 
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if xi £ C°°(M) is such that XoXi = Xi and R^{X) is the incoming scattering resolvent of the 
free semiclassical Laplacian /i^A on M"+^ (we use again the isometry £ ~ M"+^ \ B{0, R)). 

n iX 

Multiplying the last equation on the right by x~2e'i^, taking the Schwartz kernel of the 
obtained operators and considering the limit in the second variable along a line with tangent 
vector ^ as in (6.2) (using the formula (6.1) with c = 1 for the expression (6.2) for the free 
resolvent R^{X)), we obtain (6.4). 

Microlocalization of Ej^. It remains to verify assumptions (A5) and (A6). By rescaling h 
and using that Efi{X, •) depends only on X/h, we may assume that A = 1. Fix ^ and take 
X2 G C°°(M) equal to 1 near {x < Eq}, but supported inside £. Then 

X2El = Rl{X)F^, Fl := {h'A - X'){x2EI) = Fj, + [h^A,X2]El. (6.5) 

The function is supported inside {x > eo/2} and 

< Ch{l + \\Eh\\L2{{x>eo}))- 

The free resolvent Rl{X) is bounded H~l^^^ Lf^^ with norm 0{h-'^) by [Bu02, Prop 2.1]; 
therefore, for each compact set K C M, there exists a constant Ck such that 

\\EI\\l2(^K) < Ck{1 + \\Eh\\L'2{{x>eo}))- 

This shows (A5), namely that the function 

^1 K 

1 + \\Eh\\L2{{x>eo}) 

is /i-tempered. To prove (A6), we use semiclassical elliptic estimate and propagation of 
singularities (see for example [Vail, Section 4.1]). We have 

{h'A - X')El = -h, h :-- 



1 + \\Eh\\L'2{{x>eo}) 

Now, Fh is a Lagrangian distribution associated to {{m,(^) \ m £ supp((ixo)}j therefore, 

WF;,(F;,) C WFhiFh) C W^, 

with C S*M defined in (4.9). 

Take (m, u) G WF?,(^^). By the elliptic estimate, (m, i^) £ S*M. Next, if 7(i) = g\m, u), 
then by propagation of singularities, either ^{t) G WF/j(F;i) C for some f > or 7(t) G 
WF/i(£'^) for all t > 0. Now, the free resolvent R^{X) is semiclassically incoming in the 
following sense: if / is a compactly supported /i-tempered family of distributions, then for 
each (m',!^') G WF/i(i?^(A)/), there exists t > such that g^{m',i'') G supp/. This can 
be seen for example from the explicit formulas for i?^(A), see [Mc]. By (6.5) and since 
supp(F^) C {x > eo/2}, we see that for {m' ,v') G WF/i(£'^), we cannot have x{m') < eo/2 
and m' ■ u' > 0. Therefore, if j{t) for all t > 0, then 'y{t) is trapped as t ^ +oo; this 
proves (A6). 
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7. Hyperbolic near infinity manifolds 

In this section, we verify the assumptions of Section 4 for certain asymptotically hyperbolic 
manifolds. Let (M, ^f) be an (n+l)-dimensional asymptotically hyperbolic manifold as defined 
in the introduction. It has a compactification M = M U dM and the metric can be written 
in the product form (1.3): 

_ dx'^ + h{x) 
^ x"^ 

where x is a boundary defining function and h{x) a smooth family of metrics on dM defined 
near x = 0. The function x putting the metric in the form (1.3) is not unique, and those 
functions (thus satisfying |(ilog(x)|g = 1 near dM) are called geodesic boundary defining 
functions. The set of such functions parametrizes the conformal class of h{0), as shown 
in [GrLc, Lemma 5.2]. The metric is called even if h{x) is an even function of x, this 
condition is independent of the choice of geodesic boundary defining function. A choice of 
geodesic boundary defining function induces a metric on dM by taking /iq = h{0) = x'^glxaJj^ 
and therefore one has a Riemannian volume form, denoted d^, on dM induced by the choice 
of X. Any other choice x = e^x of boundary defining function induces a volume form 

d^ = e"^ode where cjq = colgM- (7-1) 

We will further assume that M has constant sectional curvature —1 outside of some com- 
pact set, even though some of the assumptions of Section 4 hold for general asymptotically 
hyperbolic manifolds. 

7.1. Geometric assumptions. Let {M,g) be an asymptotically hyperbolic manifold. The 
assumption (Gl) is satisfied. We are now going to prove a Lemma which implies directly 
that the assumptions (G2) and (G3) are satisfied, except that this only proves continuous 
dependence of .^-i-oo in {m,v) in (G3). To prove dependence in a general setting, a bit 
more analysis would be required, but we shall later concentrate only on cases with constant 
curvature near infinity, in which case the dependence is smooth (see below). 

Lemma 7.1. Let {M,g) an asymptotically hyperbolic manifold. Then there exists eo > 
such that the function x satisfies (4.1) and for any unit speed geodesic ^{t) = (m(t), i/(t)) 
with x(m(0)) < eo (^nd dtx{m{t))\t=o < 0, one has the following: dtx{m{t)) < for all t > 
and m{t) converges in the topology of M to some point S,+od £ dM . More precisely, the 
distance with respect to the compactified metric g = x'^g between m{t) and S,+oo is bounded 
by 

d-gim{t),C+^) <Ct-\ 

Proof. Consider coordinates (m, z^) = (x, y; pdx + 9 ■ dy) on T*M near the boundary dM = 
{x = 0}. The geodesic flow is the Hamiltonian flow of p/2, where p{m, v) = x^{p'^ + |^|^^); 
if dots denote time derivatives with respect to the geodesic flow, we get 

x = px^, p = -x~'^p{m,u) - x'^dxh(^^y-){e,0)/2. (7.2) 

Since dxh^^ y-j is smooth up to x = 0, there exists a constant C such that 

\x^dxh(,^y){e,9)/2\ < Cx\x^y-){e,e) < Cp{m,u). 
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Therefore, there exists Eq > such that along any unit speed geodesic, we have 

X < £0 ^ P = -2;"^ + 0(1) < -x~V2 < 0. (7.3) 

This in particular implies (4.1). 

Now, let 7(t) = {x{t),y{t); p{t),6{t)) be a unit speed geodesic and assume that 2;(0) < Eq 
and x{0) < 0. It follows from (4.1) that for t > 0, we have x{t) < and thus x{t) < sq. 
(Indeed, for each s > the minimal value of x{t) on the interval [0, s] has to be achieved at 
t = s.) It remains to show that as t — )• +00, x{t) converges to and y{t) converges to some 
^+00 £ dM. For that, note that by (7.3), p{t) < — ^/2 for t > 0; since x(0) < 0, we have 
p(0) < and thus 



Pit) < -'-^t. 

Setting u{t) := x{t) , we find u{t) = -p{t) > (^0 V2)i; therefore. 

In particular, x{t) — )• as t — ?■ +00. Now the equation for y{s) tells us that 

n 



1 



m = x^h1iy^xek = o{x{t)) = o{t' 



k=\ 

and therefore |y(t) — y(t')| < Cji! for any t > t' > 0, which implies lim(_s.oo ""^(i) = '^00 for 
some ioo G QM and \m{t) - Cool = Oiljt). □ 

The geometric assumption (G4) is a more complicated one, and we will restrict ourselves 
to asymptotically hyperbolic manifolds with constant curvature —1 in a neighbourhood of 
dM and x a geodesic boundary defining function. Let ^ G 9M, then there exists a neigh- 
borhood of C in M, and an isometric diffeomorphism ■0^ from n M into the following 
neighbourhood Vg^ of the north pole go in the unit ball B := {m G ]R"+^; \m\ < 1} equipped 
with the hyperbolic metric g^n+i (the hyperbolic space M^~^'^ is B equipped with g^n+i): 

:={gGB| |g-go| <l/4}, g^„+i = A^^—^-^ (7.4) 

where V'^(0 = I0 ^^^d | • | denotes the Euclidean length. This statement is proved for instance 
in [GuZw, Lemma 3.1]. We shall choose the boundary defining function 

"° = '(TTM) ^'-'^ 

on B. Note that we might not be able to find ifj^ such that x = ip^XQ. This would be possible 
only if the boundary is globally conformally flat and x is chosen so that the metric on dM 
induced by x'^g is flat. 

We define for each p G = dM the Busemann function on B 



= log (^) 
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The geodesic trajectory g^{q, dcl)^{q)) generated by the differential dc/)* converges (in the 
Euchdean ball topology) to p and the Lagrangian manifold 

A«:={(g,#^(g))G5*]H"+i|9e]B} 

is the stable manifold of the geodesic flow associated to p on IB . The level sets of (j)^ are 
horospheres based at p. We cover a neighbourhood of dM by finitely many V^. for some 
G dM and take a partition of unity Xj £ C°°{dM) on dM with Xj supported in V^j DdM. 
Then there exists e > such that for all j and all ^ G supp Xj ; the set 

U^:={meM\ dg{m, ^) < e} (7.6) 

lies inside V^^., where g = x^g is the compactified metric. Put 

:= {{m,i) e M xdM\meU^}. 

Define the function 

M"^) E^^-(^)'^Je,(?)(^5.(H), {rn,0 G U^. (7.7) 
j 

Since ip^^ are isometries, each function i?i>|(m) := i;^^^ (g)(V'Cj(^)) is such that d(l)''^{m) is the 

unit covector which generates the unique geodesic in M starting at m, staying in for 

positive times, and converging to ^ (therefore, the difference of any two functions (/)| for 

different j is a function of ^ only). Therefore dm4'^{rn) = Xj{0'^m(p''^i'm) is also equal to 
this unit covector; (G4) and (G5) follow. The dependence of all objects in m,^ is smooth 
here. Finally, (G6) can be reduced via ip^. to the following statement that can be verified by 
a direct computation: if p G S" and q, q' G B, then dp(j)^{q) = dp(j^{q') if and only if q and q' 
lie on a geodesic converging to p, and the matrix dpg(j)^{q) has rank n. 

7.2. Eisenstein functions and analytic assumptions. Let {M,g) be asymptotically hy- 
perbolic. The Laplacian A on {M,g) has absolutely continuous spectrum on [n?/A,oo) and 
a possibly non-empty finite set of eigenvalues in (0,n^/4). By [MaMe, Gu], if g is an even 
metric'^ the resolvent of the Laplacian 

R{s) := (A — s{n — s))^^ defined in the half plane Re(s) > n/2 

admits a meromorphic continuation to the whole complex plane C, with poles of finite rank 
(i.e. the Laurent expansion at each pole consists of finite rank operators), as a family of 
bounded operators 

R{s):x^L^{M)^x-^L^{M), if Re(s) - n/2 + iV > 0, 

moreover it has no poles on the line Re(s) = ^ except possibly s = ^, as shown by 
Mazzeo [Ma]. Let us fix a geodesic boundary defining function x on M. By [MaMc], the 
resolvent integral kernel R{s] m, m') near the boundary dM has an asymptotic expansion 
given as follows: for any m £ M fixed 

m ^ R{s;m,m')x{m')-' G C°°(M) 

■^There is a simpler proof by Guillope-Zworski [GuZw] when the curvature is constant outside a compact 
set. 
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and similarly for m' G M fixed and m — )• dM. Since we are interested in high frequency 
asymptotics, we will consider the semiclassical rescaled versions 

Rh{X) := h-^R{n/2 + iX/h), 

Notice that the physical region Res > n/2, in which the resolvent is bounded on L^, corre- 
sponds to ImA > 0, which agrees with our convention for Euclidean case. 

Definition 7.2. Let 1/2 < |A| < 2 and h > 0, then Eisenstein functions are the functions 
in C°°(-/Vf X dM) defined for any fixed ^ G dM by the following limit of the resolvent kernel 
at infinity 

Eh{X,C;m) :=^^^ lim x{m')-^/^-''^''Rf,{X;m,m'), 

C{X/h)m>^(, 

The normalisation constant in (7.8) is like the constant in (6.2) so that in B, Eh{X) is a 
horospherical wave as described below in (7.11). For any ^ E dM, we will denote by Eh{X, .^) 
the function m i— t- Efi{X, ^; m), and we observe that they solve (1.4): 

{h\A-n^/A)-X^)Eh{X,O = 0. 

One also has Efi{X, S,; m) = Eh{—X,S,;m) as an easy consequence of Rh{X)* = R{—X) for 
A G M. Prom its definition, Efi{X,^) depends on the choice of the boundary defining function 
X, but considering such a change we easily see from (7.1) that the density on dM 

{AEhiX, 0,Ehi>^, 0)l'Hm) for A G ^'^°^^{M) (7.9) 

is independent of x. 

Let us recall the decomposition of the spectral measure in terms of these functions. By 
Stone's formula, the semiclassical spectral measure is given by 

dUhiX) = -{Rh{X) - Rhi-X)) dX for A G (0, oo) 
vr 

in the sense that F(/i^(A — n^/4)) = F{X'^)dIlh{X) for any bounded function F supported 
in (0,oo). Now we can write (see [Gu]) for any m,m' 

dU„{X; m, m') = j' Eh{X, ^; m)Eh{-X, m')dC dX. (7.10) 

where (27r/i)"|C(A//i)p A" as /i uniformly in A G [1/2,2]. The assumptions (Al) 
and (A2) are then satisfied in the general asymptotically hyperbolic case (without asking the 
constant curvature near infinity). 

Outgoing/incoming decomposition. To check assumptions (A3)-(A8), we give a repre- 
sentation of the Eisenstein functions as sums of the 'outgoing' part E^ and the 'incoming' 
part El. We assume constant curvature near infinity in what follows. The expression for 
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£'^(A) on hyperbolic space IHI"^-'^ viewed as a unit ball B, defined using the boundary defining 
function xq of (7.5), is given by [GuNa, Section 2.2] 

» /I — |o|2\ n/2+iA//i 

^^(^'^''?)= (^-^) gGB (7.11) 

We thus set £'°(A,^; m) to be 

where (/>^ is the Busemann function defined in (7.7). Viewing the neighbourhood as a 
subset of one of the V^^. Vq^ where Vq^ C IB is defined in (7.4), the Laplacian in this 

hyperbolic chart pulls back to Ajjn+i . Since i;^^(m) = 0^^ ^^^('i/'^^.(m))+Cj(^) for some function 
Cj(^) independent of m, we directly have in {U^ is defined in (7.6)) 

(/i2(A-nV4)-A2)£;0(A;^,^) = 0. 

We let xo £ C°°{dM x M) be a function such that xH-, ") is supported in [7^, equal to 1 near 
^ and smooth in . Therefore we obtain 

F^(A,e) := (/i'(A - nV4) - \^)x^EI[\,0 = [h^A,xoK{X,0 (7.13) 

and we claim that 

FhiX,0 G C~(M) and ||x-1f^(A, ^k^CAf) = 0{h) 

uniformly in ^. Indeed, this is an elementary calculation since from (7.11) we see that 
^hi^^O G x'2~^^hC^{M \ {^}) and in geodesic normal coordinates near the boundary 

[A,Xo] = -x^{dlxo) - '^x{da;Xo)xda; + x'^[A^^),xo] +n{xdxXo) - ^T^^h(x)(.dxh{x))x^{da:Xo) 

is a first order operator with coefficients vanishing in a neighbourhood of We thus correct 
the error by the incoming resolvent -R/i(A) by setting 

EhiX,0-=XoEl{\()+Eli^,0, with EUX,0:=-Rh{X)Fh{X,0 (7-14) 

and this makes sense since for A G M, Rh{X) : x"L^(M) — )• x~"L^(M) for any a > and 
Fh £ xL^{M). We claim that 

Proposition 7.3. The function Efi{X, (,) of (7.14) is the Eisenstein function defined in (7.8) 
for a certain boundary defining function x. 

Proof. Let P^{X) be the resolvent of the hyperbolic space (that is, the incoming right inverse 
to /i2(Ahu+i - nV4) - A^) in the ball model and let xi G C°°(5M x M) be such that xi(?, •) 
is supported in and XiXo = Xi- Through the pull-back by ■0^^- (for each j), the operator 
i?^(A) induces an operator R\{X) on V^^; if C V^^. , then we have the resolvent identity 

Rh{X)xi=XoKWxi-Rh{X)[h'^,XoWhWxi (7.15) 

for A € M, the composition makes sense as a map x'^L^ — )• for any a > 0. Let x 

be a boundary defining function, so in V^., one has iI^^.xq = xe^^ for some function Uj G 
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dM C 



Figure 4. Illustration of (A7) for the half-plane model of 111"+-'^: the set of 
points on trajectories converging to ^ E dM^~^^ with x < and x < e is the 
triangle formed by dashed lines, lying 0(e) close to ^. For e small enough, 
this triangle lies inside the lighter shaded region, denoting the set {xo = !}• 



C°°{dMr]V^-). Then multiplying (7.15) by x "/^ on the right, and taking the restriction 
of the Schwartz kernels on M x dM, we have 

E,,{X,0 = xoEli^i) - ii^(A)[/i2A,xo]^^(A,0 

with El{\;m,i) 

~ C{X/li) limm'^g(a;(m') 2 ^ h Rj^(X-^m,m')) a smooth function of m e 
and C{X/h) the constant in (7.8). Note that the Schwartz kernels of i?^ and i?^ are the 
same on the intersection of their domains, therefore does not depend on the choice of 
j. Now, since E^(X;m,^) in (7.11) is the Eisenstein function on IB for the defining function 
Xo, we deduce that in ^7^ C V^^, one has El{X,^;m) = £;0(A, m)e(t+»^)K(«)-^^«)). Here 
™') ~ '^^j,(5)(V'Cj("^))- Since E'^[X;ni,^) does not vanish, this shows that on 
any intersection dM n V^. n V^^ of the cover of dM by the open sets V^^ n dM, we get 
u)j{(,) — Cj(^) = uJkiC) — CkiO and therefore this defines a global smooth function 6 on dM. 
In its definition, Efi{X,£,) only depends on the first jet of x at dM and thus modifying x to 
be xe^, this shows the claim. □ 

It follows that (A3) and (A4) are satisfied, with 6° = et<^e('"). Assumption (A8) is then 
checked by a direct calculation, with the measure on dM corresponding to the choice of 
the function x in Proposition 7.3. 

Assumption (A7) can be reduced, using the isometrics ip^., to the following statement: if 
{q, v) € S*W^~^^ is directly escaping in the forward direction and converging to some p G S", 
then |g — p| < Cxo(g) for some global constant C; the latter statement is verified directly, 
see Figure 4. 

Microlocalization of E^. Finally, assumptions (A5) and (A6) follow, by rescaling h and 
using that Eh{X, ) is a function of X/h, from 

Proposition 7.4. Let Kq C M be a compact set containing a neighborhood of the trapped 
set. Assume that A = 1 and define 

^^^^'^^^ = ttwtM\ — • ^'-''^ 

1 + \\Eh[X,^)\\L2^K^^'j 
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Then: 

1. Ej^(X,S,) is h-tempered in the sense of (3.2). 

2. The wavefront set WFh{Ej^) is contained in S*M. 

3. If (m, u) G S*M and g^{m, v) escapes to infinity as t ^ +oo and never passes through 
the set 

:= {{m,dm(p(,im)) \ m £ supp(5mXo)} 
for t > 0, then (m, u) WFft(^i). 

Moreover, the corresponding estimates are uniform in X £ [1/2, 2] and ^ £ dM . 

Proof. We will use the construction of [Vail]. (See also [ValO]; note however that in that 
paper L+ and L_ switch places compared to the notation of [Vail] that we are using.) Let 
Meven (called Xo, 

even in [^^H]) the space A4 with the smooth structure at the boundary 
dM changed so that is the new boundary defining function. As in [Vail, (3.5)], introduce 
the modified Laplacian 

Pi(A) := x-^x-'{l + x2)'^/4-"/8(/i2A - s(n - s))(l + x^)''/^-'/^x' , s := n/2 + iX/h. 

(The conjugation by (1 + x^)*/^""/*^ is irrelevant in our case, as s/A — n/S = iX/{Ah) is purely 
imaginary. In [Vail], it is needed to show estimates far away in the physical plane, that 
is for Res ^ 1.) Note that we change the sign of A in the conjugation (in the notation 
of [Vail], Pi{X) = Pa with a = —X/h); therefore, our resolvent will be semiclassically 
incoming, instead of semiclassically outgoing, for A > 0. The operator Pi is smooth up to 
the boundary of Meven; as in [Vail, Section 3.5], we embed Meven as an open set in a certain 
compact manifold without boundary X, and extend Pi to a differential operator in ^'^(X). 
We also consider the semiclassical complex absorbing operator Q(A) G ^'^(X) satisfying the 
assumptions of [Vail, Section 3.5]; in particular, Q{X) is supported outside of Meven C X. 
Then (Pi (A) — iQ{X))~^ : C°°(X) — )■ C°°(X) is a meromorphic family of operators in A, and 
for / G C°°(X), we have (see the proof of [Vail, Theorem 5.1]) 

^-(1 + x2)"/8-^/4(Pi(A) - iQ{X))-'f\M = Rh{X){l + x2)«/8-^/Vx2(/lj,,). 

Here RhiX) is the incoming scattering resolvent on M. In principle, depending on the choice 
of (5(A), the operator (Pi(A) — iQ{X))~^ could have a pole at A. However, as R{X) does not 
have a pole for A G [1/2,2], the terms in the Laurent expansion of (Pi(A) — iQ{X))~^ have 
to be supported outside of Meven and we can ignore them in the analysis. 

Let Fji G C°°(X) be any function such that Fh = 0(h) ttn for all N, and 

h 

f, = (i + xY/'-'/V+2(f,1a/). 

Such a function exists as x'^xqE'^ S C°° (Meven \ 0) Xo G C"^ (Meven), and 
Fh = x^+'{l + xY^^-'/^[Pi{s),xo]{l + x^yl^-'^'^x-'El 
is supported away from ^. Define the function E\ G C°°(X) by 

gi_ (Pi(A)-.Q(A))-iF, 

^ l+||i^/.(A,OllL^(i.o)' 
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Figure 5. Left: physical space picture of geodesies converging to ^. The 
darker shaded region is the support of dxo, and thus of Fh- In the hghter 
shaded region, xo = 1- Right: phase space picture near ^ after the conjugation 
of [Vail]. L_ is the sink consisting of radial points, Q is the complex absorbing 
operator, and the shaded region corresponds to the wavefront set of Fh- The 
vertical line hitting L_ is the boundary of Meven, while the horizontal line is 
the fiber infinity. In both pictures, we mark two points {m, v) satisfying the 
assumption of part 3 of Proposition 7.4 and the forward geodesies starting at 
these points. 

Then 

Consider the map l : T*M ^ T*X given by 

L{m, v) = ^m, V — Inx(m-) — - ln(l + x{m)'^)^ ^ , m G M, v G T^M; 

then for an /i-tempered u € C°°(X), 

WF;,(x^(l + x^r/'-'/\\M) = rHWFhiu)). 

Then 

WF/,((Pi(A) - iQ{\))El) n T*M C i{W¥h{Fh)) C i{W^). (7.17) 
Now, as \\EI\\l2^^Ko) < C and thus \\EI\\l2(^Ko) <C + ||^fe||L2(^g), we have 

Consider an operator Qk G xI/'=°™p(X) supported in such that (t{Qk) < everywhere 
and each unit speed geodesic ^{t) either escapes as t — )• +<x or passes through the region 
{(t{Qk) < 0} at some positive time. This is possible since Kq contains a neighborhood of the 
trapped set. Then the operator Pi (A) — iQ{\) — iQk satisfies the semiclassical nontrapping 
assumptions [Vail, Section 3.5]; therefore, by the nontrapping estimate [Vail, Theorem 4.8], 

WeIWl^^x) < Ch-'\\{P,iX) - iQ{X) - iQK)Elh2^x) 

< Ch^-^\\Fh\\L2(x) + Ch~^\\QKEl\\L2(^xy 

However, is bounded by ||£'^||^2(^g); therefore, ||-E^||/^2(x) = 0{h^^) and in particu- 

lar Ej^ is tempered; it follows that Ej^ is also tempered. This proves part 1 of the proposition; 
part 2 follows by ellipticity (note that WFh{Fh) C C S*M). 
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Now, assume that (m, z^) G S*M satisfies the assumption of part 3 of this proposition. 
Then it fohows directly from (7.17), the analysis of [Vail, Section 2.2], and the definition 
of i, that the Hamiltonian flow line of (t{Pi) starting at i{m,v) converges to the set L_ of 
radial points as t — )• +oo and does not intersect WF/i((Pi(A) — iQ{\))E^) for t > 0. In a 
fashion similar to the global argument of [Vail, Section 4.4], we combine elliptic regularity 
and propagation of singularities (see [Vail, Section 4.1]) with the radial points lemma [Vail, 
Proposition 4.5] for L_, to get i(m, v) WF/i(£'^). Therefore, (m, v) WF/j(£'^) as required. 

□ 

Appendix A. Limiting measures for hyperbolic quotients 

In this appendix, we give an explicit description of the limiting measures //^ in case when 
M is a hyperbolic quotient r\]HI"'^-'^, in terms of the group T. This is a particular case of 
asymptotically hyperbolic manifolds discussed in Section 7. 

A.l. Convex co-compact groups. Let B be the unit ball in R"+-^, and W^~^^ the (n + 1)- 
dimensional hyperbolic space, which we view as B equipped with the constant negative 
curvature metric g^n+i := 4|dmp/(l — |mp)^. The boundary §" = (9B is the sphere of 
radius 1, which is also the conformal boundary of H""*"^. A convex co-compact group T of 
isometrics of H""*"^ is a discrete group of hyperbolic transformations (i.e., transformations 
having 2 disjoint fixed points on B) with a compact convex core, and F is not co-compact. 
The convex core is the smallest convex subset in F\EI'^+\ which can be obtained as follows. 
The limit set Ar of the group and the discontinuity set Or are defined by 

Ar := {7(m) G B;7 G F}nS", J^r := S" \ Ap (A.l) 

where the closure is taken in the closed unit ball B and m G B is any point (the set Ar does 
not depend on the choice of m) . The group F acts on the convex hull of Ap (with respect to 
hyperbolic geodesies) and the convex core is the quotient space. 

An important quantity is the Hausdorff dimension of Ar 

6 := dimn Ar < n (A. 2) 

which in turn is, by Patterson [Pa] and Sullivan [Su79], the exponent of convergence of 
Poincare series: for any m G B, 

we henceforth denote by d{-, •) the distance function of the hyperbolic metric on B. Notice 
that the series (A. 3) is locally uniformly bounded in m G B. 

The group F acts properly discontinuously on fir as conformal transformations of the 
sphere and the quotient space F\ilr is a smooth compact manifold of dimension n. The 
quotient 

M = F\H"+i 

is a smooth non-compact manifold equipped with the hyperbolic metric g induced by g^^n+i , 
and it admits a smooth compactification by setting M = MU (Pyj^r), i-e. with dM = F\r2r- 
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Then M is an asymptotically hyperbolic manifold in the sense of Section 7, of constant 
sectional curvature —1. We shall denote the covering map by 

TT : IB U fir M. 

We refer the reader to [Ni] for more details and properties of convex co-compact groups. 

A. 2. Limiting measures in this setting. In constant curvature, it turns out that the 
limiting measure fi^ exists for all ^ (rather than for Lebesgue almost every ^ as in Section 4.3), 
and can be described as a converging sum over the group. We give an expression below, which 
is the same as the one obtained in [GuNa] when 6 < n/2. 

For ^ G S", we let 0^ be the Busemann function^ on the unit ball B defined by 
The map ^ defined by 

gives a diffeomorphism between the unit cotangent bundle S*M'^'^'^ and B x S", and satisfies 




if dfiL is the Liouville measure (viewed as a volume form on the unit cotangent bundle) and 
the canonical measure on S". (This is a more general version of (A8) for the considered 
case.) Any isometry 7 of IH"^"-'^ acts on both spaces by 

j.{m,u) = (7m, {d'y{m)u*y), for {m,u) £ S'*]H"+^ 
7.(m,0 = (7m,7^), for (m, ^ € B x S", 

where * denotes the map identifying T*I¥'+^ with r]H"+^ through the metric. We have 
^>(7.(m, ^)) = 7.$(m, and thus $ descends to a map r\(H"+i x S") S*{T\W+'^), which 
we also denote by 

The limiting measure in the considered case is given by 

Lemma A.l. Let M = r\]H"+i be a quotient o/M^+i by a convex co- compact group T of 
isometrics, let J- be a fundamental domain. Then the measure M7r(^) of (4.12) exists for all 
G r^r O'lT'd is described as a converging series by the following expression: if ^ £ JlrnJ-" and 
a e C^{S*M), then 

[ a dfi^i^^^ = f Y] a{m, d(/>^5(m))e"('^^«('")+i°s dvole^+i (m) 

where (j)^{m) is the Busemann function onB associated to ^ G S" and |d7(C)| is the Euclidean 
norm of d'y{(,). 

^In Section 7, we used the coordinate g £ B, p G §" for certain charts near infinity of M, and the notation 
(t>^(q) for the Busemann function on the ball. This was to avoid confusion with the coordinate m, ^ on M, dM. 
We keep in this appendix the notation ^^(m) to match the notation of the general setting of the article. 
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Proof. We can view a as a compactly supported function on the unit cotangent bundle S*J- 
over a fundamental domain J-" C B and we extend a by in S'*E["+-'^\5*J-" (the resulting 
function might not be smooth, but it does not matter here). The flow on S*M is obtained 
by projecting down the geodesic flow of the cover S*M^~^^. Let G Or H J-, then small 
neighbourhoods of 7r(^) in M are isometric through vr to small neighbourhoods of in the 
unit ball B. By the construction of the decomposition (4.6) for the asymptotically hyperbolic 
case in Section 7.2, the function £'J|(A, 7r(^); 7r(m)) is equal to e^'^/'^+^^/^)'t>i('n^) for m near ^ 
(C being fixed) and thus = e"'^e("^). One has 

a((7-*(m,d(/)5(m)))e"<^«('")dvolM(H = / 5(5-*«>(m, 0)e"'^«("') dvolen+i (m) 

M JT 

where 'a{m,v) := ^^gp a(7.(w-, z/)) is the lift to S*W"'^'^ of the function a on S* M and 
dvoljjn+i (m) is the Riemannian measure on H""*"^. Using the map $ of (A. 4), one can define 
a map g| : B — > B by 

ff*$(m,e) = ^{gl{m),i), 
this is a diffeomorphism which preserves the measure e^'t'e.ivT-) dvoljjn+i . By [GuNa, Lemma 4], 
we have e^'t'd^'^rn) ^ e""^^«('") 1^7(01", but we also have 7.$(m,0 = ^>(7m,7^). Let U+ be 
defined in (G4) and put U := {m \ (m, 7r(^)) G U^}, then U lies in a small neighborhood of 
7r(^) in M. We can identify U with a small neighbourhood U oi in and we get for 
defined in (4.10), 



d/I,(5)= [ a(g-*(m,#^(m)))e'^'^«('")dvolM( 
= / ^a(7.r*^("i,0)e"^^^"^^dvoW 



+i(,m 

r (A.5) 

/ a(<i>(755-*m,7e))e«'^«(-)dvolHn+i(m) 



7Gr 



u 



V / _ a(m,#.^5(m))e"'^s(^ '"'^ dvoWi (m). 



We now observe that for all 7 G F, limf_s.+oo 1 — tfj = 1, since [/ is a neighbourhood of ^ 
in B containing all points directly escaping to ^. This achieves the proof by recalling the 
definition (4.12) of /U^(^) and taking the limit in (A.5) and using the dominated convergence 
theorem, as there exists C,C' > such that for all m in the compact set supp(a) 

Y^^nM^-^m) ^Y^ n - '^"_'^h " < c sup e-"'^^^-^-'") < C 

by locally uniform (in m) convergence of Poincare series (A. 3) at s = n. □ 

Appendix B. The escape rate 



Let us discuss the classical escape rate in some particular cases, following the work of 
Bowen-Ruelle [BoRu], Young [Yo], and Kifer [Ki]. 
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B.l. Escape rate and the pressure of the unstable Jacobian. We consider {M,g) a 
complete non-compact Riemannian manifold which has a compact set Kq C S*M which is 
geodesically convex, that is any geodesic trajectory in S*M which leaves Kq never comes 
back: 

3*0 < h, g^°{m, u) € Kq and 5*H"^> v) £ M\KQ^^t> h, g\m, u) £ M\ Kq. 

We assume that Kq contains a neighborhood of the trapped set K. The examples we consider 
are {M,g) which are hyperbolic or Euclidean near infinity, and Kq = S*M n {x > Eq} with 
x,£q given in (G2). The trapped set from Definition 4.1 can be written as 

K=f]g\KQ) = f]g^{KQ) 

This is a compact maximal invariant set for the flow g*. We deflne the escape rate as in 
[Yo, Ki] by 

(5 : = hm sup - log /i/, (T(t ) ) , 
s>oo * 

with the Liouville measure and T{t) defined in (1.9). Note that, since Kq is geodesically 
convex, we have T{t2) C T(ti) for < ti < t2- The escape rate is clearly non-positive. 

In this section, we assume that fiiiK) = and write Q in terms of the topological 
pressure of the flow, under certain dynamical assumptions. More precisely, we assume that 
the trapped set K is uniformly partially hyperbolic, in the following sense: there exists £/ > 
and a splitting of T(S*M) over K into continuous subbundles invariant under the flow 

T^S*M = Ef e G K 

such that the dimensions of i?" and E^^ are constant on K and for all e > 0, there is to £ 
such that 

Vw E E'^, \dglv\ > e^/>|, 
\/v £ Ef, \dg\v\ < e^*|i;|. 
Let J'" be the unstable Jacobian of the flow, deflned by 

J«(z) ■=-dt{deidg\\E^)\t=o 

where dg^ : E"^ — )■ E'^t and the determinant is deflned using the Sasaki metric for choosing 
orthornormal bases in E"^. If ^ is a (7*-invariant measure on K, one has 



Vz £ K, \ft> tQ, 



I rdf^ = - [ ^A+d/z 



where A^l-z) are the positive Lyapunov exponents at a regular point z £ K counted with 
multiplicity (regular points are points where the exponents are well deflned, and this is 
set of full //-measure by the Oseledec theorem). It is also direct to see that fj^J'^dfi = 
- log det ( dfif ^ I ) d/i . 

The topological pressure of a continuous function 99 : X — t- M with respect to the flow can 
be deflned by the variational formula 

P{^):= sup (h^{g')+ [ ^dfi) (B.l) 
fi€M{K) V JK J 
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where Ai{K) is the set of (7*-invariant Borel probabihty measures and h^{g^) is the measure 
theoretic entropy of the flow at time 1 with respect to fi. In particular P(0) is the topological 
entropy of the flow. 

A particular case of uniformly partially hyperbolic dynamics is when K is uniformly hy- 
perbolic, that is when there is a continuous splitting i?^* = MHp © into flow direction (Hp 
is the vector field generating the geodesic flow) and stable directions E'^ where for t >tQ 

G E'^, \dgiv\ < e-^f^\v\. 

The set K is called a basic hyperbolic set and the flow is said to be Axiom A when the 
periodic orbits of on K are dense in K and g^lx is topologically transitive. 

It is proved by Young [Yo, Theorem 4] that if K is uniformly partially hyperbolic, then 

Q = hm jlogfiUm) = Pin- (B.2) 

In the Axiom A case, the same formula was essentially contained in the work of Bowen-Ruelle 
(using the volume lemma [BoRu, Lemma 4.2 and 4.3]). Moreover by [BoRu, Theorem 5], 
if the incoming tail r_ (which is the union of stable manifolds over the trapped set) has 
Liouville measure 0, then P( J") < 0. Thus we deduce by (5.2) 

fXL{K) = and 5* is Axiom A =^ P{r) < 0. 

Young [Yo, Theorem 4] gives a lower bound Q > P{—^j^J) which applies without any 
assumption on K (but we are more interested in an upper bound). 

B.2. Relation with fractal dimensions in particular cases. Assume first that the met- 
ric has constant curvature —1 in a small neighbourhood of the trapped set K (this includes 
the case of convex co-compact hyperbolic quotients studied in Appendix A). Then the ge- 
odesic fiow on S*AI is uniformly hyperbolic on K and has Lyapunov exponents (with 
multiplicity 1) and ±1 (each with multiplicity n). Therefore, the maximal expansion rate 
Amax from (1.11) is equal to 1, one has J^{z) = —n for all z £ K, and (see for example [Fa, 
Theorem 4]) 

P(J") = htop{K) -n = (dimi^(K) - l)/2 - n (B.3) 
where /itop is the topological entropy of the flow on K, and divuHiK) € (0,n) is the Haus- 
dorff dimension of K (which is equal to the Minkowski box dimension in this case). For 
convex co-compact hyperbolic quotients r\]HI"'"'"^ (see Section A for definition), one has by 
Sullivan [Su84] 

5 := dim^^(Ar) = htop{K) (B.4) 
where Ap is the limit set of the group F defined in (A.l). 

If g has negative pinched curvature near the trapped set, then one still has upper and 
lower bounds on P(J") in terms of htop{K) and the pinching constant. If the trapped set 
K is uniformly hyperbolic, it is also shown in [Fa] that dimffiK) < 1 + 2/itop(-f^)/Amax- 
In dimension 2 there is an explicit relation between the Hausdorff dimension dimH{K) and 
entropies of certain measures for Axiom A cases: if 

1 1 

a"(2;) = lim - log ||ci5r*|£;« 11 > 0, a^'iz) = lim - log < 

t-5>0 t t-5>0 t 
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then Pesin-Sadovskaya [PcSa] show the following formula 

dimniK) = 1 + T + with P(-ra") = P{-t'a') = 0. 

Appendix C. Egorov's theorem until Ehrenfest time 

In this section, we prove Proposition 3.9, following the methods of [BoRo], [AnNo, Sec- 
tion 5.2], and [Zw, Theorem 11.12]. Without lack of generality, we assume that to > 0. 

C.l. Estimating higher derivatives of the flow. First of all, we need to estimate the 
derivatives of symbols under propagation for long times. Consider the open set 

Ui = {{m, u) G T*M I m G [/, I - 2£e < \v\g < I + 2ee}. 

For each k, we fix a norm || • ||cfe({7j^) for the space C^{Ui) of k times differentiable functions 
on C/i. (The particular choice of the norm does not matter, as long as it does not depend on 
t.) The following estimate is an analogue of [AnNo, (5.6)]; we include the proof for the case 
of manifolds for the reader's convenience. 

Lemma C.l. Take Ai > (1 + 2ee)A 

max- Then for each k, there exists a constant C{k) such 
that for each a G C^{Ui) and each t G M, 

\\aog%,f^jr,)<C{k)e'^^'^'^\\a\\c.(^jr,y (CI) 

Proof. Without loss of generality, we assume that t > 0. We first recall the formula for 
derivatives of the composition b o ip of a function b G C°°(M°') with a map ■0 : — >■ W^: 

m 

d%bo^) = Y,^c.Adn-J^b) o i; -lld^'i^j,, (C.2) 

a,j 1=1 

where Ca,j are constants, ji, . . . ,jm G {1, . . . , d}, and ai, . . . , Om are nonzero multiindices 
whose sum equals a. We see from (C.2) that (C.l) is implied by the following estimate on 
the derivatives of the fiow (required to hold in any coordinate system): 

\a\<k^ sup < C7„el°l^^*. (C.3) 

Uing-HUi) 

The converse is also true, which can be seen by substituting cooordinate functions in place 
of a in (C.l). 

To estimate higher derivatives of the fiow, we will need several definitions from differential 
geometry. For a vector field X on Ui, define its pushforward glX by 

X{a o g') = {{giX)a) o g\ a G C°° (5* (^^) ) • 

Then 5* A is a vector field on g (Ui). In local coordinates, we have 

{gixy =J2iX'dig'^)og-'. 
I 

Note that since 5* = exp{tHp/2), where Hp is the Hamiltonian vector field of p, and since 
glHp = Hp, we have 

dtgiX = -^[Hp,gtX] = -^gl[Hp,X]. (C.4) 
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We fix a symmetric affine connection V on T*M. For vector fields X and Y, consider the 
differential operator V^y, acting on functions or on vector fields, defined as follows: for a 
function / and a vector field Z, 

VjcYf = XYf - {VxY)f, V\yZ = Vx^yZ - Vy^yZ. (C.5) 

In local coordinates, we have (using Einstein's summation convention) 

V\yf = X^Yi{dlf-T\^dif), 

{v\yzr = x'Y^idfjZ"" + r™ + r^a.z" - rf^a„z"^ 

Here r\j are the Christoffel symbols of the connection V. The advantage of V^y over XY 
is that the coefficients of this differential operator at any point depend (bilinearly) only on 
the values of X and Y at this point, but not on their derivatives. 

We now return to the proof of (C.l). The estimate (C.3) for k = 1 follows directly from 
the definition (3.17) of Amax- It is then enough to assume that (C.3) holds for some k > 1 
and prove the estimate (C.l) for k + 1. It suffices to show that for any two vector fields X, Y 
on T*M and any a G Cq°{Ui), we have the estimate 

||Xy(ao5*)||^,_,(j^) < Ce(^'+i)^^*||a||^.^.i(^). (C.6) 
The left-hand side of (C.6) is equal to \\{glXglYa) o g'^\\f-,k-i(ij^y We first claim that 

Indeed, in local coordinates 

^^llXgiY'^) ° 9' = {X^d^gi){Y^d(sg)){{dfja - r[^dia) o g'). (C.8) 

We can now apply (C.2) to get an expression for any derivative of order no more than k — 1 
of (C.8). The result will involve derivatives of orders l,...,k of (7*, but not its k + I'st 
derivative; therefore, we can apply (C.3) to get (C.7). 

Given (C.7) and (C.5), it is enough to show 

ll((V,*x5:n«) °9lc^-i(I7r) < Ce('=+i)^iia||^.(j;^). (C.9) 

The vector field Vgtx5*i^ involves the second derivatives of 5*, therefore the left-hand side 
of (C.9) depends on the k + I'st derivatives of g^ and we cannot apply (C.3) directly. We 
will instead use the method of the proof of [BoRo, Lemma 2.2], computing by (C.4) 

-^[H„gix]9iY - Vgtx[Hp,giY]) = ^g^'Zt, 
where Zt is the vector field given by 

Zt = VlxgiYHp + RviHr„gtX)igiY). 
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Here R\/ is the curvature tensor of the connection V. Then 

^glxgtY = gtiVxY) + gi-'Zs ds. (C.IO) 

We have 

It is then enough to handle the integral part of (C.IO). The field Zg depends quadratically 
on the first derivatives of g^, but does not depend on its higher derivatives; therefore, writing 
an expression for Zg in local coordinates similar to (C.8), we get for a G C^{Ui), 

Applying (C.l) for the norm (given by the induction hypothesis) and using the geodesic 
convexity of U, we get 

\\{{gl 'Zs)a)og*\\^^_,,jj-.ds= / \\{Zs{aogt ')) o g'\\^,^_,,jj-. ds 
JO 

< C r e(^'+i)^i1|a o «*-1U,,rr^ ds < C t e(^-+i)^i^e'=^i(*-^) lla|Uw7r^ ds 



Jo 

S oe l|a|lcfc(C/i) 

and the proof is finished. □ 

C.2. Proof of Proposition 3.9. The proof of Proposition 3.9 is based on repeatedly ap- 
plying the following corollary of Lemma C.l. The functions b^^^ below will be the remainders 
in the formula for the commutator [h? IS., A^^\t)\, while the functions c*--'^ will be the errors 
arising from multiplying our operators by Xi and X2. 

Proposition C.2. Take Ai > (1 + 2ee)Amax- Fix to > and let ip G C^(C/i) satisfy 
\ip\ < 1. Assume that e C°°(r*M) and for each j > 0, b^^\t) G C°°([0,to] x T*M), and 
c^^^ £ C°° (T* M) , with support contained in some j -independent compact set. For j > 0, 
define a^^^ G C°°([0,to] ^ T*M) inductively as the solutions to the equations 

a(0)(0) = ao, a(^'+i)(0) = <f ■ a^^\to) + c'^^+^^; 
dta^^\t) = ]^Hpa^^\t) + b'^^\t). 

Then for each k, and each j , we have (bearing in mind that each a^^^ is supported inside 
some j -independent compact set and thus its C*^ norm is well-defined up to a constant) 

sup \\a^'\t)\\ck^T'M) < Cik){e^'''^'^\\ao\\c. 
te[o,to] 

+ maxe(^-^)^^i*o( sup \\b^'\t)\\ck + \\c^''>\\ck)) , 
0<«<J te[o,to] 

where C{k) is a constant independent of j. 
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Proof. We can write 

a(j\t) = a(^)(0) 05*+ r 60) (s) o 5*-^ ds. 
Jo 

Since Iq is fixed, it is enough to estimate the derivatives of a^^\0). Define 

0<m<i 

applying the Leibniz rule to (p^^\ estimating each nontrivial derivative of(fog"^^o Lemma C.l, 
using that Iv'l < 1 and absorbing the (polynomial in /) number of different terms in the Leib- 
niz formula into the exponential by increasing Ai slightly, we get Hv'^-'^ H^fc = ©(e-^'^^^*"). 
Now, 

a(^')(0) = (^(^'^ • (ao o 5J'*o) + V (y^^^-*) / b^'\s) o gO-0*o-« 
3 

i=l 

Here we put 99(0) = 1. We can now apply Lemma C.l again to get the required estimate. □ 



We are now ready to prove Proposition 3.9. Fix a quantization procedure Op^ on M; our 
symbols will be supported in a certain compact set (in fact, no more than distance to to 
the set U) and we require that the corresponding operators be compactly supported. Put 
Ai = A'o. 

Let / satisfy (3.18). We will construct the operators 

A^i\t) = Ovh( . < t < to, < j < m > 0, 

^0<m'<m ^ 

(7) 

Here the symbols am will be supported in a fixed compact subset of T*M and satisfy the 
derivative bounds 

sup ||a(^)(t)||t;fc < C(A;,m)/i(i"2p,)m-p,fc^ 

te[o,io] 

with the constants C{k, m) independent on j and pj defined by (3.19). The operators Am {t) 
will satisfy the relations 

AW(o) = ^ + 0(/i°°)^_, 
Ag+^\Q) = X2AW(to)Xi + Ovh{c^I^) + 0{hn^-^, (C.12) 
hDtA^l\t) = \[h^/\,A^l\t)] + ^Op,(6(^)(t)) + 0(/i-)^-oo, 

where the symbols bm (t) and Cm are supported in some fixed compact set and satisfy bounds 

sup \\h^l\t)\\c.MiHc^ < C(A;,m)/^(l-2^^■)(™+l)-^^■^ (C.13) 

te[o,to] 

with the constants C{k,m) again independent on j. 
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We construct the symbols am iteratively, by requiring that they solve the equations 

dtali\t)=^-H,a'ii\t)-b^£,{t). 

Here A = Op!,{ao)+0{h°°)^-oo and we put 6^^] = cL^] = 0. The function e C^{Ui) is equal 
to a{Xi)a{X2)tp{\iy\), where V e C^(l - 2ee, 1 + 2ee) is such that = 1 near WFh{A). 

We use the fact that the function is invariant under the geodesic flow. The estimate (C.ll) 
follows immediately from (C.13) and Proposition C.2. As for the equations (C.12) and the 
bounds (C.13), they follow from (C.ll) and the following commutator formula: 

[h^A, Opja)] = ^ Op^iHpa) + OpJ6) + 0{h^h-.., b = Oih^-^qa\\s,)s,, 

true for any p < 1/2 and any a G 5™™^. 
Now, consider the asymptotic sums 

a(^-)(t)~^a(^)(t) 

m>0 

and define the operators A^^\t) = Opf^{a^^\t)). By (C.12), these operators satisfy 
^(o)(0) =A + 0(/i~)^-oo, A(^'+i)(0) = X2A^'\to)Xi + 0(/i°°)^-oo, 
hDtA'^^^t) = ^[/i2A,A(j')(t)] + 0(/i°°)*-oo. 

We then have 

iX2Uito)yA{U{-to)Xiy = ^«(0) + 0{h^)L2^L2. 

It remains to recall that a(')(0) G 5'^°™'' uniformly in /. The principal symbol and microlocal 
vanishing statements follow directly from the procedure we used to construct the symbols 



Appendix D. Proof of quantum ergodicity in the semiclassical setting 

In this section, we illustrate how our methods yield a proof of the following integrated 
quantum ergodicity statement in the semiclassical setting: 

Theorem 5. Let {M,g) be a compact Riemannian manifold of dimension d and assume 
that the geodesic flow on M is ergodic with respect to the Liouville measure fiL on the unit 
cotangent bundle S*M . For each h > 0, let (ej)jgf^ be an orthonormal basis of eigenf unctions 
of h'^A with eigenvalues X'j. Then for each semiclassical pseudodifferential operator A G 
^'°(M), we have 



h 



d-l 



E 



Xje[i,i+h] 



{Aej,ej)L2^M) 



liL{S*M) 



S*M 



a (A) dfii 



as /i — > 0. 



(D.l) 
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A more general version of Theorem 5 was proved in [HeMaRo] , in particular relying on the 
result of [DuGu, PeRo] on o(/i) remainders for the Weyl law when the closed geodesies form a 
set of measure zero. The purpose of this Appendix is to provide a shorter proof. Theorem 5 
is formulated here for the semiclassical Laplacian for simplicity of notation, but it applies 
to any self-adjoint semiclassical pseudodifferential operator P{h) with compact resolvent on 
a compact manifold, if the Hamiltonian flow of the principal symbol p of P{h) has no fixed 
points and is ergodic on the energy surface p~^{0) and we take eigenvalues in the interval 
[0,h]. 

The key component of our proof is the following estimate: 
Lemma D.l. Let M be as in Theorem 5. Then for each A G ^^{M), we have 

h''-' J2 Ue.Wh^M) < {C\\a{A)U2^s'M) + 0{h)f. (D.2) 

Here \\o'{A)\\i2(^g*^,j-^ is the norm of the restriction of a {A) to S*M with respect to the 
Liouville measure. The constant in 0{h) depends on A, hut the constant C does not. 

Proof. Assume first that A is compactly microlocalized. We can rewrite the left-hand side 
of (D.2) as the square of the Hilbert-Schmidt norm of /i^'^"^)/^^]!^! where n[x^i+/i] = 
l[i^(i+/i)2](/i^A) is a spectral projector. It can then be estimated using the local theory of 
semiclassical Fourier integral operators, by (3.24) (applied to the adjoint of the operator in 
interest). 

To handle the case of a general A., it remains to note that if WYh{A) r\ S*M = 0, then 
the left-hand side of (D.2) is 0{h'^), as each Aej is 0{h°°) by the elliptic estimate (Propo- 
sition 3.2; see also the proof of Proposition 4.5). □ 

Putting A equal to the identity in (D.2), we get the following upper Weyl bound: 

#{j \ XjG [1,1 + h]}<Ch^~''. (D.3) 

We can now prove Theorem 5. Take A G ^"^(M); by subtracting a multiple of the identity 
operator and applying the ellipticity estimate, we may assume that A is compactly microlo- 
calized and 

/ a{A)diiL = 0. (D.4) 

Js*M 

Define the quantum average 

{A)T = ^£u{t)AU{-t)dt. 

Here U{t) = e**''^/^ is the semiclassical Schrodinger propagator. By Egorov's theorem 
(Proposition 3.8), for any fixed T the operator {A)t lies in modulo an 0{h°°)i2^i^2 
remainder, and its principal symbol is 

a{{A)T) = {a{A))T := ^ £ <t{A) o / dt. 
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Note that for each j, we have U{t)ej = e**'^^/(^'') and thus {{A)Tej,ej) 
Cauchy-Schwarz inequahty in j and the bounds (D.2) and (D.3), we get 



{Aej,ej). Using 




1/2 



< h' 



d-i 




<C\\{a{A))T\\LHs*M) + OT{h). 



However, by (D.4) and the von Neumann ergodic theorem [Zw, Theorem 15.1], we have 
\\{a{A))T\\L^(S*M) — as T — )• oo. Therefore, for each e > we can choose T large enough 
so that the left-hand side of (D.l) is bounded by e/2 + 0{h). Then for h small enough, it is 
bounded by e; since the latter was chosen arbitrarily small, we get (D.l). 
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